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ABSTRACT 


A  method  for  visualizing  objects  immersed  in  water  is  formulated 
analytically  and  demonstrated  experimentally.  The  technique,  called  "acoustic 
holography,"  is  an  adaptation  of  Gabor's  two-step  imaging  process  known  as 
wavefront  reconstruction  or  holography.  The  hologram  is  first  formed  from 
coherent  acoustic  radiation  and  then  the  image  is  reconstructed  optically 
using  coherent  light  source.  Acoustic  holography  has  advantage  over  other 
schemes  for  imaging  ,in  optically  opaque  media  in  that  lenses  or  other  focusing 
devices  are  not  required,  and  a  complete  amplitude  and  phase  reconstruction 
of  the  scattered  field  may  be  obtained.  Since  instantaneous  amplitude  is  an 
acoustic  observable,  the  reference  field  may  be  simulated  electronically. 
Moreover,  by  resorting  to  heterodyne  or  phase  detection  the  cross  product 
term  between  object  and  reference  signals  may  be  generated  without  the  un¬ 
desired  extraneous  terms  which  occur  in  conventional  holography.  A  scanning 
technique  for  generating  acoustic  holograms  of  underwater  objects  in  the 
laboratory  is  described  in  detail.  Using  this  system,  acoustic  holograms  have 
been  recorded  which  show  angular  resolution  of  3.  6  milliradians,  approximately 
1. 5  times  the  Rayleigh  limit.  A  variable  contrast  television  display  was  used 
to  view  the  acoustic  holograms.  To  limit  the  attenuation  of  acoustic  wave  in 
sea  water  to  a  tolerable  value,  only  acoustic  signal  with  frequencies  below 
1  MHz  should  be  used.  We  used  1  MHz  signal  for  our  experiment.  Con¬ 
sequently,  the  quality  of  any  acoustic  image  is  degraded  by  poor  resolution  and 
specular  reflection.  It  is  suggested  that  diffuse  or  incoherent  illumination 
be  used  to  overcome  the  defect  that  acoustic  images  often  appear  as  diffracted 
f:  highlights  rather  than  as  extended  forms.  Thus,  it  appears  that  acoustic 

holography  is  not  necessarily  optimum  for  all  acoustic  imaging  situations. 
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CHAPTER  I 


INTRODUCTION 


Historically,  optics  lias  referred  primarily  to  the  study  of 
light;  imaging,  to  the  specific  class  of  optical  phenomena  in  which 
material  bodies  (objects)  are  represented  by  geometrically  similar 
patterns  of  light  (images).  However,  this  is  unduly  restrictive  since 
imaging  is  merely  a  consequence  of  the  wave -like  properties  of  electro¬ 
magnetic  radiation.  More  properly,  imaging  concepts  should  be  extended 
to  include  mathematically  related  phenomena  in  which  the  object  is 
"illuminated1'  with  wave  fields  other  than  light  and  invisible  "image" 
patterns  are  formed  from  the  scattered  radiation.  The  scope  of  optics 
is  thereby  enlarged  to  encompass  a  broad  group  of  physically  distinct 
but  mathematically  similar  wave  phenosena. 

Here  we  are  concerned  with  acoustic  imaging,  a  process  belonging 
in  this  more  general  category.  Ir  acoustic  imaging  the  object  scene  is 
"ii? \sinated”  with  acoustic  radiation  rather  than  light,  and  acoustic 
”i~age"  patterns  are  formed  free  the  scattered  acoustic  fields  by 
various  means.  The  acoustic  images  say  then  be  visualized  -  rendered 
visible  -  through  application  of  some  sort  of  acousto-optic  interaction. 
Since  by  definition  all  wave  fields  obey  mathematically  similar  equa¬ 
tions  cf  action  we  are  guaranteed  that,  at  least  to  first  order,  visual¬ 
ized  acoustic  images  will  bear  snse  recognizable  resesblance  to  tbs 
original  seers.  Imaging  of  this  type  may  be  extremely  useful  whenever 
the  object  is  imbedded  in  rc-dis  opaque  to  light. 
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Acoustic  imaging  differs  from  and  should  not  be  confused  with 
the  conventional  echo  ranging  techniques  conroonly  known  by  the  generic 
term  sonar  in  that  imaging  provides  displays  of  the  object  scene  lying 
in  planes  transverse  to,  as  well  as  along,  the  direction  of  propagation 
whereas  sonar  provides  primarily  range  profiles  along  the  direction  of 
propagation;  i.e.,  longitudinal  “images." 

The  reasons  for  the  study  of  acoustic  imaging  should  be  clear. 

The  quantity  of  information  transmitted  via  wave  fields  and  displayed 
as  images  is  enormous.  Indeed,  under  ideal  conditions  (infinitely 
large  appertures  and  unlimited  detector  resolution)  the  maxiausa  infor¬ 
mation  density  in  an  optical  display  is  limited  only  by  the  information 
carrier,  here  scalar  wave  fields  of  wavelength  >.  frca  which  the  image  is 
formed.  Gabor*  and  others  have  shewn  that  wave  fields  cannot  transmit 
information  abort  boundary  details  with  characteristic  spacings  such 
smaller  than  X  because  such  information  is  carried  by  evanescent  waves 
which  decay  so  rapidly  that  at  distances  beyond  a  few  wavelengths  frees 
the  boundary  the  field  amplitude  associated  with  these  details  will 
always  be  less  than  the  ambient  noise  amplitude.  For  example,  conven¬ 
tional  inages  forced  with  light  say  contain  infbrsati®  densities  ex- 
ceediijg  10"  bits/cm  of  retrievable  information-,  informatics  densities 

S  ■?  * 

as  high  as  10  hits/ca"  have  been  reccroed  by  photographic  seas. 

*  This  figure  is  representative  of  the  lodak  series  649  ssulsiass 
which  are  known  to  resolve  in  excess  of  20-30  lines/ss.  The  aver¬ 
age  black  and  white  snapshot  fils,  such  as  Sodas  Fanatomic-X  can 
record  Isnge  detail  to  100  linss/ss  or  10^  bits/cn^.  (See 
Reference  2.)  For  cccrpari son,  high  density  integrated  circuit 
<cxpeixcccuc ting  cecputer  memories  are  presently  capable  of 
storing  approedsatel}  2300  bits/css?.  (See  Reference  3.) 
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Information  densities  of  this  ragnitude  era  typical  of 
scenes*  Resolution  of  solecular  or  stosic  detail  requires  illusissticsj 
at  wavelengths  such  shorter  than  that  of  visible  light,  say  at  x-rsy  cr 
electron  wavelengths. 

Information  transmitting  capacity  of  such  magnitude  is  or  is- 
estimable  value  in  scientific  isvestigatifsj  and  exploration.  However, 
the  in&nsation  contained  in  wave  fields  is  often  vastly  greater  than 

the  capacity  of  the  eye  cr  other  device  to  receive  directly.  For  ex?®- 

* 

pie,  nature1  virtual  i rages  of  the  cosmic,  rderosccpic  and  atomic 
regions  of  the  universe,  always  of  great  interest  to  the  scientist,  con¬ 
tain  information  densities  far  exceeding  that  which  the  eye  can  detect 
and  adequately  resolve  without  ami  11  ary  apparatus.  Sbch  regions  are 
usually  regarded  as  ’’invisible.” 

Quite  apart  fror  resolution  and  field  of  view  restricticss  im¬ 
posed  by  any  re?. li cable  receiving  cr  display  process,  the  infcnstica 
which  can  be  transmitted  fras  the  inject  into  sn  optics!  display  is 


>"  'll  iiiiimiW|[IIMI»»Kn»WII«,'i  . . . "... 


although  intrinsic  absorption  by  the  medium  may  be  small ,  randan  turbu¬ 
lence  induces  scattering  which  limits  the  detail  which  can  be  resolved 
through  such  media.  For  example,  observation  of  stellar  bodies  from 
the  Earth's  surface  is  severely  hampered  by  uncontrollable  fluctuations 
of  the  terrestrial  atmosphere^. 

Hence,  from  the  time  of  Galileo  (the  telescope)  and  Leeuwenhoek 

(the  microscope)  one  of  the  great  challenges  of  experimental  science 

has  been  to  render  visible  those  regions  beyond  the  ability  of  the  eye 

to  perceive;  more  technically,  to  extend  visual  perception  into  regions 

* 

outside  the  space  bandwidth  product  of  the  eye.  The  quest  for  greater 
visibility  has  led  to  imaging  with  fields  to  which  the  eye  does  not 
directly  respond;  to  wit,  x-rays  (crystallography),  electrons  (electron 
and  field  emission  microscopy) ,  acoustic  radiation  (acoustic  holo¬ 
graphy)  ,  etc.  The  emphasis  here  is  directed  toward  acoustic  imaging 
by  which  is  meant  the  visualization  through  acoustic  means  of  objects 
imbedded  in  solid  or  liquid  media  usually  opaque  to  light. 

Interest  in  acoustic  imaging  has  arisen  in  response  to  rapidly 
expanding  activity  in  two  relatively  new  but  quite  unrelated  fields, 
biomedical  engineering  and  oceanography.  The  use  of  x-rays  in  medical 
diagnosis  has  always  been  somewhat  problematical.  Not  only  does  x-ray 
imaging  subject  the  patient  to  extremely  hazardous  radiation,  but  the 
resulting  images  often  do  not  fully  satisfy  the  requirements  of  the 
examiner.  For  example,  the  technique  does  not  readily  distinguish 
different  types  of  tissue.  Further,  because  it  is  difficult  to  focus 


This  concept  is  discussed  in  Chapter  II. 
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x-rays,  the  display  is  actually  an  x-ray  shadow  rather  than  a  truly 
three  dimensional  focused  image  of  internal  body  structure.  Acoustic 
imaging  need  not  possess  either  of  these  drawbacks. 

Acoustic  imaging  may  also  play  a  role  in  an  entirely  different 
endeavor,  the  development  of  underwater  resources.  Lack  of  visibility 
is  becoming  an  impediment  to  more  rapid  progress  in  undersea  explora¬ 
tion,  mining,  and  transportation.  Even  in  absolutely  pure  water,  light 
will  not  propagate  much  further  than  perhaps  100  meters.  Worse,  one  is 
usually  most  interested  in  operating  along  the  ocean  floor  where  the 
water  has  likely  been  churned  into  murky  clouds  of  mud  stirred  up  from 
the  bottom.  The  range  of  x-rays  in  underwater  situations  is  limited  to 
the  order  of  a  few  centimeters.  On  the  other  hand,  acoustic  fields  can 
propagate  through  water  for  hundreds  to  thousands  of  meters,  depending 
upon  frequency,  without  excessive  attenuation.  Thus,  for  oceanographic 
applications  acoustic  imaging  assumes  importance  as  a  means  of  extending 
visibility  beyond  the  range  of  electromagnetic  radiation. 

Acoustic  imaging  is  an  infant  science,  at  a  stage  of  development 
which  may  be  likened  to  that  of  conventional  optics  when  Galileo  first 
observed  the  moons  of  Jupiter,  A  multitude  of  challenges  must  be  over¬ 
come  in  order  to  make  acoustic  imaging  viable  with  other  investigative 

* 

techniques.  The  goal  is  an  acoustic  system  which  provides  in  real  time 
visualization  with  approximately  the  quality  of  inexpensive  camera 


Standard  broadcast  television  frame  rates  (30/sec)  will  satisfy 
the  real  time  requirement.  However,  this  is  probably  restrictive 
and  8irm  movie  frame  rates  (16/sec)  or  lower  (8-10/sec)  will  be 
sufficient. 


snapshots  or  standard  television  displays.  The  sophisticated  refine¬ 
ments  of  light  optics  are  not  required. 

If  acoustic  imaging  is  to  proceed  along  conventional  lines  then 
an  acoustic  lens  or  other  focusing  device  will  be  needed.  However, 
such  devices  present  formidable  development  and  fabrication  problems . 
Consider  a  typical  camera  apperture;  a  comparable  apperture  for  1  MHz 
acoustic  radiation  v/ould  be  a  thousand  times  larger,  measuring  perhaps 
100  meters  or  more  in  diameter.  Such  a  device  could  hardly  be  con¬ 
sidered  portable.  Further,  the  scale  of  the  object  scene  will  be  the 
same  whether  the  illumination  is  acoustic  or  visible.  Hence  the  F 
number  (focal  length  to  diameter  ratio)  of  ideal  acoustic  lenses  must 
be  quite  small  and  the  images  will  then  be  subject  to  the  severe 
afcberations  and  distortions  inherent  in  nonparaxial  imaging.  There¬ 
fore,  for  practical  reasons  acoustic  lenses  will  very  likely  be  ap¬ 
preciably  smaller  than  ideal  and  the  resulting  displays  will  be 
comparable  to  those  obtained  by  pinhole  imaging  with  light. 

Recently,  an  entirely  new  and  different  method  of  imaging  which 
does  not  require  lenses,  the  two  step  process  invented  by  Dennis 
Gabor^,®,^,“®'1'* and  known  as  holography,  has  become  the  subject  of 
intense  research  activity  in  visible  optics.  Because  Gabor's  technique 
does  not  require  lenses,  crude  experiments  in  acoustic  imaging  can 
actually  be  performed  more  easily  by  adopting  the  methods  of  holography 
rather  than  conventional  methods,  employing  auxiliary  focusing  apparatus. 
Ironically,  this  is  just  the  reverse  of  the  situation  in  visible  optics 
where  good  lenses  are  easily  fabricated  but  holograms  require  carefully 
controlled  laboratory  setups.  Of  course,  the  reason  is  clear. 
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Holography  requires  coherent  radiation.  Coherent  acoustic  radiation  is 
easily  generated  by  crystal  controlled  oscillators.  Conversely, 
coherent  light  usually  requires  a  laser  source  and  special  precautions 
to  ensure  that  the  illumination  is  sufficiently  coherent  to  produce 
high  quality  hologi*ams.  We  shall  present  a  study  of  acoustic  imaging 
by  holography;  i.e.,  acoustic  holography. 

In  order  to  include  holography  in  optical  discussions,  the  term 
"image"  must  be  replaced  with  the  broader  concept  of  optical  representa¬ 
tion.  By  an  optical  representation  of  any  object  we  include  not  only 
those  referred  to  conventionally  as  "images"  which  bear  some  geo¬ 
metrical  similarity  to  the  object,  but,  more  generally,  ary  field  con¬ 
figuration  which  may  result  from  an  illuminated  or  l’tminous  boundary cr 
volume.  For  example,  a  snapshot  is  a  conventional  image  but  a  hologram 
is  an  optical  representation  in  the  more  general  sense;  holograms 
usually  bear  no  resemblance  to  the  objects  from  which  they  are  formed 

although  they  actually  contain  more  information  about  those  objects 

* 

than  do  conventional  photographs.  The  entire  process  of  forming 
optical  representations  of  any  type  shall  be  referred  to  loosely  as 
imaging. 

The  plan  for  the  discussion  of  acoustic  imaging  to  follow  is 


There  has  been  some  confusion  surrounding  this  point.  It  is  often 
concluded  that  the  reconstruction  from  a  hologram  contains  more 
information  than  the  image  from  a  lens.  This  is  not  correct;  the 
lens  image  is  also  three  dimensional.  However,  the  three 
dimensional  nature  of  the  lens  image  cannot  be  recorded  on  a two 
dimensional  medium  without  resorting  to  special  methods  such  as 
holography  or  stereoscopy.  (See  Reference  13.) 
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simple.  In  Chapter  II  the  elements  of  scalar  imaging  are  presented  as 
they  apply  to  simple  acoustic  imaging  systems;  the  basic  imaging 
problem  is  to  determine  the  causal  relation  between  contrast  of  object 
boundary  field  spatial  structure  and  contrast  of  coi responding  image 
field  structure.  In  general,  image  contrast  is  a  decreasing  function 
of  spatial  frequency;  object  structure  for  which  the  corresponding 
image  contrast  vanishes  is  not  resolved.  In  acoustic  imaging  it  is 
especially  important  to  maximize  the  resolution  obtainable  from  a 
given  aperture  because  of  the  relatively  long  wavelength  of  acoustic 
radiation.  No  doubt  image  resolution  and  contrast  is  one  of  the  major 
criteria  by  which  acoustic  imaging  systems  must  be  evaluated. 

The  properties  of  the  detector  and  the  statistical  behavior  of 
the  scalar  image  field  significantly  constrain  the  contrast  transfer 
characteristics,  hence  resolution,  of  imaging  systems .  In  this  sense 
acoustic  imaging  systems  are  more  interesting  than  optical  imaging 
systems .  To  wit,  imaging  with  light  is  restricted  to  detectors  whose 
output  is  seme  long  time  average  of  the  electromagnetic  field  energy. 
Conversely,  the  piezoelectric  and  ferroelectric  compounds  form  an 
important  class  of  acoustic  radiation  detectors  for  which  the  output  is 
proportional  to  instantaneous  acoustic  pressure  ampltidue^.  This 
allows  latitude  in  the  design  of  acoustic  imaging  systems  which  is  not 
possible  for  optical  imaging  systems. 

The  significance  of  the  flexibility  inherent  in  acoustic  imaging 
with  arrays  of  piezoelectric  elements  becomes  apparent  when  the 
statistical  time  evolution  of  the  field  is  taken  into  account.  For 
detector  outputs  which  are  long  time  averages  of  the  field  energy, 


8 


spatially  incoherent  imaging  tends  to  provide  greater  ultimate 
resolution  but  less  contrast  than  coherent  imaging.  In  addition, 
incoherent  images  do  not  suffer  from  edge  ringing  and  other  inter¬ 
ference  effects.  By  electronically  synthesizing  higher  order  moments , 
say  the  fourth  order,  it  should  be  possible  to  further  improve  the 
resolution  of  incoherent  acoustic  imaging  systems.  Moreover,  by 
employing  scanning  systems  such  as  the  Covington-Drane  antenna  rather 
than  lens  systems  the  contrast  of  the  incoherent  acoustic  system  can  be 
made  to  approach  that  of  the  coherent  system  but  with  twice  the 
resolution. 

Heretofore,  acoustic  imaging  has  been  confined  to  spatially 
coherent  point  source  acoustic  ilhmination  and  using  holographic 
techniques.  It  seems  clear  that  a  new  class  of  experiments  in  acoustic 
imaging  is  called  for,  namely  using  spatially  incoherent,  extended 
acoustic  sources,  lenses  or  other  focusing  devices,  and  a  variety  of 
electronically  synthesized  detection  characteristics.  Such  experiments 
have  scientific  merit  without  regard  to  practical  applications  since 
similar  experiments  cannot  be  performed  optically.  For  example,  by 
varying  the  averaging  time  of  the  detector  output,  easily  accomplished 
electronically,  the  effective  degree  of  coherence  of  an  acoustic 
imaging  system  can  be  varied  continuously  from  completely  incoherent  to 
completely  coherent. 

The  other  major  factor  to  be  considered  in  judging  the  viability 
of  acoustic  imaging  systems  is  the  adaptability  of  the  imaging  process 
to  real  time  imaging.  Holography  is  basically  a  two-step  process 
therefore  is  not  readily  adaptable  to  real  time  imaging  with  piezo- 
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electric  detectors.  In  this  regard  it  should  also  be  noted  that  since 
incoherent  imaging  requires  finite  detector  averaging  times  the  frame 
rate  at  which  acoustic  images  are  formed  may  be  less  than  one  might 
desire . 

The  variety  of  acoustic  imaging  systems  which  can  be  devised  is 
certainly  more  extensive  than  the  list  presented  in  Chapter  II.  In 
this  context  acoustic  holography  must  then  be  regarded  simply  as  one  of 
many  possible  acoustic  imaging  schemes  and  it  seems  rather  premature 
to  assess  the  viability  of  acoustic  holography  until  ether  imaging 
methods  have  been  investigated  experimentally. 

Difficulties  of  unique  importance  in  acoustic  imaging  are  posed 
in  Chapter  III.  Specular  reflection  and  attenuation  will  significantly 
degrade  the  quality  of  images  obtained  by  long  range  acoustic  imaging 
of  objects  in  the  sea.  Specular  reflection  is  the  most  severe 
problem;  it  renders  most  objects  mere  highlights  of  extended  objects. 
The  problem  of  specular  reflection  can  be  eliminated  by  using  diffuse, 
extended  sources  of  acoustic  illumination.  Here  again  incoherent 
acoustic  sources,  either  natural  or  man  made,  may  prove  superior  to 
coherent  sources. 

In  Chapter  IV  a  theory  of  holography  is  formulated  in  terms  of  an 

eigenvalue  equation.  This  type  of  equation  represents  mathematically 

the  phenomena  of  wave  front  reconstruction.  It  is  interesting  to  note 

*  In  order  to  provide  real  time  acoustic  imaging  by  holographic 

means  and  employing  piezoelectric  detectors  it  will  be  necessary 
to  form  the  hologram  on  some  photosensitive  material  such  as  the 
heat  developing  film  manufactured  by  The  Minnesota  Mining  and 
Manufacturing  Company  which  can  be  dry-developed  essentially 
instantaneous ly . 
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that  equations  of  this  type  play  a  central  role  in  many  scientific 
disiplines.  Following  an  outline  of  the  eigenvalue  formalism,  the 
holographic  process  is  describes  in  terms  of  a  diffraction  theory 
based  upon  expanding  the  fields  in  plane  wave  spectra  rather  than  the 
more  usual  treatment  employing  a  double  application  of  the  Rayleigh- 
Sonmerf eld- Green  or  Fresnel-Kirchoff  diffraction  integrals.  Plane 
wave  expansions  are  used  because  the  situations  encountered  in  acoustic 
imaging  are  often  non-paraxial .  Special  properties  of ,  the  detector 
which  do  not  arise  in  the  theory  of  optical  holography  but  which  imost 
be  considered  in  acoustic  holography  are  included  in  the  theory.  If 
certain  simplifying  assumptions  are  made  then  the  eigenvalue  equation 
possesses  a  continuum  of  solutions.  In  the  more  genral  non-paraxial 
case  it  is  not  known  if  there  are  any  solutions:  i.e.,  if  perfect 
reconstruction  is  possible. 

Finally,  in  Chapter  V  the  results  of  experiments  are  presented 
which  demonstrate  that  images  of  underwater  objects  may  indeed  be 
obtained  by  acoustic  means  employing  the  techniques  of  holography. 

These  experiments  involve  a  scanning  procedure  which  proved  expedient 
for  the  purpose  of  demonstrating  the  feasibility  of  acoustic  holography 
but  which  does  not  bear  practical  merit.  Special  effort  was  expended 
to  make  the  experiments  representative  of  the  type  of  situations  that 
would  be  encountered  in  practical  undersea  acoustic  imaging  problems. 
For  example,  the  frequency  was  low  enough  (1  >Kz,  X  =  1.5  inn)  to 
provide  a  range  of  perhaps  50  -  100  m  in  dear  sea  water,  the  aperture 
was  a  realistic  size  (approximately  1.5  m  square),  objects  were  imaged 
in  reflection  rather  than  backlighting;  and  the  ranges  were  1.5  -  5  m. 
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CHAPTER  II 


ELBESTS  OF  SCALAR  IMAGING 

1.  IMAGES,  A  GENERAL  PROBLBi  IN  SCALAR  KAVE  FfOPAGATICK 

Cities,  in  the  broad  sense  to  be  taken  here,  refers  to  the  study 
of  any  propagating  disturbance  (not  necessarily  electrasagnetic)  de- 
scribable  by  a  real  scalar  fraction  ?(r,t)  of  position  r  and  tiae  t, 
hereafter  called  the  optical  or  field,  which  satisfies  the  hcrao- 
genecus  scalar  wave  equation 

^(r,  t) - i— 7  «  0  fll.l) 

[c(r)]z  3t* 

subject  to  Cauchy  boradary  conditions  esi  open  space- tioe  boundaries*5. 
The  disturbance  propagates  thresh  a  medium  M  with  velocity  c(r) ,  the 
position  dependence  of  c  denoting  an  infcooogeneous  saterial.  More 
generally,  the  material  Kill  also  be  anisotropic  and  dispersive  and 
wave  propagation  therein  aore  cooplex  than  that  isplied  in  Equation 
(II.  1).  Nevertheless,  the  extensive  and  isportant  class  of  optical 
phenoDena  known  as  scalar  imaging  is  contained  in  Equation  (II. 1). 

Many  quantities  of  diverse  physical  origin  obey  the  scalar  wave 
equation,  at  least  approximately,  hence  qualify  as  optical  fields. 

Aacng  others  we  nay  list  snail  amplitude  acoustic  pressure,  velocity, 
and  displacement  fields  in  non-viscous  fluids,  and  the  components  of 
any  electromagnetic  field.  Indeed,  isider  certain  conditions,  electro¬ 
magnetic  radiation  in  the  visible  spectrum  (light)  is  fully  represented 
by  a  single  scalar  field  V(r,t)  referred  to  as  the  light  solitude16. 
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Imaging,  the  process  of  forming  optical  representations ,  is  an 
important  and  venerable  branch  of  optics  which,  historically.  Las  been 
associated  almost  entirely  with  light.  In  this  chapter  the  sore  impor¬ 
tant  aspects  of  scalar  imaging  theory  will  be  extended  to  enctspass 
related  acoustic  phenomena  as  well. 

A  general  imaging  situation  is  depicted  schematically  in  Figure 

XI.  1  .  A  sisply  or  multiply  connected  region  0,  the  object  space,  is 

’'illuminated”  by  a  wave  field  ?^(r,t)  radlatisig  free  luminous  source 

region  L.  The  object  scatters  a  field  fg(r,t)  which,  in  general ,  will 

propagate  in  all  directions.  After  passing  through  the  intervening 

* 

medium  M(r,t)  a  portion  of  ray  happen  to  fall  upon  sense  arrangesaeat 
of  material  F(rj ,  the  imaging  or  focusing  device.  The  inaging  device 
transforms  the  field  ^  in  the  space  0  into  a  new  field  in  the 
space  I,  the  image  space,  immediately  adjacent  to  F.  In  the  image 
space  Yj  forms  various  spatial  patterns  H(r,t)  which  are  optical  repre¬ 
sentations  for  the  object  0.  These  patterns  say  be  recorded  for  use  at 
a  later  tire. 

Fcr  any  given  object,  the  exact  configuration  of  H  has  consider¬ 
able  latitude.  It  may  be  geometrically  similar  to  the  object  in  which 
case  it  is  known  as  an  image  and  is  recorded  as  a  photograph.  However, 
in  one  important  situation,  holography,  H  is  an  interference  pattern 
bearing  little  or  no  resemblance  to  the  object;  a  record  of  this  inter¬ 
ference  pattern  is  called  a  hologram.  The  emphasis  in  the  remaining 
chapters  is  upon  holography.  However,  before  taking  up  this  study, 
certain  elementary  properties  characteristic  of  any  imaging  system 
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must  be  provided;  in  particular,  basic  relations  governing  resolution 
will  be  formulated. 

2.  THE  OBSERVABLES  OF  ACOUSTIC  IMAGING 

We  begin  by  defining  the  observables  of  any  imaging  system.  The 
general  propagation  problem  is  to  determine  ¥(r,t)  satisfying  Equation 
(II. 1)  and  specified  boundary  conditions  in  space- time.  The  imaging 
problem  requires  two  additional  constraints:  (1)  only  certain  func¬ 
tionals  D(r,t)  of  optical  fields  are  observable;  i.e.,  can  be  measured, 
and  (2)  ¥(r,t)  is  a  stochastic  function;  i.e.,  a  function  which  repre¬ 
sents  a  physical  quantity  exhibiting  random  fluctuations  in  space- time. 

Condition  (1)  merely  reflects  the  properties  of  the  apparatus 
used  to  detect  optical  fields,  namely  that  the  output  of  the  detection 
elements  may  not  respond  directly  to  the  field  amplitude,  rather  to 
various  physical  quantities  D  called  image  observables  which  are  func¬ 
tionals  of  the  field  amplitude.  Condition  (2)  is  a  statement  that  any 
physical  process  giving  rise  to  optical  fields  is  inherently  statis¬ 
tical  in  nature.  The  manner  in  which  the  field  statistics  will  influ¬ 
ence  imaging  depends  upon  the  characteristics  of  the  detection  process, 
the  stuffy  of  which  we  now  take  up. 

A  general  class  of  image  observables  may  be  written  as 

D(r,t)  =  j  f[T(r,t')]w(t-t*)dt'  (II. 2) 

-00 

where  f  represents  some  function  of  the  field  and  w(t-t')  is  a  temporal 
weighting  factor.  For  example,  the  only  detectors  known  to  be 
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sensitive  to  light  respond  to  incident  field  energy  flux  rather  than  to 
field  amplitude.  Moreover,  the  response  time  of  any  known  light  sensi¬ 
tive  detector  (photographic  emulsion,  bolometer,  photocell,  retina, 
etc.)  is  much  longer  than  the  period  of  oscillation  of  light,  hence  the 
detector  output  D  cannot  follow  the  instantaneous  oscillations  of  the 
light.  Accordingly,  Skinner  has  shown  that  only  the  time  averaged 
intensity  of  light  can  be  measured  and  the  appropriate  image  observable 
is 


D(r)  -«  <^['?(r,t]2^>  =  <^(r,t)/(r,t)^> 


(H.3) 


where  angular  brackets 


indicate  an  infinite  time  average 


*  lim 
T-h» 


(  )  dt 


(II. 4) 


and  superscript  denotes  complex  conjugation. 

Here  we  have  adopted  the  convention  that  the  complex  valued 

mathematical  fields  are  to  be  taken  as  the  analytic  signal  representa- 
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tions  of  the  real  valued  physical  fields  T  .  It  is  understood  that  at 
the  end  of  a  calculation  the  physical  quantity  is  regained  as  the  real 
part  of  the  analytic  signal,  thus  T  =  Reij>.  The  infinite  time  average 
is  a  mathematical  convenience.  It  should  be  interpreted  as  meaning  an 
average  over  times  long  compared  to  the  period  of  oscillation  of  the 
radiation  and  the  mean  time  between  statistical  fluctuations  cf  field 
amplitude  and/or  phase.  Devices  obeying  response  (II. 3)  are  known  as 
energy  or  square  ’aw  detectors.  Because  all  phase  information  is  lost 
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in.  the  time  averaging  process,  they  belong  to  the  general  class  of 
incoherent  detectors. 

Since  the  observable  is  intensity  I  not  field  it  is  desirable 
to  formulate  imaging  directly  in  terms  of  the  former  rather  than  the 
latter.  However,  this  is  not  possible  because  the  intensity  in  any 
region  is  not  uniquely  related  to  the  intensity  on  the  boundaries 
surrounding  that  region.  It  is  possible  to  introduce  a  new  function 
called  the  mutual  coherence  function^ ’ ^  r(rj ,  r^,  t^,  t^), 

rCl*.  r2,  tr  t2)  =  , t^t) / (r2 , t2+t)^>  ,  (11.5) 

which  contracts  to  the  intensity 

rCli»  1Li>  tx»  t2^  *  ICl»t)  (II. 6a) 


as 


II 


-2 


t,  t. 


(II. 6b) 


and.  for  which  the  propagation  problem  is  unique.  Moreover ,  because  r 
is  in  the  form  of  a  correlation,  the  statistics  of  the  field  are  con¬ 
veniently  included  in  this  function.  If  the  temporal  dependence  of  r 
is  only  through  the  time  difference  t  =  t2  -  tp  V  is  stationary  in 
time  and  we  may  then  write 


r(rr  x_2 *  *2^  1’  — 2*^  ~  ^12^ 


(II. 7) 
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The  equations  of  motion  for  r^CO  ate  well  known  ,  we  merely  repeat 
than  here 


2  i  **  ^12^ 

7iri2(T)  '  ~7  - 2 - =  0 

C  3t 


(II. 8) 


72F12^ 


1  9rl2^ 

7 


=  0 


Solutions  for  may  be  found  by  the  standard  techniques. 

Because  acoustic  fields  are  a  mechanical  phenomena,  electro¬ 
mechanical  processes  are  required  for  their  detection,  and  acoustic 
observables  may  have  properties  quite  different  from  the  intensity 
observable  of  light.  The  most  consnon  acoustic  detectors  operate 
through  the  phenomena  of  piezoelectricity  (quartz  or  any  other  crystal 
lacking  a  center  of  syrrmetry)  or  ferroelectricity  (barium  titanate, 
etc.).  We  shall  not  enter  into  a  detailed  discussion  of  these  mech¬ 
anisms  here  but  merely  point  out  that  with  either  phenomena  is  as¬ 
sociated  the  generation  of  an  electric  potential  linearly  related  to 

21 

the  incident  acoustic  stress  or  strain  .  Provided  the  detector  is 
suitably  terminated,  the  frequency  of  acoustic  oscillations  is  so  low 
that  the  potential  across  the  crystal  will  follow  the  instantaneous 
acoustic  field.  Accordingly,  for  acoustic  fields  the  fundamental 
observable  is 


D(r,t)  «  ¥ (r , t)  .  (II. 9) 

Response  law  (II. 9)  belongs  to  the  class  of  coherent  detectors 
and  is  especially  important  because  many  other  response  laws  can  be 
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electronically  synthesized  from  it.  Among  others  we  distinguish  the 
square  law  detector,  Equation  (11,3),  and  the  correlation  detector 

D(r,t)  «  ¥(r,t)^tS(r,t)  (11.10) 

where  denotes  time  correlation  according  to  the  prescription  given 
in  Appendix  B,  S(r,t)  is  a  "reference"  against  which  the  signal  is  ccxn- 
pared,  and  '*'(r,t)  is  assumed  to  be  of  finite  duration.  A  special  case 
of  (II. 10)  is  the  time-matched,  filter  for  which 

D(r,t)  cc  ¥(r,t)*t¥(rft)  .  (II. 11) 

Finally  higher  order  detection  processes,  for  example 

D(r)  «  <jnr,y)]4>  ,  (11.12) 

are  possible. 

3.  TEMPORAL  STRUCTURE  OF  WAVE  FIELDS 

Mathematically,  imaging  is  simply  finding  solutions  of  the 
propagation  problem,  Section  (1),  for  the  image  observables,  Section 
(2).  Since  instantanecus  amplitude  is  observable  in  acoustic  imaging, 
a  formal  analysis  directly  in  terms  of  field  rather  than  mutual  coher¬ 
ence  is  valid.  In  view  of  the  linearity  of  Equation  (II. 1),  the 
imaging  process  can  be  formulated  as  a  superposition  using  appropriate 
Green's  functions.  To  see  how  the  superposition  should  be  accomplished 
we  must  consider  the  temporal  structure  of  the  fields. 
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Wave  fields  are  either  polychromatic  or  quasi -monochromatic; 
purely  monochromatic  fields  are  inpossible.  For  quasi-mono  chroma  tic 
radiation  the  temporal  spectrum  5»(r,v)  is  confined  to  a  narrow  fre¬ 
quency  range  Av  about  a  central  line  v  such  that  Av«v,  and  may  be 
written 


$(r,v)  =  i(r,v)*t6(v-«) 


(11.13} 


where 


£(r,v)  =  0  for  jv-v]  ^  Av/2  , 


(11.14) 
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6  is  the  Dirac  delta  function  ,  and  *t  denotes  time  correlation  ac¬ 
cording  to  the  prescription  given  in  Appendix  B.  Acoustic  fields  with 
sufficient  spectral  purity  to  satisfy  the  quasi-monochromatic  approxi¬ 
mation  are  easily  generated;  indeed,  spectral  purities  exceeding  1  part 
in  106  are  connton.  Unless  otherwise  noted,  the  quasi-monochromatic 
approximation  is  hereafter  assuned  valid. 

Applying  relation  (B.7)  for  the  inverse  Fourier  Transform  of  a 
convolution  yields 


❖(r,t)  =  5(r,t)e 


-27rivt 


(11.15) 


the  analytic  signal  representation  for  quasi-monochrcmatic  fields. 

It  is  a  simple  matter  to  show  that  £(r,t),  a  function  varying  slowly 
with  time,  obeys  the  relation 


v2c(r,t)  +  k2C(r,t)  =  0 


(11.16) 
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where  k2  =  -  .  It  is  then  clear  that  quasi-monochrmatic  fields 

c 

can  be  treated  as  purely  monochromatic  if  the  temporal  variation  of  £ 
is  small  during  some  time  interval  x  characteristic  of  the  specific 
problem  to  which  (11.16)  is  applied.  The  time  depend !.-\ce  of  £  can 
then  be  omitted  in  Equation  (11.16)  and  the  result  recognized  Sis  the 
Helmholtz  equation 


A(r)  +  k2£(r)  =  0 


(11.17) 


long  familiar  as  the  equation  of  propagation  for  purely  monochromatic 
radiation. 

Now  the  only  time  spans  of  significance  in  imaging  are  those 
required  for  an  optical  disturbance  to  propagate  over  the  possible 
optical  paths  from  object  to  image.  Calling  the  maximum  optical  path 
length  L,  the  characteristic  time  interval  for  imaging  problems  is  no 
greater  than 


_  L 
T  c 


(11.18) 


We  then  require  that  the  time  x’  for  £  to  undergo  significant  change 
be  restricted  so  that 


x*  »  4 


(11.19) 


An  approximate  value  for  t  '  can  be  obtained  by  considering  T  as 
isisting  of  wave  trains  of  duration  x  * .  Then 
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£  =  Const 


|tj  <  x/2 


(11.20) 


=  0  i  tj  >  t/2 


(11.21) 


(11.22) 


Therefore,  in  order  to  treat  the  fields  as  purely  monochromatic  we 
demand 


Av  «  ~  or  L  «  .  (11.23) 

L  Av 

For  example,  if  the  radiation  is  a  1  Miz  acoustic  signal  with  spectral 
purity  of  1  part  in  106  then  L  is  approximately  1500  meters  in  water. 

The  foregoing  has  been  concerned  only  with  determining  the  proper 
wave  equation.  To  complete  the  discussion  requires  additional  comment 
on  the  statistical  behavior  of  the  field,  in  particular  on  the  compara¬ 
tive  time  evolution  of  the  field  at  two  arbitrary  points  r. ,  and  r^  in 
space.  The  parameter  which  measures  this  is  the  two  point  temporal 
cross  correlation  or  mutual  coherence  function  r^CO  previously 
defined,  or  a  related  quantity,  the  complex  degree  of  coherence  y-^W 
which  is  normalized  so  that 
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(11.24) 


9 


0  ~  -  1 


By  definition 


Yi2Ct)  - 


ri2Cx) 


s  1/2 


(11.25) 


[riiCoDr22ccor 

If  thu  two  points  r,  and  are  allowed  to  contract  to  a  single 
point  in  space  the  normalized  autocorrelation  y-q(~)  of 

the  field  at  that  ooint;  autocorrelation  is  a  measure  of  the  temporal 
coherence  properties  of  the  field.  If  the  time  evolution  of  the  phase 
of  the  field  is  perfectly  linear  and  deterministic  then  the  time 
structure  of  the  field  is  correlated  for  all  time  differences  t  and 


|y1x(t) j  =  1 


(11.26) 


We  then  speak  of  a  temporally  coherent  field.  This  condition  can 

exist  only  if  the  field  is  purely  monochromatic.  For  quasi- 

monochromatic  fields  the  autocorrelation  vanishes  beyond  sane  interval 

x_  called  the  coherence  time;  for  auasi-monochromatic  fields 
c - 


iYn(t)i  »  i 


=  0 

we  identify  x'  with  xc. 


x  «  x 


X  >x_ 


(11.27) 


Temporal  coherence  plays  an  important  role  in  interference 
phenomena,  particularly  those  involving  the  self- interference  of 
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optical  beams.  For  example,  the  speckled  images  characteristic  of 
coherently  illuminated  objects  occur  when  the  scattered  illumination 
interfers  with  the  incident  illumination.  This  can  happen  as  long  as 
the  optical  path  length  difference  AL  satisfies 


Al  «  — 
Av 


(11.28) 


Since  path  length  differences  rather  than  absolute  path  lengths  are 
involved  in  (11.28)  the  latter  is  less  restrictive  than  the  condition 
for  validity  of  the  Helmholtz  equation.  Thus,  we  usually  assume  that 
quasi -monochromatic  radiation  is  temporally  coherent. 

The  other  limiting  form  for  the  complex  degree  of  coherence  is 
of  more  importance  in  imaging.  If  the  time  difference  t  is  allowed  to 
vanish,  Y^2(t)  becomes  the  normalized  mutual  intensity  Yj,2(Q)  ;  mutual 
intensity  is  a  measure  of  the  spatial  coherence  properties  of  the  field. 
If  the  time  evolutions  of  the  phase  of  the  field  at  two  points  in  space 
are  in  unison  then  the  time  structure  of  the  field  at  the  two  points  is 
correlated  and 


|V12(0)!  =  1,  all  rx  and  r2  .  (11.29) 

We  then  speak  of  spatially  coherent  fields.  If  the  two  time  evolutions 
are  independent  then  the  time  structure  of  the  field  at  the  two  points 
is  uncorrelated  and 


!y12(o) !  "  0,  *  r,  . 

We  then  speak  of  spatially  incoherent  fields. 


(11.30) 


25 


*>)"  •,W>fihWWlTl ^}i'<¥gjfl*nyi»V»flrtint  ;»  »^n;  iV'-'iw'-i*  '  <i'  -n« '  1  '*<:< 'w\ 1 »?  ■>  >m'>  *, 1  ii'  t|  >i  n  t."HI  * .  -mi-  IiY" ‘jT'Mit'Hii't  ll1"i»1’  ^i1'0  v*'  <)l 


It  is  physically  impossible  to  obtain  a  field  which  is  either 

entirely  coherent  or  entirely  incoherent.  The  former  obtains  only  for 

purely  monochromatic  fields.  For  quasi-monochrcmatic  radiation 
17 

Skinner  has  denonstrated  that  an  arbitrarily  narrow  power  spectnsn 

does  not  guarantee  spatial  coherence  anywhere,  a  somewhat  surprising 

24 

result  in  view  of  the  van  Cittert-Zemike  Theorem  .  Conversely, 

25 

Beran  and  Parrent  have  shown  that  any  field  for  which  Y^CO)  =  0 
everywhere  does  not  propagate,  hence  cannot  exist  except  possibly  as 
rapidly  decaying  evanescent  waves.  Hence,  for  all  quasi-mcncchrcmatic 
fields  there  will  be  some  region  over  which  spatial  coherence  will 
obtain.  Thus,  in  general 


where  the  condition  a  >  1  is  necessary  in  order  that  the  field  radiate. 
If  the  field  is  spatially  coherent  o  »  1;  if  spatially  incoherent, 
aSl. 

4.  WAVE  THEOPY  OF  IMAGING 

Although  neither  completely  coherent  nor  completely  incoherent 
fields  are  physically  realizable,  there  are  many  situations  of  practical 
importance  in  which  the  fields  are  approximately  coherent  or  approxi¬ 
mately  incoherent  over  the  spatial  region  of  interest.  We  now  study 
the  influence  of  the  field  statistics  upon  imaging  processes  for  these 
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limiting  cases.  Hie  aotr-'atien  for  this  investigation  arises  £rtss  a 
need  to  maximize  acoustic  image  quality. 

Acoustic  imaging  is  beset  with  a  ^altitude  or  lisdtaticss  of 
practical  origin,  most  being  attributed  to  tbs  relatively  locg  kss- 
length  of  acoustic  radiation.  Tie  result  is  acoustic  images  of  infe¬ 
rior  quality,  often  barely  discernible.  However ,  sizes  the  frssdagssBtai 
acoustic  abservablc  is  instantaneous  field  ampii-uide  ratter  than  field 
intensity,  acoustic  imaging  is  not  limited  to  energy  detectics  pro¬ 
cesses  as  is  light  inaging.  Otter  processes  say  be  synthesized 
electronically,  he  wish  to  determine  if  image  quality  ray  be  Isprevsd 
over  that  which  can  be  obtained  with  energy  detection. 

to  the  basis  of  Equation  (11.17)  the  fcrmatics  of  an  image  wits 
quasi-aonochrcfsatic  radiation  can  be  gives  as  a  spatial  superposition 
of  fields.  To  wit.  see  Figure  II- 2 


where  Cj(rj,t)  is  the  image  field  asplittee,  ^fr^.t)  is  the  object 
field  amplitude,  and  G(r.  ,r^)  is  the  Green*  s  risseties  for  the- 
process;  i.e. ,  the  field  solitude  at  image  coordinate  r.  is  image 
plane  I  resulting  frat  a  point  source  at  object  coordinate  r^  is  object 
plane  0.  For  simplicity  assume  that  the  image  formation  process  is 
stationary  in  spatial  coordinates;  i.e. . 


ctej.io)  *  ccrficP  • 


(11.33} 


Then 


Clfept)  *  f  ^(rg.tJGtrj-^d2^  (11.34) 

0 

which  is  valid  for  isoplanatic  image  formation. 

To  proceed  further  requires  -hat  the  process  be  stated  in  terms 
of  image  observables.  For  light  the  observable  is 

D(r)  =  4(r,tK  Cr,t)>  =  I(r)  .  (11.35) 

Substituting  from  (11.34)  and  canmuting  the  spatial  integration  with 
the  time  average  yields 


■Ml 5-/ /  ®(£.I-SD)G*(rI-iJ)  <«0(r6,t){o(^,tJ>  d^/rj  . 


(11.36) 


At  this  point  the  statistical  nature  of  the  optical  field  must  be 
introduced.  We  distinguish  two  cases.  If  the  light  is  spatially  in¬ 
coherent  the  time  evolutions  of  the  fields  ^(r^t)  and  CQ(r^#t)  are 
statistically  independent  and 


~  6vJo“^  <\^o^Io’t)^0^^0,ti^>  *  cn*37) 


Then 
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^(Ij)  =  = 

//w  drio^V^  <\^o^^o,t-*  ^o^,t^/>  d  rod  ro 

0 

/l^-Xo)!  <\^(lo,t^^,t^/>d  rO  *  (II. 38) 

0 

Equation  (11.38)  may  be  written  in  the  compact  form 

Ij(r)  =  I0(r)  *  |G(r)|2  ,  (11.39) 

the  isoplanatic  imaging  law  for  incoherent  optical  fields  when  the 
image  observable  in  the  output  of  a  square  law  detector.  Clearly,  in 
this  case  the  image  formation  process  is  linear  in  intensity.  By 
employing  suitable  electronics,  acoustic  imaging  systems  may  be 
devised  which  also  obey  Equation  (11.39). 

If  the  field  is  spatially  coherent,  the  time  evolutions  of 
VSb  ,t)  and  ^(r^.t)  are  identical  except  possibly  for  a  uniform  trans¬ 
lation  r(r)  in  time.  To  describe  this  synchronization  it  is  convenient 
to  reference  all  fields  to  sane  standard  Ss(t)  which  represents  the 
characteristic  structure  of  the  time  evolution  and  which  is  normalized 
so  that 

<?sCt)Cg(t)>  =  1  .  (11.40) 
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Then 


Soteo,t:)  =  ?0{^5S^t+T^^  (11.41) 

£0(r£,t)  =  ^(r^ytn^)]  ,  (11.42) 

where  Cq^q)  and  ^(Tq)  are  real  constants  representing  the 
"magnitudes"  of  the  fields.  To  find  the  £q(t)  multiply  the  appropriate 
defining  equation  by  its  complex  conjugate,  time  average  the  resulting 
product,  and  recognize  that  infinite  time  averages  are  invariant  under 
time  translations,  thereby  obtaining 

Cq(I<))  =  %/  (11.43) 

Cq(£q)  =  J  '(^(lo » £q(£q  » •  (11.44) 

The  correlation  between  fields  then  becomes 

/  *  v  (H.45) 

*  \%tt+TCzb)35s[t+T(^6)]/  * 

By  suitable  change  of  variables  we  obtain 

<^s  [ t+T  (Iq)  3  4  [  t+t  (Iq)  f>  =  <^s  Ol)  4 1 n+T  (hD3  * 1  do3 C11*463 
26 

Edwards  and  Parrent  have  shown  that  the  normalized  autocorrelation  of 
a  coherent  quasi-monochromatic  function  may  be  written 
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<5{r,tK*(r,t+Tj>  .  ei2™ 


(11.47) 


Thus 


iZirlrCr^D-TCiQ)] 


(11.48) 


Substituting  (11.48)  into  (11.35)  we  obtain  for  coherent  imaging 


W ' 


/gCIj-Io)/  <5oCib- t)5o(lo.t)>  a'Z,,1T(lo)d2r0 


to 


/<? <!*■$> J< «0(iJ.t)«o(TS,t)>  e2,'1Tti:0)d2r6 


to 


(11.49) 


Note  that  Ij(r.j)  has  been  separated  into  independent  and  r^  integra¬ 
tions.  Thus 


W  ■ 


/GHrio  )  /<^T> 


2\  '2irixCio)  ? 


dr. 


(11.50) 


Defining  an  equivalent  time  independent  complex  field  amplitude  as 


U(r)  =  /<j?(r.t)|2>  <f2™©  . 


(11.51) 


Equation  (11.50)  may  be  written  in  the  compact  form 


life)  =  i«bfe)*cfe)r 


(11.52) 
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the  isoplanatic  imaging  law  for  coherent  optical  fields  when  the  image 
observable  is  the  output  of  a  square  l?w  detector.  Clearly,  in  this 
case  the  image  formation  process  is  nonlinear  in  intensity;  in  fact, 
image  intensity  is  not  even  directly  related  to  object  intensity, 
rather  to  a  complicated  functional  of  the  complex  valued  amplitude. 

By  employing  suitable  electronics,  acoustic  imaging  systems  may  be 
devised  which  also  obey  Equation  (11.52). 

5.  RESOLUTION  AND  FIELD  STATISTICS 

We  now  compare  the  characteristics  of  coherent  (11.52)  and  in¬ 
coherent  (11.39)  imaging.  When  the  radiation  is  light  there  are  many 
ways  to  make  the  comparison,  depending  upon  the  application,  and  it 
is  not  possible  to  formulate  a  general  statement  classifying  one  type 
of  illumination  as  always  "better"  than  the  other.  However,  when  the 
radiation  is  acoustic  there  is  no  doubt  concerning  tl«.  choice  of 
criteria  by  which  image  quality  is  to  be  judged.  Restricted  apperture 
size,  long  wavelengths,  and  coarse  detector  resolution  severely  limit 
the  resolution  of  acoustic  imaging  systems.  Clearly  then,  the  problem 
of  maximizing  the  degree  to  which  details  of  the  object  structure  are 
resolved  as  distinct  and  similar  details  in  the  image  structure  must  be 
considered  of  prime  importance  in  acoustic  imaging.  Thus,  it  is  ap¬ 
propriate  to  examine  here  the  relation  between  resolution  and  field 
statistics. 

To  compare  the  resolution  of  coherent  and  incoherent  imaging 
systems  begin  by  taking  the  Spatial  Fourier  transforms  of  Equations 
(11.39)  and  (11.52)  (see  Appendix  A).  From  the  former 
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(11.53) 


life)  *  I0(s)(G(s)*G*(s)] 
and  from  the  latter 

IxCs)  =  [U0(s)G(s)3  *  [U*(s)G*(s)3  .  (11.54) 


The  physical  interpretation  of  Equation  (11.53)  is  clear.  If  the 
radiation  is  incoherent,  each  spectral  component  Iq(s)  of  the  object 
boundary  intensity  has  a  corresponding  spectral  counterpart  Ij(s)  in 
the  image  intensity  except  for  multiplication  by  the  complex  factor 
G(s)*G  (s)  =  H(s);  H(s)  is  called  the  spatial  frequency  transfer 
function. 

The  zero- frequency  spectral  component  1^(0)  represents  an 
intensity  uniformly  distributed  in  space,  and  requires  special  con¬ 
sideration.  Since  Ij(r)  is  non-negative,  lj(0)  >  0  unless  Ij(r)  =  0 
everywhere,  a  trivial  case.  Thus,  there  will  always  be  a  background 
intensity  Ij(0)  distributed  uniformly  over  any  image  formed  with 
incoherent  radiation.  It  is  convenient  to  normalize  each  spectral 
component  of  the  image  intensity  to  this  background.  Defining 


jig) 

lj(0) 


i0(°) 


^  G(0) 


m.  =  nr 


|^=H(s)  (11.55) 
H(0) 


we  obtain  the  normalized  form  of  (11.53) 


IjCsj  =  Ic(5)[G(s)*G*(s)3  =  I0(s)H(s) 


(11.56) 


The  normalized  component  Ij(s)  is  the  ratio  of  intensity  component 

•*  A 

I J  c S )  to  the  uniform  background  intensity  lj(0);  accordingly,  H(s)  is 
called  the  contrast  transfer  function. 

Equation  (11.56)  can  be  used  to  determine  the  resolution  limit 
for  incoherent  images.  Clearly,  unlimited  resolution  demands 


H(s)  >  0 

for  all  £ 

G(s)  >  0 


(11.57) 


If,  in  addition,  the  image  is  to  be  aberration- free  the  constraint 


G(s)  =  constant 


(II. 5S) 


is  also  required.  Of  course,  no  system  with  these  characteristics  has 
ever  been  derised.  T5*e  first  condition  can  never  be  realized  because 
object  details  separated  by  less  than  a  wavelength  do  not  excite  pro¬ 
pagating  waves.  It  is  unlikely  that  the  second  condition  will  hold 
over  more  than  a  small  angle  about  the  axis  of  the  imaging  system 
(paraxial  approximation)  because  of  path  length  differences.  More 
realistically,  we  assume  that  for  reasonably  good  imaging 

■ n  (fe) = "  (i) "  (t)  t!I-59) 

where 


n(x)  =  1  jxj  <  1 


=  0 


>  1  . 


(11.60) 
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Then 


HCS)  .  A  (£)  -  A  A  (j^j 

where 

A(x)  =  1  -  jxj  Jxj  <  1 

=0  SX1  >  1 


(11.61) 


(11.62) 


A  A 

These  idealized  forms  for  G(s)  and  ri(s)  are  shown  for  positive  spatial 
frequencies  in  Figure  II. 3-  We  see  that  all  spatial  frequency  com¬ 
ponents  of  the  object  boundary  intensity  lying  below  s  =  2sr  are  re¬ 
solved  in  the  image  although  the  contrast  near  the  limit  2sm  will  be 
very  poor  and  the  image  barely  discernible  against  the  uniform  back¬ 
ground. 

The  physical  interpretations  of  Equation  (11.54)  is  net  as  simple 
as  the  interpretation  of  Equation  (11.53).  Whereas  incoherent  imaging 
is  linear  in  intensity,  coherent  imaging  is  nonlinear  in  intensity. 
Consider  an  object  bounJary  field  Uq  containing  spatial  frequency 
components  up  to  s  =  s^.  The  spatial  frequency  spectrum  Uq(s)  may  then 
be  written 


U0(s)  =  Uq(s)  (11.63) 

where  Uq  represents  the  "strength”  of  each  component.  Since  the  zero 
frequency  spectral  component  Uq(0)  may  vanish,  Equation  (11.54)  cannot 
be  normalized  in  a  meaningful  way.  Assuming  again  the  ideal  form 
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(11.64) 


G(s)  «  G(0)n 


Equation  (11.54)  may  be  rewritten  as 


[(s)  =  %(s)n  (s)n  jG(0)}2 

■  (t)!*^"  ©] i5(o) 


sa>so 


(11.65) 


s  <  s~ 
m  0 


Relations  (11.65)  indicate  that  for  spatially  coherent  fields 

the  image  intensity  may  contain  spatial  frequency  coqxments  up  to 

2sr  when  s  >  s„  and  up  to  2s  when  s  <  srt.  However,  those  image 
u  mu-  m  my 

components  in  the  range  s  >  sQ  in  the  first  case  and  s  >  sq2i  the 
second  case  do  not  arise  from  object  components  lying  within  correspond¬ 
ing  ranges  of  the  object  spectrum,  hence  do  net  represent  resolution 
of  object  spectra  above  sQ  in  the  first  case,  above  s^in  the  second 
case.  The  spurious  image  spectra  is  simply  a  manifestation  of  the 
interference  which  occurs  as  the  coherent  field  propagates  from  object 
to  image.  Thus,  the  coherent  imaging  law  (11.54)  is  interpreted  to 
mean  that  all  spatial  frequency’  components  of  the  object  boundary  field 
lying  below  s  =  s  are  resolved  in  the  image. 

»3 

Comparing  the  form  of  6(s)  with  the  fora  of  H(s)  it  is  evident 
that  the  contrast  of  coherent  images  is  greater  than  that  of  in¬ 
coherent  images.  However,  because  of  the  interference mencsena 
alluded  to  above  coherent  images  suffer  from  edge  ringing,  edge  dis¬ 
placement,  speckling"",  etc.,  effects  which  are  not  present  in  in¬ 
coherent  images. 


(kt  the  basis  of  the  preceding  analysis  it  is  scrretlres  daises! 
that  the  resolution  of  object  detail  through  a  sell-corrected  insgisg 
system  operating  with  incoherent  radialire  ?s  twice  the  resolution  of 
object  detail  through  the  smr»  system  operating  with  coherent  radia¬ 
tion.  However.  this  sratesmt  recu ires  ssl  if  icaties.  >fcre  properly, 
if  aa  optical  system  operating  with  coherent  radiation  resolves  spatial 
frequency  cccponents  of  the  object  hrsrfsry  field  less  than  s_  «  s„, 
then  the  same  system  opera  ting  with  incoherent  radiation  will  resolve 
frequency  carecnsnts  of  the  object  feetsdary  invests ity  less  than 
£  =  2s„-  Therefore,  is  order  tist  resolution  be  determined  is  a  sesst- 
isgful  way  we  must  distinguish  betweei  rise  structure  of  the  object 
boundary  surface  C(r) ,  fine  structure  of  the  ccrxesrctci^  object 
boundary  field  U^Cr) ,  and  fine  structure  of  the  catTcspssssiisg  object 

_  y 

boundary  iatsasit)-  iU^Cr}-^.  In  zest  practical  situations  the  three 
say  be  considered  essentially  identical*  ftKeser.  when  operating  unde- 
diffraction  lisited  renditions  the  distinction  is  important. 

Assne,  for  sinplicity,  a  planar  test  object  (lying  is  the  xy 
plane)  constructed  such  that  sene  physical  property'  ?(x,y)  of  tie 
object  naterial-  for  example  the  blacsEiirg  of  a  photographic  emulsion, 
varies  over  the  reject  surface  in  a  ree-dlmessireal  cosinusoidal 
pattern.  Tas. 

?{x,y)  =  P(x)  =  ere  2xs„x  *  B  (11.66) 

where  s^  is  the  spatial  freresxy  of  the  cosinusoidal  variation  and  B 
is  the  unperturbed  (quiscest)  state  of  P.  Assure  further  that  with  ? 
nay  be  associated  sore  optical  characteristic  of  the  aaterial  such  that 
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when  the  object  is  uniformly  illuminated,  the  resulting  boundary  field 
UQ^y)  at  the  objecc  surface  will  be  proportional  to  P.  Typically, 
for  small  variations  of  P  about  the  nominal,  the  boundary  field  will 
assume  one  of  the  following  characteristic  forms.  Either 

U0(x,y)  =  2  +  cos  2tfsxx  (11.67) 


or 


UQ(x,y)  =  | cos  2ir$xx| 


(11.68) 


or 


UQ(x,y)  =  cos  2ttsxx  .  (11.69) 

Forms  (11.67)  and  (11.68)  represent  objects  which  do  not  introduce 
phase  variations  into  the  field  (non-phase  objects),  for  example  when 
the  object  is  a  photographic  transparency  or  a  thin  me. ,1  plate.  Form 
(11.69)  represents  objects  which  do  introduce  phase  variations  into  the 
field  (phase  objects),  for  example  a  suitably  designed  phase  grating. 

In  many  instances  the  object  will  be  of  type  (11.67).  Certainly,  this 
is  true  of  the  acoustic  targets  used  in  the  experiments  detailed  in 
Chapter  V. 

If  the  boundary  field  is  given  by  (11.67)  then  the  boundary 
intensity  I^(x)  is 
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IqOO  =  [|U0(x)|]2  =  (2+cos  2ttsxx)2 


(11.70) 


9  l 

j  •>  4  cos  2ttsxx  +  j  cos  4ttsxx 

Thus,  if  the  physical  boundary  has  a  spatial  frequency  component  s  , 
the  intensity  at  the  boundary  will  have  an  additional  component  2s  . 

A 

However,  this  spatial  harmonic  will  be  only  1/8  the  amplitude  of  the 
fundamental  component,  hence  can  be  neglected  in  determining  resolution. 
Therefore,  for  non-phase  objects  of  the  type  (11.67)  the  spectral  com¬ 
positions  of  object,  object  field,  and  object  intensity  are  substan¬ 
tially  the  same  and  to  resolve  one  is  to  resolve  the  others. 

If  the  boundary  field  is  given  by  either  (11.68)  or  (11.69)  then 
the  boundary  intensity  is 

Iq(x)  =  ~  (1+  cos  4ttsxx)  .  (11.71) 

Again  for  non-phase  objects  the  spectral  composition  of  intensity  and 
field  are  the  same.  In  contrast,  for  phase  objects  the  intensity 
contains  the  first  spatial  harmonic  2sx  plus  the  usual  uniform  back¬ 
ground,  whore  is  the  field  contains  only  the  fundamental  spatial  fre¬ 
quency  sx-  Thus  the  resolution  of  phase  objects  is  different  from  the 
resolution  of  non-phase  objects,  and  any  statement  comparing  coherent 
and  incoherent  resolution  must  include  this  fact. 

In  acoustic  imaging  we  shall  have  occasion  to  consider  orders  of 
the  field  higher  than  n  =  2.  In  particular,  consider  Uq.  If  the 
boundary  field  is  given  by  (11.67)  then 
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(11.72) 


[UQ(x)]4  =  +  38  cos  2ttsxx  + 

-P-  cos  4tts  x  +  2  cos  6tts  x  +  4  cos  Stts  x 

4  X  X  O  X 

If  the  boundary  field  is  given  by  either  (11.68)  or  (11.69)  then 

[Uq(x)]4  =  -•  +  ~  cos  4ttsxx  +  ^  cos  8rrsxx  .  (11.73) 

Again,  phase  objects  are  characterized  by  a  predominance  of  higher 

4 

order  harmonics  in  the  spatial  frequency  spectruii  of  UQ  whereas  for 
non-phase  objects  the  fundamental  is  still  predominate.  These  results 
are  compiled  in  Table  II. 1. 

We  can  now  compare  coherent  and  incoherent  resolution  directly 

in  terms  of  the  object  boundary  itself.  If  the  phase  of  the  field 

does  not  vary  across  the  object  boundary  surface  (non-phase  object) 

then,  a.*:  we  have  shown,  the  spectral  compositions  of  object,  object 

field,  and  object  intensity  are  substantially  the  same  and  incoherent 

resolution  is  indeed  twice  that  of  coherent  resolution.  Conversly,  if 

the  phase  of  the  field  is  not  constant  across  the  object  boundary 

surface,  (phase  object)  then  the  predominate  spectral  component  of  the 

object  intensity  occurs  at  twice  the  spatial  frequency  of  the  object 

field.  In  this  case  incoherent  resolution  of  object  detail  is  no 

greater  than  that  of  coherent  resolution  of  object  detail.  These 

results  are  compiled  in  Table  II. 2.  Note  however,  that  resolution  with 

coherent  radiation  is  never  greater  than  s  whereas  resolution  with 

incoherent  radiation  is  .iever  less  than  s  .  For  this  reason  incoherent 

m 

radiation  may  be  of  considerable  utility  in  acoustic  Imaging. 
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Predominate 


Component 


Object 

U0 

U0US 

cos  2nsx 

sx 

sx 

2sx 

2sx 

|cos  2irsx| 

2sx 

2sx 

2sx 

2sx 

2  +  cos  2t!sx 

sx 

sx 

sx 

sx 

TABLE  II. 2 

Maximum  Resolvable  Object  Spatial  Frequency  Component 

n=2  n=4 


Coherent 

Incoherent 

Coherent 

Incoherent 

cos  2ttsx 

sm 

sm 

m 

2s 

m 

(cos  2ttsx1 

sm/2 

V2 

2  s 
m 

2  +  cos  2ttsx 

sm 

2  s 

m 

sm 

4sm 

m 

6.  ACOUSTIC  IMAGING  WITH  FOURTH  ORDER  DETECTORS 

As  we  have  stated,  the  fundamental  observable  for  acoustic  radia¬ 
tion  is  instantaneous  field  amplitude;  other  observables  may  be  syn¬ 
thesized  electronically.  An  analysis  of  acoustic  imaging  based  upon 
field  intensity  as  the  synthesized  observable  has  just  been  completed. 
We  now  consider  another  observable 

D(r)  =  <([V(r,t)]4)>  .  (11.74) 

Interest  in  4th  order  detection  for  imaging  follows  from  the 
analysis  of  the  preceding  section  wherein  it  was  demonstrated  that 
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resolution  with  spatially  incoherent  radiation  is  often  greater  than 
with  spatially  coherent  radiation.  Mathematically,  improved  resolution 
arises  from  the  convolution  G(s)  *  G*(s)  occurring  in  the  spatial  fre¬ 
quency  transfer  function.  The  convolution  occurs  because  the  field  is 
spatially  incoherent  and  because  the  image  observable  is  intensity. 
Thus,  one  is  naturally  led  to  consider  imaging  with  higher  orders  of 
the  field,  say  the  4th  order.  We  now  inquire  into  the  optical 
properties,  particularly  resolution  of  an  imaging  system  when  the 
observable  is  <^[7(r,t)  ]4^>  rather  than  ^¥(r,t)]^. 

First,  the  time  average  of  the  4th  power  of  the  field  must  be 
calculated.  To  do  this  recall  that  for  quasi -monochromatic  Tadiation 

Y(r,t)  =  Re  [S(r,t)e2lTm]  .  (II. ?5) 


Thus 


nr.t) 


2 


(11.76) 


Substituting  (11.76)  into  <^[Y(r;,t)]4^>  ,  expanding,  and  supressing  the 
r  and  t  dependence  of  £(r,t)  yields 


<b(r,t)]4>  - 


16 


.  *4  -iSmvtX 

+  \S — ® -  \ 

16  / 


+  <  1  C5C%i4^t>  +  ^3e'i4^V  (11.77) 


+  <  3  2  *2^ 

X 
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New  consider  the  first  term  (11.77).  Expressing  the  time  depend 


ence  of  K  in  terms  of  temporal  Fourier  transforms  and  commuting  time 
and  frequency  integrations  we  obtain 


lim 

T-x» 


mi  icvpicv^icv^icv^ 


(ii.?8: 


^  /*  2iii(v1+v2+v3+v4+4v)t 

zr  J  e  dt 

dv.dv~dv.dv. 

12^4 

•  -T  J 

Carrying  Gut  the  time  integration  gives 


^4  i87Tvt\ 


16 


lim 


Jiff  i(v1)c(v2)i(v3)c(v4) 

-CD 

sin  2tt(v1+v2+v3+v4+4v)T 


2tt(v1+v2+v3+v4+4v)T 


dv.dv0dv,dv. 

12^4 


(H.79) 


For  quasi-monochrcmatic  radiation  (11.79)  is  readily  evaluated, 
sin  2t7(\>.+v_+v-,+v4+4v)T 

As  T  becomes  large  - - -  becomes  an  increasingly 

2h(v^+v2+v3+v4+4v)T 

narrow  spike  centered  at  v^+v2+v3+v4+4v  =  0.  Further,  recall  that  for 
quasi -monochromatic  radiation 


Hence,  in  the  limit  of  time  averages  over  durations  2T  long  compared 
with  the  mean  period  1/v  of  the  radiation,  no  combination  of  the  xn 
exists  for  which  both  the  sinx/x  term  and  all  F.(xn)  are  simultaneously 
non- zero.  Therefore, 


0 


(11.81) 


Similar  results  obtain  for  the  second,  third,  and  fourth  terms  in 
(11.77)  leaving  as  the  final  expression  for  the  fourth  order  acoustic 
image  observable 


<[nr,t)]^>=  |  <^C(r,t)32[CX(r,t)]2>  =  D(r)  .  (11.82) 

In  passing  note  that  by  applying  similar  techniques  we  may  show 

^(r.t)]2)  <5 (r , t) (r , t)^>  ,  (11.83) 

long  familiar  as  the  expression  for  the  intensity  of  purely  mono¬ 
chromatic  radiation,  and 


(p(T,t))Zn+1y  =  0  , 


(11.34) 


demonstrating  that  time- averaged  odd  order  acoustic  image  observables 
are  not  possible. 

Substituting  (11.82)  into  (11.34)  and,  as  before,  invoking 
ccRsautability  of  time  and  spate  integrations  obtain 
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DjClj)  =  f f f  j G  (Ij  G  (ij  “Is)  G  til  'I3)  G  (lj  - 


u' 


(11.85) 

<Jo !'— 1 » ^0  (—2  ’ r)  ^0  ^—3 f  t^0  (^4  *  t}/>  ^  rl^  r2^  r3^  r4 


4  ■ 


To  proceed  further  the  statistical  properties  o£  the  optical  field 
must  be  introduced.  The  situation  for  4th  order  detection  is  mere 
ccctplex  than  for  2nd  order  detection  because  of  the  need  in  the  former 
to  consider  4th  order  as  well  as  2nd  order  time  averages.  In  general, 
the  two  will  not  be  related;  i.e.,  <f^(r , , t) (r 2 ,tKQ(r3 , t) 5*  (1^ , t)> 
is  not  related  to  <^(r^  » t)  ^q(I?  » t))>  »  and  images  formed  according  to 
(11.85)  will  provide  additional  information  about  the  object  region, 
hopefully  in  the  form  of  greater  resolution  of  object  detail. 

If  the  field  is  spatially  incoherent  then  when  the  are  widely 
separated  the  time  evolutions  of  the  corresponding  £(r. ,t)  are  statis- 
tically  independent  and  the  fourth  order  time  average  may  be  decomposed 
into  products  of  second  order  time  averages  thus 


<S)di  *  t^i'C'(-j  ’  t)^o(-4c,t)  ^ol-i 1 1^/>  =  >t)/> 


(11.86) 


or  any  permutation  of  (11.86).  Reran  and  Parrent25  have  shc*n  that 
for  any  radiating  field 


<  aA 


meaning  that  as  the  separation  between  two  field  points  becomes  less 
than  ciA(a>l)  the  corresponding  time  evolutions  for  radiating  fields 
must  become  correlated*  The  time  average  <(?(r^,t)£  (r ^ , t)^>  may  be 
written  in  many  different  .'orms.  For  convenience  we  assume 

■  (II. 88) 

Equations  s.II.86)  and  (11.88)  are  suitable  expressions  for  the 
4th  order  time  average  only  when  at  least  one  of  the  r^  is  widely 
separated  from  the  other  three.  When  all  the  lie  within  seme  small 
distance  6A  of  each  ocher  (Bf'a, $>l,o>l)  it  rust  then  be  assumed  that 
the  decomposition  (11.86)  is  no  longer  valid.  In  this  case  we  take 

<|3(r1,tj4(l2,t)e0(r.,t);;(r4,t)>  - 

(11.89) 


Thus,  it  is  necessary  to  write  the  4th  order  time  average  in  different 
forms  depending  upon  whether  the  field  points  are  separated  or 
clustered  together. 

Using  Equations  (II. 86),  (II. 88),  and  (II. 89),  Equation  (II. 85) 
becomes 
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the  iscpisnatic  issuing  la*  for  spatially  Incoherent  acoustic  fields 
%&en  the  usage  observable  is  the  output  of  a  4th  Isa  detector.  If 
Dq(Tq)  is  statistically  independent  of  I^Cr-O  then  an  isage  foised 
according  to  (11.92}  say  exhibit  better  resolution  of  object  boundary 
detail  than  one  fenced  according  to  (11.59) . 

In  order  to  satisfy  siruir^necusly  conditions  (II. S3)  and  (II- S§) 
it  is  necessary  that  the  modulus  as  veil  as  the  phase  of  the  field 
amplitude-  be  renders  functions  of  tine.  To  see  this  write 
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aL,t)  =  jCCr^je1*^ 


(11-33) 


«hers,  in  geseial,  the  zzodalus  j£(r,t)}  3*3  fssse  e(r,t)  are  statisti¬ 
cally  independent  randan  functions  of  tine.  Then 


<jCli  ,t)5*(rjftf>  =  ^ICd^Ol^Xe 


iCo(ri,t)-d(r.,t)]> 


(11.54) 


and  it  is  clear  that  a  spatially  incoherent  field  in  the  sense  (II .§5) 


can  be  synthesized  sirply  by  ensuring 


<e  *  1  >  =  0 


rtf 


vhich  is  possible  if  o(£.t)  is  a  suitably  renders;  function  of  tine. 
However,  if  j£(r,t)]  is  deterministic  then  and  In  are  ret  indepen¬ 
dent  which  violates  condition  (II. S9) .  Therefore,  iu  order  that  4th 
order  detection  provide  additional  informatics  vr  desand  that  the 
nodules  and  phase  of  the  field  be  statistically  Independent  rszdaa 
runctix-jss  of  tire. 

Fields  of  this  type  are  called  fluctuating  fields.  Fluctuating 
fields  ray  be  synthesized  by  employing  large  arrays  suitably  drives  or 
they  say  occur  naturally  as  a  result  cf  turbulence  and  scattering  is 
the  a&diun  interposing  object  and  inage,  a  css ass  occurrence  is  tie 


To  test  the  cent  entice  that  increased  resolutics*  of  object 
boundary  detail  Is  obtained  vlth  4th  order  detection  take  the  Spatial 
Fourier  Transfers  of  (11.92)  obtaining 


*  »■  Ml ''.iPlilffilWlllMIW!!l!llll4 


Dj(s)  =  2*W  |[|G(r)|i]I0(^[  .  |[|G(r)hl0(s) 


A4X4  [G(r) 


:i)!2<t50(£,t)l2;  .  tG(r)!2<tE0C£,t)J^ 


♦  nWi|G(r)|4]D0(s) 


(11.96) 


As  before  the  higher  spatial  frequency  components  of  the  first  and 
second  terras  in  (11.96)  do  not  represent  increased  resolution  of  the 
object.  However,  the  third  terra  presents  interesting  posibilities. 
With  no  loss  of  generality  assume  G(s)  =  G(-s).  Then 


|C(r)j4  =  G(s)  *  G(s)  *  G(s)  *  G(s)  . 


(11.97) 


Computation  of  the  four- fold  convolution  yields 


|8(r)l4  " 


s3  -  2s  s2 
rax  2  mx  x , 


.  (Hs34-  2s  s2) 
\  3  ray  2  my  yf 


H  s3  +  2s  s^+8s2  s 
3  rax  6  mx  x  mx  x 


.  3  y  +  2  ^  g  2  \ 

U  my  6  my  y  ray  y  J 


which  is  plotted  normalized  in  Figure  II. 3.  Note  that  the  spatial 
spectrum  of  |G(r)|4  extends  over  a  range  four  times  broader  than  that 
of  G(£). 


2sm  <  sx,sy  <  2sra 


(11.98) 


lsxl'lsy!<  4Sb 
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Referring  to  Table  II. 1  it  is  apparent  that  for  non-phase  objects, 

spatial  frequency  Sq,  the  predominant  spatial  frequency  component  in 
4 

*J0(r)  is  also  Sq,  although  spatial  harmonics  through  the  third  order 

will  be  present  as  well.  Maximization  of  object  resolution  does  not 

4 

require  the  total  spatial  spectrum  of  Uq  be  imaged,  only  the  spectnmi 
of  the  object  boundary  itself.  Therefore,  if  the  object  boundary  com- 
ponent  Sq  in  UQ  is  transferred  to  the  image  and  is  not  lost  in  the  back¬ 
ground  noise  then  the  corresponding  object  boundary  detail  will  be 
resolved.  Thus,  in  principle  the  resolution  of  any  optical  system 
operating  with  4th  law  detectors  and  spatially  incoherent  fluctuating 
fields  may  be  greater  than  the  same  system  operating  with  coherent 
radiation. 

Since  only  the  |G(r)j4  term  in  Equation  (11.92)  provides  useful 
imaging  it  is  desirable  to  eliminate  the  other  two.  In  principle  this 
may  be  accomplished  simply  by  subtraction.  Even  so,  the  |G(r)|^  term 
cannot  be  too  small  or  else  it  will  be  buried  in  the  noise.  An  ap¬ 
proximate  worst  case  evaluation  of  the  relative  size  of  die  terms  can 
be  found  by  comparing  contributions  to  the  uniform  background. 
Straightforward  calculation  gives  the  ratio  of  contributions  to  the 
uniform  background  from  the  first  and  third  terms  of  (11.98)  as 
approximately 


(11.99) 
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where  s  =  s  =  1/yX ,  y  >  I.  It  should  be  possible  to  arrange  the 
edc  my 

field  statistics  (a  and  B)  so  that  (11.99)  is  oi  the  order  10  1  to  1 
for  typical  y's.  It  may  be  that  the  desired  effect  will  be  too  small 
to  be  readily  applicable  to  real  tine  imaging  systans.  However,  4th 
order  detection  could  be  useful  for  seismic  investigations,  such  as 
fault  location,  where  real  time  imaging  is  of  little  value. 

In  order  to  understand  physically  why  resolution  of  object 

boundary  detail  shot  Id  depend  upon  field  statistics  we  n&ist  distinguish 

between  resolution  and  total  quantity  of  information  transferred  frcsn 

27 

object  space  to  image  space.  The  two  are  not  necessarily  the  same  . 

The  general  imaging  process  involves  information  transfer  frcsn  an 
object  to  some  optical  representation  of  that  object  via  a  field  which 
is  not  completely  deterministic.  The  information  is  transmitted  as 
spatial  variations  of  the  modulus  and  phase  of  the  optical  field 
amplitude;  and  is  received,  in  the  case  of  real  time  visible  imaging, 
by  some  device  which  senses  only  time  averages  of  various  functions  of 
the  optical  field.  Because  of  the  time  averaging  characteristic  of  the 
receiver  the  phase  information  will  be  lost,  except  for  those  situa¬ 
tions  outlined  in  the  following  paragraph,  and  complete  reconstruction 
of  the  object  scene  is  impossible.  However,  resolution  of  a  boundary 
detail  in  the  detected  optical  image  requires  merely  that  sane  intensity 
pattern  in  the  image  can  be  isolated  and  uniquely  associated  with  a 
corresponding  boundary  detail  of  the  object.  In  principle,  the  phase 
of  the  field  is  of  no  consequence  except  to  enhance  contrast,  in  the 
case  of  spatially  coherent  fields.  Thus,  resolution  does  not  require 
the  entire  information  carried  by  the  optical  field. 
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For  convenience  we  have  classified  optical  fields  into  two 


rather  broad  categories  cased  upon  idealized  statistical  properties ; 
spatially  coherent  and  spatially  incoherent  fields,  clearly  special 
*•  ses  of  tie  '%?<■'-  \  partially  coherent  field.  These  categories 

are  useful  because  of  ti<eir  markedly  different  behavior  under  tiss 


averaging.  Either  class  of  field  contains  phase  and  amplitude  infor¬ 
mation.  However,  because  of  the  long  response  time  of  the  detection 
process  all  phase  information  will  be  averaged  out  unless  interference 
situations  (such  as  holography)  can  be  arranged.  Interference  with 
coherent  radiation  is  easily  arranged;  interference  with  incoherent 

7C 

radiation  entails  considerable  difficulty”  ’  .  Nevertheless,  since 

resolution  does  not  require  phase  information  we  cannot  conclude  that 
the  resolution  of  incoherent  systems  should  be  less  than  of  coherent 
systems. 

The  question  remains,  why  can  incoherent  resolution  be  greater 
than  coherent  resolution?  The  answer  is  that  whereas  the  thase  and 
amplitude  relation  between  object  points  coherently  illuminated  is 
independent  of  time,  the  phase  and  amplitude  relation  between  points 
incoherently  illuminated  is  a  rand  ran  function  of  time  and  position. 
Thus,  for  incoherent  radiation  the  detector  output  is  a  composite  or 
ensemble  over  many  of  the  possible  illuminated  configurations  of  the 
object  and  we  may  expect  that  time  averaging  will  provide  additional 
information,  possibly  in  the  form  of  enhanced  spatial  resolution  of 
object  detail.  For  coherent  radiation  spatial  coherence  severely 
restricts  the  number  of  illuminated  configurations  contained  in  any 
time  average;  hence  no  new  information  is  obtained.  Clearly,  the 
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averaging  tics  nist  be  long  compared  with  the  aean  tire  between  randcn 
fluctuations  of  the  field  parameters  or  else  the  detected  incoherent 
field  will  appear  as  if  it  were  a  coherent  field.  Ke  expect  thr_t  there 
is  a  relation  between  resolution  and  averaging  time  for  any  specified 
statistical  beh^  ior  of  the  field.  This  relation  could  be  of  use  in 
acoustic  imaging. 

7.  ISO  JOINT  RESOLbilGff 

The  resolution  ccsaparisons  given  in  Sections  5  and  6  are  mathe¬ 
matical  concepts  based  upon  transformation  of  an  object  spatial  fre¬ 
quency  spectiua  into  an  image  spatial  frequency  spectrin.  Resolution 
was  caspared  strictly  in  tenss  of  the  cutoff  frequency  of  the  transfer 
function  without  regard  for  its  shape.  Although  in  principle  the  eye 
should  be  able  to  resolve  to  the  lisvt  of  the  transfer  function,  be¬ 
cause  of  noise  and  because  of  the  phys  iologicai  processes  involved  in 
vision  there  is  no  guarantee  that  this  will  ever  be  possible.  Thus  a 
nor**  appropriate  test  of  an  optical  system  is  based  upon  the  ability  of 
the  eye  to  resolve  object  detail  in  the  h age.  It  is  obvious  that 
contrast,  hence  the  shape  of  the  transfer  .  unction,  must  play  an  im¬ 
portant  role  in  determining  visual  perception  of  an  image.  However, 
since  visual  perception  is  a  subjective  concept  it  is  not  amenable  to 
precise  calculation.  The  best  we  can  do  is  construct  mathematical 
models  which  approximate  the  process  by  which  the  eye  recognizes  object 
detail  as  being  distinct. 

One  such  criteria  is  the  Sparrow  two  point  resolution  condition^®. 
Consider  two  idealized  point  sources  symmetrically  disposed  about  the 
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origin  along  the  x*  axis  and  separated  by  distance  2a.  By  the  Sparrow* 
criterion,  the  images  of  two  such  points  will  be  considered  distinct 
(iesclved)  when  the  separation  2a  between  object  points  is  such  that 


S'V£>  1 
~^Tl 


(11.100) 


Ke  first  apply  (II. 100}  to  square  law*  detection  of  incoherent 
radiation.  For  simplicity  consider  line  sources  parallel  to  the  y* 
axis  rather  than  point  sources.  Then  we  have 


IjW  =  |G(x)j"  *  I0(x) 


i0(x*)  -  6(x'-a)  ♦  5(x'+a) 


and  desire  that  separation  for  which 


(11.101) 


.2.  , 

3  IT  ^ 

i  _  a 

- 1  “  ~1 

*  x=0  * 


\  J  |G(x-x')  j2  I0(x’)dx* 

A  J  X: 


=  f  5  I  (x')dx;  =  0 

{  x=0 


(11.102) 


Calculation  of  (II. 115)  is  facilitated  by  using  the  identity 


f(x)  =  f(s 


2d 


(11.103) 


whence  after  some  manipulation 
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£^252 


^-i.isa^  (-is)-  H(s)  ds  =  G 


Substituting 


fiCs)  ■  K)5  K) *  H)B  C  v) 


(11.105) 


into  (11.104)  yields  the  following  condition  upon  the  spacing  of  in¬ 
coherent  line  sources  if  they  are  to  he  resolved  by  square  Isk  detec¬ 


tors  according  to  the  Sparrow  criterion; 


?  ? 

(3-8r  a“s:jccs  4ras  +  3ras„  sin  4ras„  =  3 
£r  3  a  s 


(11.106) 


It  is  relatively  staple  to  find  a  numerical  solution  to  (11.106).  We 


obtain  the  ninisus  spacing  for  Sparrow  resolution  as 


„  _  .415  _  .83  _  2.61 

2“- V  ’W 


If  the  radiation  is  coherent  then 


(11.107) 


Ijto  -  |U0(x)  *  G(x)j2  , 


U.(x')  =  <$(x’-a)  ♦  6(x'+a) 


(11.108) 


Escicying  s?di3ds  siirilar  to  tkss  used  to  cotam  (1 1. 104}  givss 


3%  ! 


3x*  ! 


«x=0 


(' 


■?  \  $  t 

n  *  HE\  1 

^  i?y %  ^  'j ! 


x=0 


ji  *  Ji  i 

'TO*  5x  I 


r=0 


All 


:2Jias,„-2riss,K23is,j2g(s.) 


.  f.Jria(s's'-*e'2“2!s's'5lG*fs-s,Ms 


te 


]G  (s-s')d sds* 


(II. ICS) 


y ^  ^2-ias* 


+e-2sias,)(2sis,)5(s.) 


{ g2sia(s-s  ’ } +g- 2sia{s-s  * )  1  ^  (s_. ,  ^  (s_s ,  )dsds , 


=  0 


asstsing  G(s)  -  G(-s).  ^instituting 


G(s)  =  (II. 110) 

into  (11.109)  yields  the  following  condition  upon  the  spacing  of 
coherent  line  sources  if  they  are  to  be  resolved  by  square  law- 
detectors  according  to  the  Sparrow  criterion; 

4«s  a  cos  2as  a  +  f(2its  a)^-2}  sin  2«s  a  =  0  .  (11.131) 

n  m  in  n 


59 


and  the  %ssrraf  resolution  ccssditicsi  then  beocnes 


-2t  <  5  -  t 
t  <  s  <  2t 

(11.115) 


-t  <  s  <  0 
0  <  s  <  t. 


vhere  t  =  2s  into  (11.114)  yields  the  following  condition  upon  Us 
spacing  of  incoherent,  fluctuating  line  sources  if  they  are  to  be 
resolved  by  4th  order  detectors  according  to  the  Sparrow  criterion; 


60 


8. 


REALIZABILITY 


a.  Averaging  Time 

A  few  comments  of  practical  import  are  in  order.  A  fluc¬ 
tuating  field  will  exhibit  properties  characteristic  of  spatial  inco¬ 
herence  only  if  the  detector  averages  over  times  2T  much  longer  than 
the  mean  time  At  between  fluctuations  of  the  field.  For  a  signal  of 
bandwidth  Av 

2T  »  At  -  .  (11.118) 

Conversely,  real  time  imaging  demands  that  a  certain  minimum  number  N 
of  frames  be  formed  per  unit  time.  Thus 

2T<i  or  Av  »  N  .  (11.119) 

For  acoustic  imaging  at  the  standard  television  frame  rate  (N  -  30 
frames/sec)  the  required  bandwidth  is  Av  »  30  Hz,  probably  Av  =  300  Hz 
is  sufficient.  However,  for  this  bandwidth  the  monochromatic  length  is 
reduced  (5  meters  for  v  =  1  MHz)  hence  the  validity  of  the  mono¬ 
chromatic  assumption  must  be  reexamined  for  imaging  over  long  ranges. 

b.  Space  Bandwidth  Theorem 

In  the  preceeding  analysis  image  quality  was  presumed 
limited  only  by  the  Green's  function  for  transformation  of  wave  fields 
from  object  space  into  image  space.  Resolution  and  size  of  the 
detector  was  assumed  not  to  affect  image  quality.  Such  systems  are 
termed  diffraction  limited.  We  now  consider  the  constraints  which  must 
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be  Imposed  upon  the  detector  in  order  that  the  system  may  properly  be 
considered  as  diffraction  limited.  This  is  of  particular  importance  in 
acoustic  imaging  because  the  wavelength/ aperture  ratio  is  large. 

The  constraint  is  known  as  the  space  bandwidth  theorem,  a 
statement  of  the  conservation  of  che  space  bandwidth  product  (SBP) . 

The  SBP  of  any  optical  field  or  functional  of  the  field  is  simply  the 
product  of  its  maximum  spatial  extent  and  maximum  spatial  frequency  in 
any  plane  of  interest.  Clearly,  the  SBP  characteristic  of  any  detector 
must  be  related  to  the  SBP  of  the  optical  fields  which  are  to  be 
detected. 

The  conservation  of  SBP  between  detected  field  and  detector 
is  simply  a  consequence  of  the  scaling  properties  of  Fourier  Transform 
pairs;  i.e.,  if 

A(x)  =  A(s)  (11.120) 

then 

A(ax)=|A^|j  .  (11.121) 

If  A  extends  only  over  a  spatial  range  ix|  <  and  A  only  over  a 
spatial  frequency  range  |sj  <  sm  then  the  SBP  of  A(x)  is  xmsm.  If  'the 
detector  extends  over  a  spatial  range  jx|  <  x^  and  can  resolve  over 
spatial  frequency  range  |s|  <  s^  then,  in  principle,  if  x^s^  >  x^s^ 
the  field  can  be  accommodated  to  the  detector  by  use  of  an  ideal 
magnifier. 
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Heretofore,  it  was  assumed  that  the  optical  system  was 
space- invariant.  However,  in  order  to  calculate  the  field  of  view  of 
any  instrument *  it  is  obvious  that  a  space- invariant  formulation  is  not 
useful  since  such  assumptions  imply  unrestricted  field  of  view.  Thus 
we  must  reformulate  certain  parts  of  the  preceding  theory  to  account 
for  the  finite  field  of  view  of  the  instrument,  In  general  for  a 
space-variant  optical  system 

=  J  (Iq)  d2rQ.  (11.122) 

0 

For  each  position  r^,  G(rj.,i*Q)  may  be  written  as  a  function  or 
see  Figure  11.4  ,  but  this  functional  dependence  will  vary  with  Tq» 

We  denote  this  condition  be  written 


GtlpJb)  " 


(11.123) 


as  the  Green's  function  for  a  space-variant  optical  system. 

In  all  practical  situations,  especially  in  acoustic  imag¬ 
ing,  the  Green's  function  will  vary  continuously  across  the  object, 
becoming  less  symmetrical  in  r^-r^  near  the  extremes  of  the  field  of 
view.  Restricting  attention  to  those  systems  which  have  been  stepped 
down  so  that  the  Green’s  function  does  not  vary  across  the  field  of 
view  then 


(11.124) 
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where  denotes  the  edge  of  the  field  of  view.  Substituting  (11.124) 
into  (11.122)  gives  for  the  image 


[(-P  "  11  {r..')  r0 

n  L  \  r/  j 


(II.12S) 


which  is  in  the  same  form  as  the  unlimited  field  of  view  case  except 
that  the  original  object  field  £n(r„)  must  be  replaced  with  an  equiva- 

(jj)W]  • 


lent  space  limited  object  field 


The  relations  for  incoherent  imaging  then  become 


Ij(r)  =  |G(r) j2  * 


n  fe)  *0© 


(11.126) 


and 


Ij(s)  -  H(s) 


!0(s)* 


(sin  2nsxrfx^  /sin  Zhs^  ■ 


2ms. 


2tts. 


(11.127) 


1=5 


For  coherent  imaging 


Ix(r)  = 


G(r)* 


n(|)  uo^' 


and 


ijCD  - 


[.I  /sin  2"V&)  /sln  2"syrfvV 

Gt£)  h>®*  v-Hsrv  1“^  /. 

/sin  2ns  r  \  /sin  2ms  r  ' \ 


(11.128) 


(11.129) 


Consider  first  incoherent  imaging.  The  Green's  function  is 


such  that 
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ecrj-ip)  >  o 


iii-4-1  rg 


(11.130) 


=  0 


otherwise 


For  conventional  lenses  r  is  the  spot  size.  Frcm  (11.126)  it  is  then 

o 

clear  that  the  image  intensity  extends  over  an  area  (r  +r£X) (r  +r£y) 
in  the  image  plane.  From  (11.127)  it  follows  that  the  maximum  spatial 
frequency  component  in  the  image  intensity  is  2s^  since 

H(a)  =  A  U±-\  .  (11.131) 

\  / 

Thus,  the  space  bandwidth  product  of  the  image  intensity  is 


SBP  *  4s  s  (r  +Tr  1 (r  +r-  ) 
mx  my'  gx  fxJ  gy  iy' 


(11.132) 


and  in  order  that  the  image  be  diffraction  limited  the  detector  area 
XY  and  resolution  Sy  and  must  be  such  that 


XYS  S  ^  4s  s  (r  +r  ,  ) (r  +r-  ) 
x  y  mx  my*-  gx  fx' v  gy  fyJ 


(11.133) 


For  coherent  imaging  similar  considerations  lead  to 


OTSxSy  a  2smxV(V*r&UVV  (II'134) 

as  the  condition  upon  the  detector  in  order  that  the  imaging  be 
diffraction  limited. 

To  reduce  that  space  bandwidth  expressions  further  requires 
that  the  type  of  imaging  system  be  known.  For  example,  for  a  well 
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corrected  lens  of  large  aperture  r^  «  r^  and  to  good  approximation 
(11.134)  reduces  to 


XYS  S  >  2s  s  tc  rr 
x  y  mx  my  fx  fy 


(11.135) 


On  the  other  hand,  this  assumption  may  not  be  valid  for  holographic 
systems,  particularly  acoustic  systems,  he.  ce  the  full  expression 
(11.134)  must  be  retained.  In  Chapter  IV  the  space  bandwidth  theorem 
appropriate  for  acoustic  holograms  is  discussed. 

9.  IMAGING  KITH  NONLINEAR  ANTENNAS 

In  the  preceeding  we  have  considered  imaging  systems  in  which  the 
field  scattered  from  an  object  scene  is  imaged  by  passing  through  a 
suitable  focusing  aperture.  An  alternative  scheme  is  to  collect  the 
field  scattered  from  the  object  scene  by  a  highly  directional  antenna, 
scanning  the  object  scene  with  this  antenna  to  map  out  the  object  scene. 
Various  configurations  have  been  devised  for  producing  a  narrow  antenna 
receiving  pattern.  Of  special  interest  in  this  connection  are  the  non¬ 
linear  antennas,  in  particular  the  Ccvington-Drane  antenna^®* 

Because  instantaneous  pressure  amplitude  can  be  an  acoustic  observable 
the  Covington-Drane  antenna  is  readily  adapted  to  acoustic  imaging. 

Consider  two  receiving  antennas,  not  necessarily  identical,  which 
respond  to  instantaneous  field  amptlitude.  The  complex  instantaneous 
voltage  output  E  for  the  two  antennas  may  be  written 

El,2^-I,t:^  =  f  ^^Io,t^Gl,2^-I,^d  r0  (11.136) 

0 
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under  the  same  assumptions  of  temporal  behavior  as  before.  The  antenna 
patterns  ^  are  Measured  with  respect  to  a  common  coordinate  system 
and  Kq  is  the  object  boundary  field  distribution. 

If  the  instantaneous  outputs  of  the  two  antennas  are  multiplied 
and  a  long  time  average  is  taken  of  the  resulting  instantaneous  product 
we  obtain,  assuming  spatial  stationarity. 


^1  » t^2  “  J 

0 

•Gi^rio)G2(ii'^o)d  r0d  r0  * 


(11.137) 


As  we  have  already  shown  if  the  fields  are  illuminated  with  spatially 
incoherent  radiation  then 


(11.138) 


and  we  may  write  for  the  combined  antenna  outputs 

(—1  * t)E* (—1  * t-*/> =  /  I0(3b)GicIrIo)G2cIi'^b)d2ro  *  (n*139> 
0 
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For  convenience  degine  an  equivalent  spatial  frequency  transfei 
function  H^Cs)  for  the  composite  system  as 

H12(s)  =  GjCs^Cs)  .  (11.141) 

Then  it  is  clear  that  (11.140)  is  of  the  same  form  as  (11.53),  the 
corresponding  relation  for  conventional  incoherent  imaging  with  square 
law  detectors,  except  that  whereas  H(s)  is  an  autoconvolution  H-^Cs)  is 
a  crossconvolution.  This  is  an  important  distinction  because  one  has 
considerably  more  latitude  in  the  choice  of  forms  for  crossconvolutions 
than  for  auto  convolutions.  For  example,  the  Covington-Drane  antenna 
coasists  of  a  two  element  interferometer  and  a  uniform  array  of  equal 
length  placed  end  to  end.  If  a  uniform  array  as  large  as  the 
combination  has  a  response 

C-(s)  =  n(s/sm)  (11.142) 

then  the  response  of  the  Covington-Drane  system  can  be  shown  to  be 

S1?Cs)  =  H12(s)/H12(0)  =  n(s/2shP  (11.143) 

which  is  plotted  in  Figure  (I I. 3). 

Many  other  responses  can  be  synthesized  by  proper  choice  of  G^ 
and  G^.  The  important  point  here  is  to  note  that  since  amplitude  is  an 
acoustic  observable,  the  use  of  incoherent  acoustic  radiation  allows 
one  to  cascade  responses  in  a  way  which  provides  considerable 
flexibility  in  the  composite  system  response.  The  cascading  is  some¬ 
what  different  than  usually  found  in  linear  or  optical  systems  because 
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it  is  in  parallel  rather  than  series.  It  may  prove  fruitful  to  apply 
this  concept  to  imaging  with  more  conventional  devices  such  as  lenses. 

10.  SUMMARY 

The  most  serious  barrier  to  useful  acoustic  imaging  is  the  need 
for  physically  enormous  apertures.  It  is  difficult  to  secure  acoustic 
detectors  and  focusing  elements  in  the  large  sizes  required  for  these 
apertures.  Moreover,  the  ideal  acoustic  detector  should  be  in  a  form 
analogous  to  photographic  film.  Such  detectors  have  not  yet  been 
developed.  Thus,  maximizing  image  quality  and  resolution  for  a  given 
aperture  size  is  of  primary  importance  in  the  design  of  acoustic 
imaging  systems. 

In  this  chapter  we  have  considered  the  effect  of  field  statistics 
and  detector  characteristics  upon  image  resolution  and  contrast. 
Considerations  of  this  type  are  particularly  important  because  the 
choice  of  field  statistics  and  concomitant  detection  process  is  the  first 
step  in  the  design  of  an  imaging  system.  (Until  the  advent  of  the 
laser  there  was  usually  no  choice  for  optical  systems.)  Ke  have  shown 
that  for  acoustic  imaging  there  my  be  advantage  to  spatially  incoherent 
fields  over  i.oherent  fields  if  diffraction  limited  performance  for 
acoustic  lenses  and  mirrors  can  be  achieved.  Until  experimental 
investigation  of  incoherent  acoustic  imaging  systems  is  carried  out  it 

*  It  remains  to  be  determined  whether  achieving  diffraction 
limited  performance  is  primarily  a  design  or  primarly  3 
fabrication  problem. 
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1-  INnttXCIICK 


In  the  preceding  chapter  certain  basic  properties  of  a  genera! 
scalar  imaging  systea  vere  formulated;  aastic  isaglsg  «as  classified 


as  involving  acoustic  rather  than  ’risible  ”illS3 


I ef  of  the  cbjact 


scene.  Betore  investigating  bow*  a  practical  acoustic  13 


systais 


sight  be  devised,  cei 


peaslsrities  or  acoustic  ^csve  crupagsr-cn 


vhich  significantly  effect  the  quality  of  acoustic  rragss  should  be 
considered.  In  Chapters  IV  and  ¥  details  of  a  partimlar  acoustic 
imaging  sdsEse  sail  be  presented. 
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MathscaticalJy,  assmarg  that  scalar  aeve  theory  provides  valid 


representatits  of  tbs  fields. 


the  mnnent  neglecting  sttsssz- 


tica,  the  only  dirferetsce  between  acoustic  and  visible  electr-estegnetic 
radiatiesi  is  the  vsveleagth  X  =  c/f  icere  f  is  the  freqasry  and  c  is 


toe  velocity  or  isve 
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velocities  range  betveso  20u  r^sec  (ducrise  gas  at  standard  tascera- 
rc-re  and  pressure)  are  rt/sec  (quarts) ,  acoustic  frequencies 
between  20  Ha  sad  10  Git”1.  Has,  acoustic  vavelecsfctis  approach  those 
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coper  racge  of  tie  acoustic  spectrait.  hb&eerer. 


because  of  absorption  due  to  viscous  caspir^  and  other  Interact  acaisr 
effects,  long  range  acoustic  imping  of  undersea  objects  is  restricted 


to  radial 


cies  less  than  1  JSz,  {see  Section  8),  hsnoe  to 
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wavelengths  greater  than  1.5  mm;  i.e.,  approximately  3  *  10  longer 
than  visible  wavelengths. 

The  significance  of  the  restriction  to  relatively  long  wave¬ 
lengths  for  acoustic  imaging,  is  that,  in  many  situations  of  practical 
importance,  the  wavelength  cannot  be  assumed  small  with  respect  to 
characteristic  dimensions  (such  as  the  standard  deviation  of  surface 
irregularities)  of  the  object  surface  texture.  Thus,  most  objects  must 
be  considered  acoustically  "smooth”  even  though  they  are  visually 
"rough."  This  relation  is  the  cause  of  an  effect  known  to  illumi nation 
engineers  and  commercial  photographers  as  "high- lighting"  and  has  pro¬ 
found  influence  upon  the  quality  of  acoustic  images.  When  viewed  by 
acoustic  imaging  most  objects  will  appear  as  diffracted  highlights 
rather  than  extended  figures. 

All  object  surfaces  may  be  classified  according  to  the  following 
characteristic  forms: 

1.  smooth  (curved  or  planar) 

2.  multifaceted 

3.  rough. 

Description  of  the  surface  texture  as  either  rough  or  smooth  is  a  con¬ 
venient  categorization  of  the  texture  microstructure  relative  to  the 
wavelength  of  illumination  under  which  the  surface  is  viewed.  However, 
when  referring  specifically  to  the  concomitant  reflecting  properties  of 
a  given  surface  texture  it  is  more  customary  to  speak  of  either 
"diffuse"  or  "specular"  reflection,  the  former  to  be  associated  somehow 
with  rough  surfaces,  the  latter  with  smooth  surfaces. 
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TOE  CURVED,  SMOOTH  SURFACE 


% 

The  degree  to  which  an  object  will  be  highlighted  depends  on  the 
reflecting  properties  of  the  object  surface,  and  the  spatial  and  tempo¬ 
ral  structure  of  the  illumination.  A  simple  example  will  serve  to 
demonstrate  the  severity  of  highlighting  in  acoustic  imaging.  Consider 
the  problem  of  viewing  through  a  small  aperture  (width  2a  centered  at 
the  origin  y,  z  =  0)  a  specularly  reflecting  cylindrical  surface  (radius 
r  centered  at  y  =  0,  z  =  d)  illuminated  by  a  monochromatic  plane  wave 
traveling  in  the  z  direction.  An  elementary  analysis  assuming  geomet¬ 
rical  optics  will  show  that  the  extremal  rays  which  pass  through  the 
aperture  emanate  from  a  very  small  patch  (highlight)  on  the  object. 

Referring  to  Figure  III.l,  one  finds  the  following  relation  for 
the  extremal  rays. 

-  r  -  -  a  sin  <f>  -  d  cos  <j>  + 

(III.l) 

f~2  2  2  ! 

cos  (j)  V  a  +d  +r  -2ar  sin  $-2dr  cos  <J>  , 

which  is  valid  providing  r  »  X.  Assuming  the  paraxial  approximations 

sin  4>  cos  4>  «  1,  d  »  a  (III. 2) 

we  obtain 

$  =  |  |  ,  (III. 3) 


75 


Aperture  width 
Object  diameter 


which  is  also  the  ratio  of  apparent  diameter  2r&  to  actual  diameter  2r. 
Note  that  the  highlight  ratio  is  independent  of  object  curvature,  simply 
a  consequence  of  the  geometry.  Equation  (III. 3)  is  accurate  to  within 
-  201  even  when  d  £  a.  It  is  also  a  simple  matter  to  show  that  the 
power  density  P  received  from  a  highlight  is 

2 

P  cc*  (III. 4) 

d^ 

as  expected. 

It  is  apparent  from  (III .3)  that  the  highlight  ratio  will  be  very 
small  unless  the  aperture  is  enormous.  In  Chapter  V  we  present  the 
results  of  an  experiment  which  dramatically  illustrates  the  severity  of 
highlighting  in  acoustic  imaging  and  verifies  Equation  (III. 3). 

The  problem  imposed  by  highlighting  may  be  considered  from  a 
different  point  of  view.  In  principle,  even  though  objects  appear  as 
highlights,  object  size  can  nevertheless  be  determined  approximately 
from  highlight  size  provided  range  d  is  known;  range  can  be  determined 
by  range  gating  or  by  using  the  focusing  characteristics  of  holograms 
or  lenses.  For  example,  such  techniques  might  be  useful  for  measuring 
size  distributions  of  a  species  of  fish.  However,  the  range  over  which 
such  techniques  are  effective  is  severely  reduced  due  to  highlighting. 

According  to  the  Rayleigh  theory  the  angular  resolution  of  an 
apperture  of  width  2a  is 


(III. 5) 
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Thus,  beyond  a  distance  d* , 


(III. 6) 


it  will  not  be  possible  to  determine  object  diameters  less  than  r  by 
spatial  measurements  and  it  will  be  necessary  to  resort  to  intensity 
measurements  (Equation  I I I. 4).  Conversely,  if  the  entire  object  disk 
were  visible  the  corresponding  maximum  distance  of  resolution  would  be 


d  = 


2ra 

.61X 


(III. 7) 


and 


d’  1  /761X 
d  “  2  V  r 


(III. 8) 


From  (III. 8)  it  is  clear  that  the  effective  spatial  resolution  of  high¬ 
lighted  objects  is  much  less  than  the  Rayleigh  limit. 

4.  THE  PLANAR  .SMOOTH  SURFACE 

From  the  preceding  discussion  it  is  clear  that  satisfactory  acous¬ 
tic  imaging  of  curved  surfaces  is  difficult  to  obtain.  As  we  shall  now 
show  planar  objects  present  somewhat  different  although  equally  diffi¬ 
cult  inaging  problems.  Whereas  at  least  to  first  order,  the  highlight¬ 
ing  of  a  curved  surface  is  independent  of  its  transverse  position  rela¬ 
tive  to  the  axis  of  aperture,  the  transverse  position  of  a  planar  object 
is  critical.  Consider  the  situation  depicted  in  Figure  III. 2  showing  a 
planar  object  (width  2b)  parallel  to  Hut  offset  f  from  an  aperture 
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Highlight  of  a  Planar  Surface 


(width  2a)  and  illuminated  with  a  monochromatic  plane  wave  traveling  in 
the  z  direction.  If  f  +  b  <  a  then  the  entire  extended  object  will  be 
imaged.  If  f  -  b  >  a  then  only  the  field  diffracted  by  the  object  edges 
will  intercept  the  aperture,  hence  the  object  will  appear  as  an  outline 
rather  than  a  solid;  i.e. ,  the  image  looks  like  the  gradient  of  the 
object. 

This  can  be  seen  analytically.  From  Appendix  C  the  field  reflec¬ 
ted  by  the  object  in  the  plane  of  the  imaging  aperture  can  be  written 


(III. 9) 


where 


(III. 10) 


and  is  the  Fourier  Transform  of  the  field  in  the  object  plane.  If 
b  »  X  then  it  is  permissible  to  take 


and  (I I I. 9)  becomes 


-  CO 


Expanding  the  exponential  yields 


(III. 11) 


(III. 12) 
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where  C  is  a  corxstant.  Thus,  the  only  contributions  to  the  image  of  an 
offset  planar  object  cane  from  the  space  varying  parts  of  the  object; 
i.e. ,  the  edges. 

S.  MULTIFACETED  AND  ROUGH  SURFACES 

The  multifaceted  surface  is  a  randan  distribution  of  curved 
facets  disposed  about  some  smooth  median  surface.  For  example,  the 
barnacle  encrusted  hull  of  a  ship  might  appear  to  be  multifaceted  for 
long  acoustic  wavelengths.  An  acoustic  image  of  such  a  surface  will  be 
a  collection  of  highlights,  one  for  each  facet.  As  one  recedes  from  a 
multifaceted  surface  the  highlights  will  tend  to  merge  into  a  solid 
form.  As 
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(III. 14) 
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where  1  is  the  average  spacing  between  facets,  the  multifaceted 
surface  becomes  a  rough  surface. 

6.  ILLUMINATION 

Of  course  highlighting  would  present  considerably  less  problem 
if  diffuse  acoustic  illumination  sources  were  available.  Unfortunately, 
typical  acoustic  transducers  tend  to  be  point  or  collimated  sources 
rather  than  diffuse  sources.  In  visible  optics  such  common  sources  of 
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illumination  as  light  bulbs  tend  to  be  point  sources.  However,  there 
is  usually  enough  background  diffused  illumination  to  avoid  serious 
highlight  problems. 

fine  is  thus  led  to  ask  whether  there  is  sufficient  acoustic  back¬ 
ground  in  the  se3  to  render  a  similar  function  for  undersea  acoustic 

imaging.  Figure  I II. 3  shows  measurements  of  average  acoustic  back- 

32 

ground  noise  in  the  sea  as  a  function  of  frequency  .  The  two  non- 
biological  curves  represent  upper  and  lower  bounds  for  noise  due  to 
mechanical  sources  such  as  wave  motion,  rain,  thermal  agitation,  etc. 
Noise  due  to  biological  or  man-made  sources  is  considerably  stronger 
and  is  typically  found  in  the  range  indicated  near  the  top  of  the 
figure. 

The  dashed  line  dividing  the  figure  is  the  threshold  level  for 

piezoelectric  detection  of  acoustic  radiation;  this  threshold  is  of 

significance  because  piezoelectric  arrays  will  probably  play  an 

14 

important  part  in  any  future  undersea  acoustic  imaging  work  . 
Piezoelectric  detectors  are  attractive  because  of  the  very  low 
threshold  and  because  these  devices  can  withstand  the  large  hydro¬ 
static  pressures  which  exist  even  at  moderate  depths. 

It  appears  from  these  preliminary  data  that  noise  due  to 
mechanical  sources  will  not  be  useful  but  that  the  noise  due  to 
biological  sources  warrants  further  investigation.  For  example,  the 
acoustic  noise  above  shrimp  beds  is  quite  large  and  animals  such  as 
the  dolphin  are  known  to  be  the  source  of  considerable  sonic  activity. 
Simple  man-made  acoustic  sources  might  include  patterns  of  small 
explosive  charges  or  perhaps  streams  of  resonant  air  bubbles. 
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BIOLOGICAL  AND/OR  MAN-MADE 


7 


9 


Figure  I I I. 3.  Average  Background  1 


In  general,  the  problem  of  securing  sufficiently  diffuse 

illumination  appears  most  tractable  if  incoherent  rather  than  coherent 

acoustic  radiation  is  used  for  illumination.  First,  although  it  is 

possible  to  obtain  diffuse  illumination  free  coherent  sources,  the 

design  and  fabrication  of  the  necessary  diffusers  and  extended  coherent 

sources  for  acoustic  radiation  present  considerable  difficulty.  On 

the  other  hand,  as  explained  above  it  may  be  relatively  simple  to 

obtain  extended,  incoherent  sources  from  which  diffuse  acoustic 

illumination  can  be  generated.  Second,  of  necessity  incoherent 

17 

radiation  East  diffuse  as  it  propagates-1  . 

7.  ACOUSTIC  ATTENUATION  IK  KATER 

Neglecting  boundary  effects  due  to  the  ocean  bottom  and  surface, 
the  primary  source  of  acoustic  attenuation  in  water  is  viscous  losses 
and  inhomogeneities  such  as  sediment  and  air  bubbles.  For  pure  sea 
water  the  relation  between  loss  and  frequency  is  shown  in  Figure  III. 4 
for  various  levels  of  salinity'5^.  On  the  basis  of  this  data  it  is 
concluded  that  long  rage  acoustic  imaging  c-f  underwater  objects  will 
be  restricted  to  frequencies  no  greater  than  1  Miz;  i.e»,  to 
attenuation  less  than  .5  db/ra  in  pure  sea  water. 

In  addition  to  the  intrinsic  absorption  of  water  one  must  also 
consider  the  effect  of  dispersion  of  foreign  material  in  the  water. 

Such  dispersions  can  significantly  alter  the  attenuation  and  spatial 
frequency  transfer  function  of  the  water  medium.  For  example,  resonant 
air  bubbles  may  have  an  effective  scattering  cross  section  several 
orders  of  magnitude  larger  than  geometrical  diameter.  Thus,  even  a 
relative).)*  tenuous  dispersion  of  microbubbies,  such  as  might  exist  in 


the  wake  of  a  ship  many  hours  after  passage,  could  significantly 
degrade  the  quality  of  acoustic  images. 

8.  SUMMARY 

We  have  shown  that  highlighting  will  be  a  difficult,  perhaps 
impossible  problem  to  overcome  in  undersea  acoustic  imaging.  It 
appears  that  approaches  to  the  problem  in  terms  of  incoherent  acoustic 
illumination  are  the  most  promising.  On  the  basis  of  attenuation 
undersea  acoustic  imag.ng  will  probably  be  limited  to  frequencies  less 
than  1  MHz,  hence  to  resolutions  exceeding  1.5  mm. 
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CHAPTER  IV 


HOLOGRAPHY 


1-  INTRODUCTION 

In  Chapter  II  we  presented  a  general  formulation  of  scalar  imag¬ 
ing  emphasizing  acoustic  applications.  In  Chapter  III  we  presented 
certain  of  the  special  problems  which  arise  when  acoustic  fields  are 
used  to  image  objects  inmersed  in  a  liquid  medium  such  as  the  sea.  In 
this  chapter  we  turn  to  a  specific  imaging  technique  applicable  to 

acoustic  imaging*,  a  two  step  process  invented  by  Dennis  Gabor  and  called 

* 

wave  front  reconstruction,  the  hologram  process  or,  more  usually, 
holography.  We  will  concentrate  on  a  formulation  of  the  wavefront  re¬ 
construction  process  valid  for  large  diffraction  angles  and  including 
features  peculiar  to  the  acoustic  case. 

It  is  not  intended  that  this  be  an  exhaustive  account  of  holo¬ 
graphy;  only  those  features  of  immediate  application  to  acoustic  imag¬ 
ing  are  considered.  For  compact  presentations  of  the  entire  range  of 
holographic  research  the  reader  is  referred  especially  to  an  article  by 
Leith  and  Upatnieks28  and  a  book  by  DeVelis  and  Reynolds5**.  Additional 
reference  to  the  current  research  will  be  cited  in  the  text  as  appro¬ 
priate. 

Holographic  imaging  is  basically  a  simple  process  and  the 
rudiments  are,  by  now,  familiar  to  most.  An  arrangement  typical  of 
those  for  recording  holograms  with  light  is  shown  in  Figure  IV. 1. 


Hologram  was  coined  by  Professor  Gabor  from  the  Greek  work  HOLOS 
meaning  whole  or  entire.  It  is  appropriate  since  a  hologram 
records  both  amplitude  and  phase  of  wave  fields. 
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HOLOGRAM  RECORDING 


HOLOGRAM  READING 


Figure  IV.  1.  Typical  Arrangement  for  Wavefront  Reconstruction 
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Briefly,  a  hologram  is  recorded  by  adding  the  field  scattered  from  a 
coherently  illuminated  three  dimensional  object  scene  to  another 
coherent  field  (the  reference)  and  allowing  their  sum,  a  complicated 
standing  wave  or  interference  pattern,  to  expose  a  photographic 
emuision.  Lenses  are  not  required.  A  single  spatially  coherent  source 
provides  both  reference  and  object  illumination  so  chat  the  two  will 
interfere.  Usually,  the  reference  is  planar  or  spherical,  but  there 
are  situations  in  which  it  is  diffuse  or  may  even  be  part  of  the  object 
scene.  The  transparency  (called  the  hologram)  which  appears  upon 
development  of  the  exposed  emulsion  contains  a  mapping  of  the  object- 
reference  interference  pattern.  A  complete  record  of  the  object  field 
amplitude  and  phase  is  stored  in  this  pattern. 

Three  dimensional  images  which  are  faithful  replicas  of  the 
original  object  scene  are  obtained  from  the  hologram  by  illuminating  it 
with  a  duplicate  of  the  reference.  Usually  two  images  of  the  object 
scene  will  be  formed;  one  real,  the  other  virtual.  The  virtual  image 
corresponds  to  exact  replication  of  the  original  object  field,  hence 
the  term  wavefront  reconstruction,  and  is  viewed  by  looking  through 
the  hologram  as  if  it  were  a  window.  The  real  image  may  be  observed  as 
an  object  suspended  in  mid- space  between  observer  and  hologram  or  it 
can  be  projected  directly  onto  a  viewing  screen. 

Although  light  was  assumed  for  purposes  of  illustration,  the 
procedures  outlined  above  are  applicable  to  non-visible  optical  fields 
as  well,  for  example,  to  acoustic  or  microwave  fields.  Strict  preserva¬ 
tion  of  the  analogy  would  demand  that  acoustic  holograms  reconstruct 
acoustic  fields,  microwave  holograms  reconstruct  microwave  fields. 
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etc.^.  However,  holograms  o£  invisible  fields  which  reconstruct 
geometrically  similar  visible  fields  are  more  useful,  therefore 
discussion  is  restricted  to  holograms  of  this  type.  For  example,  by 
acoustic  holography  we  shall  mean  visualization  of  acoustic  fields  by 
converting  object-reference  acoustic  interference  patterns  into 
geometrically  similar  photographic  transparencies  which,  in  turn, 
serve  as  holograms  to  reconstruct  visible  fields  geometrically  similar 
to  the  original  acoustic  fields.  In  this  way  we  "see"  acoustically  by 
holography. 

The  motivation  for  holography  in  acoustics  is  not  the  same  as  in 
visible  optics.  In  visible  optics  the  need  is  to  record  information 
of  a  three-dimensional  nature  on  a  two-dimensional  photographic 
emulsion  in  a  form  suitable  for  easy  recovery  whenever  desired.  In 
this  application  holography  is  not  a  real  time  imaging  process;  lenses, 
photoelectric  devices,  and  the  eye  adequately  serve  that  purpose. 

Conversely,  in  acoustics  the  need  is  for  real  time  imaging  using 
acoustic  fields  as  the  "illumination”  but  with  the  resulting  acoustic 
images  somehow  rendered  visible  so  that  we  may  "see"  in  otherwise 
opaque  media.  Simple  imaging  devices  such  as  lenses  of  adequate 
quality  are  not  readily  available.  More  serious,  as  yet  there  are  no 
practical  acoustic  sensitive  materials  which  are  analogous  to  the 
retina.  Holography  attains  viability  in  acoustics  as  an  imaging  tech¬ 
nique  which  does  not  require  lenses,  hence  at  least  eliminates  the  first 
stumbling  block.  Moreover,  in  certain  limited  applications  real  time 
acoustic  imaging  by  holography  has  been  demonstrated0  . 
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2. 


THE  HOLOGRAM  EIGENVALUE  EQUATION 


The  remarkable  effects  attributed  to  holography  are  well  known, 
probably  the  most  dramatic  being  the  generation  of  virtual  inages  of 
diffusely  reflecting  or  diffusely  illuminated  three  dimensional  scenes 
from  two  dimensional  photographic  transparencies.  Upon  viewing  these 
images  with  the  unaided  =yc  one  observes  what  appear  to  be  the 
original  scenes  in  total  three  dimensional  detail  including  fore¬ 
shortening.  parallax,  and  depth  of  focus  as  well  as  the  stereoscopic 
sense  of  depth.  Holographically  generated  scenes  are  thus 
indistinguishable  from  original  scenes. 

This  phenomena  is  explained  in  the  following  manner.  The  holo¬ 
graphic  scene  cannot  be  distinguished  from  the  original  scene  because 
the  hologram,  when  properly  illuminated,  diffracts  an  exact  duplicate, 
both  in  amplitude  and  phase,  of  the  complex  wave  field  scattered  from 
the  original  scene.  We  may  say  that  the  hologram  process  (recording, 
developing,  reconstruction,  etc.)  operates  on  a  field  in  such  a  way  as 
to  recreate  this  field  at  another  time  and  place. 

This  relation  admits  of  a  simple  mathematical  representation 
which,  because  of  its  elegance,  has  appeal.  Denoting  the  field 
originally  diffracted  from  the  scene  by  &Q(r,t)  and  the  process 
wherein  the  hologram  P(r)  is  recorded  by  an  operator  (r , t)  which 

operates  upon  holographic  recording  of  a  field  may  then  be  repre¬ 
sented  as 

tfl(r,t)'4>0(r,t)  =  P(r)  ,  (IV. 1) 

the  first  step  in  holographic  imaging. 
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represents  some  general  time  averaging  process  of  the  form  (II. 2) 
by  which  the  holographic  record  is  formed.  We  shall  assume  it  reduces 
to  a  time  average  of  the  form  (I I. 3). 

Similarly,  denoting  the  reconstruction  process,  whatever  it  may 
entail,  by  an  operator  #2(r,t2)  acti-n2  on  the  hologram  P(r)  we  may 
write  for  the  second  step  in  holographic  imaging 

S2(r,t2)P(r)  =  «^(r,t2)  (IV.2) 

where  i|>j(r,t2)  is  the  field  diffracted  by  the  hologram. 

Combining  and  and  defining  the  hologram  operator  ff(r,tj,t2) 
as 

3(L,tvt2)  =  ^(r^p^Cr.t^  (IV.  3) 

we  obtain  for  the  entire  holographic  imaging  process 

^(l»tl»t2)^0(r,t1)  =  ^yfr,^)  .  (IV.  4) 

For  perfect  reconstruction 

^(r.t,)  =  h*0(r,t2)  (IV.  5) 

where  h  is  a  complex  constant.  Thus  for  perfect  reconstruction  (IV. 4) 
becomes 


BQ  -  hi 


(IV.  6) 
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Note  that  the  coordinate  system  is  fixed  with  respect  to  the  hologram 
because  reconstruction  should  not  depend  on  translations  of  the  holo¬ 
gram  after  recording. 

Condition  (IV. 5)  is,  in  a  certain  sense,  unduly  restrictive  since 
the  field  \p  is  usually  not  directly  observable.  When  only  functionals 
D  of  the  field  are  observable  then  as  long  as  the  original  and  recon¬ 
structed  observables  are  identical  the  reconstruction  will  appear  to  be 
perfect.  In  this  case  (IV. 5)  may  be  replaced  with  the  less  stringent 
condition 

Dj(r,t2)  =  h'DgCr,^)  .  (IV.7) 

Equation  (IV.  6)  is  immediately  recognized  as  a  form  occurring 
frequently  in  mathematical  representations  of  physical  processes,  for 
example,  in  quantum  mechanics36.  For  a  given  operator  ff(r,t^,t0)  only 
certain  functions  ♦ndfc  ,t),  the  eigenfunctions,  will  satisfy  this  equa¬ 
tion  and  only  for  special  corresponding  values  h^  of  h,  the  eigenvalues. 
In  general,  to  each  eigenvalue  correspond  several  eigenfunctions,  the 
number  being  the  degeneracy'  d^  of  the  eigenvalue  hn.  The  eigenvalues 
may  be  descrete,  continuous  over  a  range,  or  a  combination  of  the  two. 

Equations  (IV. 6]  and  (IV.7)  may  be  regarded  as  describing  a 
general  two  step  process  in  which  information  is  permanently  recorded 
at  time  tp  not  necessarily  in  recognizable  form,  to  be  recovered  later 
at  time  in  original  form.  Unless  5(r,t^,t2)  has  special  properties, 
an  arbitrary  field  will  not  be  perfectly  reconstructed.  Indeed,  even 
the  best  choice  of  H(r,t^,t,)  probably  not  guarantee  perfect 
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reconstruction  for  all  possible  object  fields,  but  only  some  (hopefully 
large)  class  of  them. 

At  this  point,  the  description  of  holography  as  embodied  in 
Equations  (IV. 6)  and  (IV.7)  is  purely  phenomenological;  it  describes, 
approximately  at  least,  the  central  feature  of  holography  but  it  does 
not  necessarily  contain  any  essential  physics  of  the  process.  Before 
concluding  that  (IV. 6)  or  (IV.7)  has  any  deeper  significance  or 
utility,  the  processes  operating  in  holography  must  be  studied  in 
detail  in  order  to  obtain  explicit  expressions  for  the  hologram  opera¬ 
tor.  Then  the  implications  of  the  eigenvalue  formalism  and  the  degree 
to  which  this  formalism  is  valid  and  useful  can  be  assessed  in  greater 
depth. 

3.  UNIQUENESS  OF  HOLOGRAPHIC  RECORDINGS 

Consider  the  general  image  recording  situation  depicted  in 
Figure  II. 1  where  the  image  space  I  is  composed  of  two  parts,  the 
recording  or  detecting  medium  P  (photosensitive  emulsions,  alkali 
halide  crystals,  etc.  for  visible  radiation)  and  the  image  forming 
device  F  (lenses,  prisms,  the  eye,  etc.).  The  simplest  imaginable 
situation  would  then  be  when  I  is  composed  entirely  of  recording 
material  with  no  auxiliary  processing  apparatus.  Holography  is  co¬ 
herent  imaging  with  precisely  such  elementary  configurations  of  the 
imaging  region. 

Assume  object  region  0  is  illuminated  with  monochromatic,  spa¬ 
tially  coherent  radiation  p(r.t)  =  ^(r)e^''1— ^  amanating  from  luminous 
source  region  l,  £ (r)  then  obeying  the  Helmholtz  Equation  (11.17), 
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Elliptic  equations,  such  as  the  Helmholtz  equation,  have  unique  and 

stable  solutions  for  either  Dirichlet  or  Neumann  boundary  conditions 

37  38 

over  closed  boundaries  *  .  That  is,  when  the  modulus  and  argument  of 

£(r)  or  -gpjC(r)  is  specified  at  all  points  on  sane  boundary  enclosing  a 
source-free  region,  the  solution  5(r)  to  the  elliptic  equation  every¬ 
where  inside  that  boundary  is  unique.  Referring  to  Figure  IV. 2,  con¬ 
sider  the  region  M>  bounded  by  sane  portion  S  of  the  surface  of  P,  the 
plane  q  extending  fran  the  extremities  of  S  to  the  sphere  at  infinity, 
and  the  half  sphere  R>  at  infinity.  R  shall  be  known  as  the  exterior 
region  and  the  field  inside  this  region  as  the  exterior  field  ^q>. 
Accordingly,  if  S  is  of  size  sufficient  to  intercept  all  of  the  field 

scattered  from  the  object  =  0  on  q),  and  taking  £n(r)  «*•  U,  then 

a  r  -  R 

the  boundary  fields  C^)  or  jjCCr,.)  will  bear  unique  relation  to  ?0>. 

A  similar  relation  does  not  hold  for  the  interior  field 
inside  the  interior  region  bounded  by  the  surface  S,  the  plane  q, 
and  the  half  sphere  R<  at  infinity  beca  use  the  source-free  assumption 
is  violated.  Physically,  the  lack  of  uniqueness  for  ^is  clear.  In 
general,  there  will  be  evanescent  as  well  as  propagating  fields 
associated  with  Cq<;  the  former  cannot  be  uniquely  determined  fran 
C^O-g).  Second,  only  thosi  regions  of  which  radiate  toward  q  will 
contribute  to  ^(tg) . 

Thus,  regardless  of  the  size  of  P  it  is  necessary  to  partition 
the  scattered  field  £q  into  an  interior  portion  ^  and  an  exterior 
portion  Eq>  according  to  the  prescription 

5q(I)  =  Co<Cl)[l  -  u(r  -  rg)]  «■  $0>(r)u(r  -  tg)  (IV. 8) 

where 
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Figure  IV* 2. 

A  General  Boundary  Configuration  for  Holographic 
Imaging,  Recording 
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fIV.9) 


Next,  suppose  that  exposure  to  radiating  irssxes  irreversible 

V 

alteration  of  sene  physical  property  of  tbs  recording  material  inside  I, 
thereby  imbedding  in  the  caterial  a  permanent  record  of  the  field. 

Fraa  the  arguments  presented  in  the  proceeding  paragi3x&-  if  rats 
record  retains  both  the  amplitude  and  phase  of  SyjClg)  or  ra^n 

it  will  be  a  unique  and  complete  represen  tavicr.  o’  w ,  bur  only  a 
partial  representation  of 

The  uniqueness  properties  given  above  for  recorded  optical  in¬ 
formation  can  be  summarized  as  a  statement  csreiahg  the  resonstrsa  or 
optical  instruments .  To  wit,  resolutiesi  is  5 imply  a  consequence  of  the 
uniqueness  properties  of  elliptic  equations-  If  the  boundary  S  is  not 
large  enough  to  intercept  all  -*f  the  field  £q  scattered  free  the  tfjject 
scene  then,  by  the  uniqueness  tfceorea,  SgCtg)  is  to  lesser  uniquely 
related  to  ^(r) .  Further,  as  the  boundary  S  decreases  in  site  the 
laxsber  of  possible  configurations  of  vfe)  which  coincide  with  a  gives 


Cjj(rg)  Gr  j^(rg)  increases-  Bence,  as  boundary  sire  decreases, 
uncertainty  in  the  relation  between  object  and  image  increases;  i.e. , 


resolution  is  reduced.  This  agrees  with  intuitive  notices  and 
quantitative  theories  governing  optical  resolution. 

As  object  fine  structure  becomes  mere  intricate  ^  will  be 
scattered  through  correspondingly  greater  angles  thereby  requiring,  a 
larger  surface  for  complete  interception  ci  the  scattered  field  and 
total  resolution  of  the  object.  However,  when  the  -object  structure 
contains  spatial  frequency  components  exceeding  I/lerea  a  very  large 


aperture  may  not  provide  complete  resolution  because  the  fields 

scattered  in  response  to  these  very  high  spatial  frequencies  will  be 

evanescent  rather  than  propagating.  Thus,  there  will  be  irreversible 

attenuation  of  the  high  spatial  fr  'quency  information  transmitted  from 

object  to  recording  surface.  Depending  upon  the  sensitivity  of  this 

surface  the  information  may  be  completely  lost  or,  at  the  very  least, 

severely  attenuated.  Because  of  the  relatively  long  wavelengths 

associated  with  acoustic  radiation  this  effect  is  of  considerable 

importance  in  acoustic  imaging.  Other  investigators  have  discussed  the 

problem  of  recovering  high  frequency  details  in  microwave  imaging,  and 

39 

the  partition  of  scattered  fields  into  interior  and  exterior  portions  . 
Even  if  evanescent  propagation  is  absent  the  recording  medium  will  not 
respond  equally  to  high  and  low  spatial  frequencies,  nor  to  large  and 
small  field  amplitude,  nor  to  all  directions  of  propagation.  These 
problems  are  also  of  import  ance  in  acoustic  holography. 

A  general  observation  on  the  nature  r'f  the  reconstruction 
process  will  facilitate  proper  formulation  of  the  eigenvalue  equation. 
Consider  the  reconstruction,  Figure  Tv.3,  from  the  recording.  Figure 
IV. 2.  The  hologram  P,  illuminated  by  coherent  source  C,  will  in¬ 
variably  be  one  of  two  possible  types;  either  predominantly  reflecting 
or  predominantly  transmitting.  For  reflection  holograms  the  recon¬ 
structed  fields  will  lie  in  the  interior  region  NP.,  for  transmission 
holograms  the  reconstructed  fields  will  lie  in  the  exterior  region  M1  . 
Thus  both  fields  C0<  and  Cq>  will  be  reconstructed  in  the  same  half 
space  with  respect  to  the  hologram  P,  even  though  they  were  originally 
in  opposite  ha]f  spaces. 
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Figure  IV.  3.  A  General  Boundary  Configuration  for  Holographic 

Imaging,  Reconstruction  from  Transmission  Holograms 
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Thus,  a  peculiar  feature  of  holographic  imaging  becomes  evident 
and  the  necessity  for  partitioning  the  scattered  object  field  £q 
according  to  Equation  (TV. 8)  is  clear.  Since  Cq<  and  Kq>  originally 
propagated  in  opposite  half  spaces  with  respect  to  the  hologram  but  are 
reconstructed  in  the  same  half  space,  obviously  the  holographic  process 
operates  on  interior  fields  differently  from  exterior  fields.  For 
example,  referring  to  Figure  IV. 5  ,  consider  that  £q  is  everywhere 
diverging;  i.e.,  Cq<  and  Kq>  are  both  diverging,  hence  represent 
virtual  images  of  object  0.  Exact  reproduction  of  Cq<  everywhere 
requires  that  over  some  surface  O',  geometrically  similar  to  object  0 
(can  be  scaled  in  size  or  rotated  in  orientation) ,  an  exact  duplicate 
of  the  object  boundary  field  £0<(£q)  be  formed.  Since  propagation  of 
the  reconstructed  field  is  clearly  from  left  to  right  this  formation 
represents  a  real  rather  than  a  virtual  image.  Similar  reasoning 
applied  to  SQ>  shows  that  it  reconstructs  as  a  virtual  image. 

There  are  other  unusual  effects  associated  with  reconstruction 

40 

of  the  interior  field.  The  most  obvious  is  pseudcscopic  inversion 
of  the  real  image  when  it  is  viewed  directly  rather  than  projected  on 
a  screen.  However,  this  is  not  a  fundamental  distortion  of  the  image. 
More  serious  is  the  absence  of  very  high  spatial  frequency  details  in 
the  reconstructed  object  caused  by  evanescent  propagation  of  these 
details  from  object  0  to  hologram  during  recording  and  from  hologram  to 
image  0*  during  reconstruction.  High  spatial  frequency  details  are 
missing  from  the  exterior  fields  so  their  absence  in  the  reconstructed 
exterior  fields  is  not  significant.  Of  course,  many  other  situation*5 
may  occur  but  the  foregoing  example  is  representative  of  the 


100 


?WIPS 
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fundamental  differences  which  may  exist  between  holographic  imaging  of 
internal  and  external  fields,  and  serves  to  illustrate  that  the  normal 
mode  formalism  may  only  be  valid  for  reconstruction  of  the  exterior 
field.  Therefore,  anticipating  that,  at  the  very  least,  interior 
fields  will  require  different  mathematics  than  exterior  fields,  the  two 
are  separated  in  the  formalism. 

One  further  comment  deserves  consideration.  Clearly,  if  the 
above  statements  are  true  for  the  surface  S,  they  are  also  true  for 
the  coterminal  surface  S'  or  any  other  cotermir.  1  surface  lying  wholly 
within  the  material  P.  Thus,  a  volume  recording  medium  provides 
redundancy  in  the  recording  of  optical  information.  It  has  been 
demonstrated  that  this  redundancy  can  be  used  to  increase  image  bright¬ 
ness.  For  example,  thick  emulsion  (multilayer)  holograms  have  a 
theoretical  efficiency  of  100%  whereas  conventional  thin  emulsion 

* 

(diffraction  grating)  holograms  have  a  maximum  efficiency  of  only  25%. 

A  concomitant  use  for  the  redundance  is  to  relax  the  constraints  on  the 
reconstruction  illumination  E„,  in  effect  making  the  decoding  process 
for  thick  emulsion  holograms  even  simpler  than  Gabor’s  original 
method4  * . 

4.  GENERAL  FORM  OF  TIE  HOLOGRAM  0PERA10R  t 

In  Section  2  a  hologram  was  defined  as  a  generalized  optical 
record  which,  when  properly  interrogated,  reconstructs  a  visible 


*  The  classification  of  holograms  as  being  multilayers,  diffraction 
\  gratings,  or  Fresnel  zone  plates  is  due  to  Professor  Nicholas 

|  George  of  Caltech. 

| 

| 

£ 
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replica  of  the  original  optical  field,  hopefully  at  all  points  in  space 
(Equation  IV. 6).  Geometrical  scaling  and  reorientation  of  the  recon¬ 
structed  field  relative  to  the  original  field  is  acceptable  and  because 
the  final  result  is  to  be  a  visible  field,  the  quality  of  the  recon¬ 
struction  can  be  judged  only  in  terms  of  field  intensity  (Equation 
IV.  7). 

The  most  convenient  method  of  interrogating  any  optical  record 
is  by  reillumination  of  the  record  with  visible  radiation  ^c(r)e27,ivt. 
For  example,  the  information  contained  in  a  photograph  is  conveyed 
quite  simply  to  the  observer  by  visual  inspection  of  the  developed 
emulsion,  which  requires  only  that  the  photograph  be  well  lighted.  In 

this  situation  incoherent,  diffuse  sources  are  superior  to  all  others. 

7  8  9 

Gabor  ’  ’  recognized  that  reillumination  could  also  be  applied  to 
holograms,  in  a  process  he  called  reconstruction,  and  stated  conditions 
on  the  illumination  required  for  this  purpose.  In  contrast  to  con- 
ventioital  photography  the  illumination  must  be  coherent  and  the  phase 
contours  must  be  carefully  controlled.  However  the  reconstruction  is 
three-dimensional . 

If  we  demand  that  the  interrogation  process  be  one  or  reillumina¬ 
tion  then  in  view  cf  the  discussion  presented  in  Section  3  we  can  state 
that  any  optical  record  P(r)  of  the  holographic  type  must  be  such  that 
upon  illumination  with  ^ (rj  a  replica  of  the  original  boundary  field 
£g(Tg)  is  produced  on  the  surface  S  of  the  hologram.  It  seems  quite 
likely  that  if 

Pdg)  «  50(rs)  (IV.  10) 
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the  required  replication  will  be  accomplished.  More  generally,  the 
hologram  will  be  a  weighted  record,  the  simplest  form  being 

P(rs)  *  €0(rs)£p(rs)  (IV. 11) 

* 

where  is  the  weighting  factor.  If  the  hologram  is  thick  then 
(IV.  11)  will  hold  throughout  the  volume  of  the  developed  hologram;  if 
the  hologram  is  thin  (IV.  11)  holds  in  the  plane  of  the  hologram.  In 
the  latter  case  the  operator  will  be  of  the  form 

Vis)  =  f  (  )4ci)s3Cl-Is)d3r  (IV.  12) 

M 

where  the  empty  parenthesis  inside  the  integral  indicates  that  the 
operator  has  meaning  only  when  it  operates  on  some  field  £q  which  is 
placed  in  the  parenthesis,  and  63  is  the  three  dimension- 1  delta 
function. 

The  hologram  may  be  considered  thin  whenever  its  thickness  is 
imich  less  than  the  minimum  fringe  spacing  to  be  recorded  .  Although 
the  validity  of  this  approximation  mist  be  carefully  weighed  when 
dealing  with  visible  holograms^,  it  is  almost  certainly  true  for 
acoustic  holograms  for  two  reasons.  First,  acoustic  sensors  used  to 
generate  acoustic  holograms  are  usually  no  more  than  one  acoustic  wave¬ 
length  thick.  Second,  the  emulsion  of  an  acoustic  hologram  trans¬ 
parency  is  usually  much  thinner  than  the  minimum  fringe  spacing 
recorded  thereon.  Thus  (IV.  12)  applies  to  the  recording  of  acoustic 
holograms. 
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Exact  calculation  of  the  field  £j(r)  diffracted  by  a  developed 

thick  hologram  transparency  P(r)  illuminated  by  a  field  ^Cr)  is  a 

formidable  task  involving  the  scattering  of  electromagnetic  radiation 

from  three-dimensional,  semi- transparent  volumes.  Rigorous  solution 

requires  a  multiple  order  scattering  theory  including  Bragg  reflection 
44 

phenomena  .  However,  for  thin  planar  holograms,  particularly  those 

A 

with  binary  format  SjCtj)  may  be  written 


5 1  J  p  (Is) G  (ij 


(IV.  13) 


where  according  to  the  Rayleigh-Saanerfeld  theory 

.iklirisM 


G(M-Is)  *  "  Sn 


-  ik|r 
e 

IW 


I  J 


(IV.  14) 


_l  U_  ,  V  [drls^l 

271  lij'isl  \  III -Is! '  iij'isl  ~  J 


and  n  is  the  unit  vector  normal  to  S  pointing  into  the  half  space 
45 

containing  Tj  ;  or  according  to  the  plane  wave  expansion  theory  (see 
Appendix  C) 


CO 

G tel -is1  *  /e 


ik,ter2s)  «-(rrIs), 


d  k 


(IV.  15) 


and  k_  and  <  are  the  longitudinal  and  transverse  wave  vectors  respec 
tively. 


Binary  formats  are  either  completely  transmitting  or  completely 
opaque;  there  is  no  gray  scale.  Acoustic  holograms  tend  to  be  of 
this  type  although  such  holograms  may  not  provide  optimum  image 
fidelity. 
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The  hologram  reconstruction  operator  may  then  be  defined  as 

=  /  (  )Cc(rs)G(rI,rs)d2rs  (IV.  16) 

S 

and  the  hologram  operator  H  =  3jH7  becomes 

“111)  •  f  f(  CIV.  17) 

M  S 

so  that  perfect  wavefront  reconstruction  may  be  written  as 

=  =  h_/* J 

M  S 


(IV.  18) 


dVrS 


It  is  important  to  determine  over  what  region  (IV.  18)  is  valid. 
Consider  a  single  point  source  test  object  located  at  Tq.  Then 
ikIrs-roi 


(IV.19J 


If  a  hologram  of  this  test  object  is  illuminated  by  a  field  with 

e  lvt  time  dependence  then  in  the  hologram  plane  the  field  will  be 

ik!r0-rr 

r 

which  is  in  the  same  form  as  a 


proportional  to 


e'l5te 


-S 

27T!iS'io* 


spherical  wave  diverging  from  the  hologram  plane  into  the  exterior 
region  Thus,  (IV. 18)  can  only  describe  reconstruction  of  the 
exterior  field  Cq>.  Since  by  definition 
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(IV.  20) 


~  ^Ck  ' 


Equation  (IV.  18)  is  more  properly  written 


^o>(lp  "  h  J  f  ^^rO^C^-S^G^-I*-S^5  rod~ 


M  S 


(IV.  21) 


If  the  factor  is  replaced  with  its  complex  conjugate  5Q€R 
then  upon  reconstruction  the  field  in  the  hologram  plane  will  be 

-iJt  iklls-^l 

P  A 

proportional  to  - 2iT{r  -r  j —  which  is  in  the  same  form  as  a 

spherical  wave  converging  from  the  hologram  plane  to  a  point  in  the 


exterior  region.  In  this  case  (IV. 18)  describes  reconstruction  of  the 
interior  field  £q<  except  that  whereas  the  former  converges  to  the 
point,  the  latter  diverges  from  the  point.  Equation  (IV. 18)  should 
then  be  written 

=  h  J J  ^0< -*  ^C^S^-I  »Is^  ^-0~-S^ d  rOd  rS  * 


M  S 


:iv.22) 


For  any  particular  tidjisverse  plane  r^  Equations  (IV.  21)  and 
(IV. 22)  are  of  the  form 


?(lj)  =  X  j <i(rs)K(rI,_rs)d?i 


(IV.  23) 


Now  in  general  5(rj)  f  C(lg)  so  that  (IV. 23)  is  a  Fredholm  equation  of 
the  first  kind  with  kernel  K(rj,rg)^d.  In  general  K  will  depend  on  z 


I 
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and  Zg  so  (IV. 24)  represents  a  family  of  equations,  one  for  each  image 
plane.  Thus,  (TV. 23)  does  not  appear  to  be  a  tractable  formulation  of 
holography. 

Moreover  (IV.  23)  is  incomplete .  Unique  reconstruction  of  CqQt) 
in  14  requires  not  only  knowledge  of  the  boundary  field  CqC^)  but  the 
equation  of  motion;  i.e. ,  k,  as  well,  This  is  particularly  important 
in  acoustic  holography  wherein  the  k  for  recording  is  much  smaller 
than  the  k  for  reconsti action;  only  the  k  for  recons truct ion  is  con¬ 
tained  in  (IV. 23). 

Nevertheless,  in  an  important  number  of  applications  the  two 
wave  numbers  are  the  same  and  it  is  then  apparent  from  (IV. 23)  that  we 
must  have 

4(rs)^(rs)  =  const  .  (IV.  24) 

It  is  important  to  investigate  the  significance  of  the  relation 
(IV.  24)  between  reference  and  reconstruction  illumination.  Expressing 
the  fields  in  terns  of  modulus  |£j  and  arguement  i  (IV. 24)  becomes 

URirs)je  K^|Cc(rs)|e  =  const  .  (IV.25) 

Clearly,  if 

|SR(rg)j  =  const  (IV. 26) 

then  by  choosing  -  £R  condition  (IV. 24)  is  satisfied  identically. 

The  most  obvious  of  the  forms  (TV. 26)  is  a  single  plane  wave  propagating 
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in  any  direction.  Another  is  uniform  diffuse  illimination  character¬ 


ised  by  random  variation  of  the  phase;  again  we  mist  take 

The  foregoing  are  special  but  important  situations  frequently 
encountered  in  holography,  he  now  consider  core  general  reference  and 
reconstruction  illumination  including  the  so-called  "ghost  imaging"  of 

47 

Collier,  et  al  ' .  By  writing  the  reference  aid  reconstruction  in 
terms  of  general  plane  wave  expansions  one  obtains 

i(iCTC/)  *£g 


5rCIs)5c(Is>  */ f 


(IY.27) 


*  e 


i[k,r<;-k2(<*)]2  , 

d  kd*k  * 


where  hRUj  and  Wc(>c)  are  the  complex  amplitudes  of  the  cth  directed 
plane  wave  components  of  £R  aid  £c  respectively,  without  loss  of 
generality  we  say  take  =  0.  Changing  variables  according  to  the 
prescription 


*L*  "  <  +  8 

?  ? 
d”jf'  =  d“S 


(rv.2S) 


then  gives 


CD  CD 

iR(lsKc(ls)  =  j  j  h'c(<)hp(^-5)d2, 


lg*£e  7 

1 


(IV.  29) 
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Taking  $  =  2rs  vs  see  trst  (IV.  29}  is  in  tie  loss  of  an  isrerse 

* 

Spatial  Fourier  Transfers  so  is  related  to  &  Spatial  Fourier 

Transfers  of  ^  4  Comities  (TV.  24)  thsi  beeesas 


J  *ct*)xR*t<*Bd2<  =  0 

0 


13!  >  £ 


lSf  <  s 


(•¥.3} 


vfjsre  e  «  1/a  and  2a  is  the  holcgraa  sire.  If  ^  and  5^  are  sssne- 
lated  spectral  distributions;  i.e. ,  ^  and  5-  are  safficissthr  dirfose, 

l.  K 

then  (IV. 53)  will  vanish  for  all  3.  In  order  that  the  correlation  rot 
vanish  everywhere  vs  ant  take  a^.  =  Thus,  if  a  reference 

is  used  the  identical  reference  isst  be  used  for  tie  iHurinsticn  sir¬ 
ing  reconstruction. 

Jiate  that  in  general  a  discrete  sue  of  plans  vanes  as  reference 
cannot  be  race  to  satisfy  (IV.30).  The  cross  sodet  tans  sill  not 
vanish  and  there  will  be  a  background  noise  in  the  reconstruction., 

Xben  he  reference  is  spherical  there  vill  be  ^ieraticcs  of  the 

-  4S,49,53 

sags  . 

5-  CSSL4L  HS-S  GF  IS  HXCCSSI  SSaSIR,  II 

It  is  new  clear  that  any  forruiacicn  of  holographr  alceg  da 
lines  of  Bjatiss  (W.6)  ant  be  in  terns  of  an  integral  eq^tica. 
kfcst  of  the  irssci  theca  y  of  integral  equations  pertains  to  Fxessla 
sontios  of 

tens  of  equations  of  this  type  it  is  necessary  to  include  the  preesga- 
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sgeena  gia  .  In  order  to  represent  fjDlcgraprry  m 


tlcsi  of  the  field  frees  coject  plane  r^ 
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recording.  As  mentioned  in  the  preceding  section  by  so  doing  we 
uniquely  specify  the  reconstruction  problem. 

To  avoid  confusion  the  following  convention  will  be  adopted 
subsequently.  Acoustic  fields  will  be  denoted  A,  electronic  signals 
will  be  denoted  E  and  visible  fields  (light)  will  be  denoted  V.  Also 
we  shall  adopt  the  plane  wave  expansion  formalism  (IV. 15)  for  a  variety 
of  reasons.  For  example,  the  analysis  of  acoustic  holography  must 
depart  from  the  usual  small  angle  theory  of  visible  optics  for  two 
reasons.  First,  there  are  likely  to  be  "near  field"  or  'Vide  angle" 
situations  for  which  the  paraxial  approximations  oi  conventional 
Fresnel  and  Fraunhoffer  diffraction  theory  are  invalid,  the  experi¬ 
mental  situations  covered  in  Chapter  V  probably  falling  in  this  cate¬ 
gory.  However,  conventional  paraxial  theories  will  remain  useful  for 
calculating  focal  lengths  since  the  concept  of  focus  is  meaningful  only 
in  first  order  theory. 

Second,  in  visible  optics  no  account  is  taken  of  the  directional 
sensitivity  of  the  molecules  making  up  the  grains  of  the  photosensitive 
emulsion  since  this  effect  is  integrated  over  many  molecules  and  grains 
randomly  oriented.  However,  for  the  bulk  of  acoustic  work  (including 
the  experiments  presented  in  Chapter  V)  a  different  situation  exists 
because  acoustic  fields  are  detected  by  piezoelectric  transducers  A 
correct  analysis  of  this  case  must  include  the  directional  character¬ 
istics  of  the  transducer  and  the  sampling  properties  of  the  scanning 
procedure.  The  directional  nature  of  the  receiver  is  determined  by 
measuring  its  response  to  an  acoustic  plane  wave  as  a  function  of  the 
angle  of  incidence.  Obviously,  any  factor  which  depends  upon  wa.e 


number  is  most  easily  incorporated  if  the  plane  wave  expansion 
formalism  is  employed. 


Finally,  the  most  elementary  holograms  (holographic  diffraction 
gratings)  are  those  made  by  the  interference  of  two  plane  waves.  These 
elementary  holograms  are  the  basic  building  bloc?.s  from  which  more 
complex  holograms  are  constructed.  Experimental  study  of  the  proper¬ 
ties  of  such  holograms  provides  the  basic  data  for  holographic 
research5*.  Thus,  it  seems  natural  to  formulate  the  theory  of  holo¬ 
graphy  in  a  way  which  emphasizes  this  fact  and  which  is  related  to  the 
most  basic  holographic  measurements. 

The  acoustic  field  scattered  into  the  hologram  plane  £g  from 
the  object  plane  r^  is 


A°tis)  =  ui)2  ff*° 


A  0  * 

(IV.  31) 


In  the  S  plane  a  hologram  is  formed  by  causing  A^Crg)  to  interfere  ei¬ 
ther  acoustically  or  electronically  (see  Section  V.4)  to  form  an  inter¬ 
ference  pattern  which  is  imperfectly  converted  into  an  electror  .c  sig¬ 
nal  by  a  piezoelectric  detector  with  plane  wave  response  D(k^) •  The 
desired  interference  signal  looks  1ike 


0ClX(ls>  ’  Ke  ^  Jf  WDC^in  &r) 


(2it) 

m  m . 

u  / 


f.lS-To>  iklA(zS'zO>  2,  .2 
e  e  d  k.q  r. 


(IV.  32) 
2 

A“  r0 
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where  III(x)  denotes  the  infinite  array  of  evenly  spaced  one  dimen¬ 
sional  delta  functions 

CO 

III(X)  =  ^  <$ (x-n)  ,  (IV. 33) 

—CO 

and  Kg  is  a  conversion  constant. 

To  produce  a  hologTam,  the  field  over  the  entire  hologram 
aperture  must  be  detected  and  recorded,  either  by  covering  the  aperture 
with  sensors  or  else  by  mechanically  scanning  one  or  several  over  the 
aperture.  An  array  involves  two  dimensional  sampling,  linear  scanning 
only  one  simensional  sampling.  The  delta  function  arrays  account  for 
the  sampling  procedure.  Vertical  scanning  (y  directed  lines  evenly 
spaced  d  apart)  was  adopted  for  the  experiments  presented  in  Chapter  V. 

The  interference  signal  (IV. 32)  is  next  converted  imperfectly 
into  a  visible  transparency  (the  acoustic  hologram)  according  to 

P(tg)  =  %  J S(rg,ryE0(mr^)Ej(nn^)d2r*  (IV.34) 

§ 

where  S(rg,rg)  represents  the  point  spread  function  (spot  size)  for 
the  optical  display  and  Kp  is  a  conversion  constant.  For  an  ideal  sys¬ 
tem  S(rg,rg)  =  $(rg-rg)  but  this  is  never  realized  in  practice.  For 
example,  in  the  experiments  presented  in  Chapter  V  the  acoustic  inter¬ 
ference  pattern  in  the  hologram  plane  was  displayed  on  a  cathode  ray 
tube  (CRT)  which  had  a  minimum  spot  size  of  0.2  mm  diameter. 

The  factor  m  accounts  for  a  linear  reduction  of  transverse 
acoustic  field  dimensions  by  m  to  those  of  the  hologram  transparency. 

A  linear  reduction  is  required  so  that  the  resulting  hologram  will  focus 
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inside  the  coherence  volume  of  the  light  used  for  reconstruction  and 
within  a  finite  sized  laboratory. 

The  acoustic  hologram  (IV.  34)  is  then  illuminated  with  coherent 
light  in  order  to  obtain  a  sc  newhat  imperfect  visible  replica  (recon¬ 
struction)  of  the  acoustic  field  A^.  The  reconstructed  replica  field 
in  plane  r.j  is  given  by 


V, 


co 

(ll)  ■  //pfes)e 


iv(^r^_lkzv(zrzs)  2  ,2 


d  Kyd  ig  . 

(IV. 35) 


Substituting  (IV. 34)  and  (IV.  32)  yields 


W  “5 *  // fff  VW^vcfr-S>D^ 


/mr'-x\  /mrl-y  \ 


(IV, 36) 


ikzA(W  W5 S}  ik2V(zl'zS) 
e  e  e  e 


,2  ,2  ,2  ,  ,2  ,2 

d  r0d  KAd  rSd  V*  rS 


Defining 
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Vii-V  Tjjt  ff'ffw PP) "' (¥) 


ijs/tirle}  ik7v(zT'zP  2  2  2  2 

^  1  Se  zV  1  S  dK^r’dy  rg 


(IV.  37) 


Equation  (IV. 36)  may  be  rewritten  as 


w 

VI'-P  =  /  V^Wd  r0 


(IV.  38) 


Now  holography  will  only  be  useful  if  it  is  possible  to  select 
pairs  of  object  planes  such  that  in  these  planes  the  positional 
dependence  of  Vj  and  Aq  is  geometrically  similar.  That  is  we  desire 


hVj(r)  =  A^r) 


(IV. 39) 


Tnen  (IV.  38)  becomes 


CO 

Vii>  ” h  / 


**1fci>*b,V*<Pd‘ro 


(TV. 40) 


Now  (IV.  58)  can  be  reduced  to  (IV. 40)  only  by  selecting  special 


combinations  of  zT  and  zn;  i.e., 
i  u 


ZI  = 


(IV.  41) 
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where  f  represents  a  functional  dependence  yet  to  be  determined.  The 
relation  (IV. 41)  is  known  as  a  focusing  condition.  When  (IV. 41)  is 
satisfied  reduces  to  a  form  satisfying  (IV. 40).  However,  even  when 
the  focusing  condition  is  satisfied  (IV. 40)  may  not  be  valid  for  any 
field  configuration  Aq,  only  certain  special  ones  A^.  Hence  (IV. 40) 
is  more  properly  written 

co 

-  h„  /  Vii’Jo>Wi(PdZro  f"-42) 

-CO 

where  h^  represents  the  relative  "strength"  of  the  corresponding 
configuration  A^  in  the  reconstruction. 

Equation  (IV. 42)  is  immediately  recognized  as  a  homogeneous 
Fredholm  equation  of  the  second  kind  with  kernel  K (r  j ,  i^) ,  eigen¬ 
functions  Aqh ( j_qO  and  eigenvalues  hn. 

Over  a  period  of  years  numerous  authors  have  contributed  to  the 
theory  of  Fredholm  equations.  Certain  of  the  properties  of  Fredholm 
equations  of  special  significance  to  holography  are  stated  below  with- 
out  proof-'  . 

1.  Corresponding  to  the  homogeneous  Fredholm  equation 

£(r)  =  X  J  K(r,r'K(r’)d2r’  (IV.  45) 

I 

construct  the  inhomogeneous  Fredholm  equation 

C(r)  =  g(r)  -  J K(r,r’K(r')d2r'  (IV. 44) 

I 
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where  g(r)  is  arbitrary.  Then  either  the  inhomogeneous  Equation 
(IV.  44)  has  a  unique  solution  for  any  g(r),  or  the  horr'veneous 
Equation  (IV. 43)  has  at  least  one  nOn- trivial  solution  £^(r)  with 
eigenvalue  X. 

2.  Any  function  A(r)  which  can  be  represented  as  a  "source” 
distribution  in  which  K(r,r_')  nlays  the  role  of  a  Green's  function; 
i.e. , 

A(r)  =  j  K(r,r*)B(r')d2r*  (IV.4S) 

where  B(r')  is  arbitrary,  can  be  expanded  in  terms  of  the  eigen¬ 
functions  £n(r)  of  K(r,r* ) . 

Although  it  is  not  known  if  solutions  to  (IV. 42)  (perfectly 
reconstructed  field  configurations)  exist  in  the  general  case,  property 
1  suggests  that  such  solutions  should  exist.  The  second  property  of 
Fredholm  equations  indicates  the  possible  use  of  such  configurations 
as  expansion  functions. 

6.  APPROXIMATE  SOLUTIONS  TO  THE  EIGENVALUE  EQUATION 

In  order  to  reduce  the  eigenvalue  equation  to  a  form  which  is 
readily  solved  we  naist  determine  the  focusing  condition  (IV. 41) .  This 
is  most  easily  accomplished  when  the  fields  are  paraxial.  Then  the  k_ 
ray  be  expanded  in  a  power  series  and  only  the  lower  order  terms 
retained.  Although  the  recording  of  an  acoustic  hologram  is  not 
generally  a  paraxial  problem  a  simple  arguement  will  show  that  the 
reconstruction  from  an  acoustic  hologram  is  always  a  paraxial  situation. 

Consider  an  acoustic  plane  wave  travelling  in  the  8^  direction 
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(wave  number  k^,  transverse  wave  number  where  sin  0^  -  xr^/k^) 
generated  from  a  grating  of  spatial  frequency  k^/2tt.  If  the  grating  is 
reduced  in  size  by  m  and  then  illuminated  with  light  at  normal 
incidence  (wave  number  ky  =  yk^)  then  the  grating  will  diffract  a 
plane  wave  in  a  new  direction  6V  given  by 


<y 

Sin  6V  =  k^  = 


m  .  . 

Uk7  =  u  Sin6A 


(IV.  46) 


Typical  numbers  for  acoustic  holograms  (at  1  J-f-iz  in  water)  are  m  =  30, 
y  -  3*10^,  m/y  =  10  Thus,  even  though  recording  an  acoustic 
hologram  is  likely  to  involve  large  angles  8^,  the  corresponding 
reconstructed  angles  will  always  be  small,  hence  paraxial. 

In  the  paraxial  approximation 


x2)1'2  * 


k 


1 

'  ir 


(IV.  47) 


For  simplicity  and  without  loss  of  generality  take  the  origin  of 
coordinates  to  lie  in  the  hologram  plane  so  that  z^  =  0.  Then  the 
integration  in  Equation  (IV. 37)  becomes 


(ky/2~Zp)  jmrg 


*  ?:=£.s 


(IV.  48) 


where  we  have  used  the  relation 


(IV.  49) 


In  the  more  general  case  the  paraxial  assumptions  are  invalid 
for  describing  the  recording  of  an  acoustic  hologram  and  instead  of 
(IV.  48)  ve  imst  write 


K. 


D(^) 


■(mi’ 


zA20 


.2 

a  k.  = 
A 


G(mrg  -  r^)*p[(rar4  -  r^/2-] 


(IV.  50) 


where 


G(r)  =  e 


ik^zo 


(IV.  51) 


Applying  the  paraxial  procedure  to  the  Cy  integration,  assuming 
that  the  exponential  factors  in  vary  slowly  in  r£  conpared  to  the 
display  spread  function  S(r<.,rg),  and  assuming  display  spot  site  raid 
shape  is  independent  of  position  (not  true  for  electrostatic  CRT 
displays  unless  restricted  to  the  center  portion)  finally  yields 


i-  fr  _  ,  _  VVVV2^  *  kAz(P  i"(VfezI)rI 


*0*1 


CD 

•  J  n(xs/X)I!(ysA')i?<j(rs)ER(5^s)Vc(rs)*  (I\r-  52) 


i-tKj/2-7.^ -STg  -  rQs 


*p; (mr,  - 


CV-I^  2-ils'C^/^i)  .2 

e 


a  r~ 
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where  the  finite  size  2X,  2Y  of  the  hologram  aperture  is  included  in 
the  function  n(Xg/X)II(yg/Y) ,  and  where 

/?dCrs)  =  S(rs)*[III(xs/d)III(ys/d)3  (IV.  S3) 


is  the  raster  function  representing  the  composite  effect  of  sampling 
and  finite  display  spot  size. 

To  evaluate  (IV. 52)  we  must  choose  and  V^.  Usually,  sources 
of  acoustic  radiation  tend  to  be  point  rather  than  planar  sources. 
Moreover,  it  will  prove  advantageous  to  use  spherically  converging 
wavefronts  during  reconstruction.  Accordingly,  we  consider  the 


special  but  important  case  wherein 
Er^)  =  Er  e  e 

i(k../2zr)  |r.  -  rj ' 

VC{TS)  =\?ceVLe^C:>~C  ; 


(IV.  53) 
(IV.  54) 


\r  is  a  spherical  wave  with  focus  located  at  E^  appears  to  be  a 

sphericao  wave  with  focus  located  at  r^.Zj,.  E^  say  originate  frcn  an 
actual  acoustic  source  or  it  cay  be  an  electronic  simulation  of  an 
acoustic  source.  (Electronic  simulation  is  usually  confined  to  planar 
references.)  Kith  respect  to  the  hologram  plane.  Focal  points  z  <  0 
represent  diverging  waves ,  focal  points  z  >  0  represent  converging 
waves. 

Substituting  (IY.55)  and  (IV. S4)  into  (IV. 52)  yields  upon 
commuting  the  convolution  with  the  r^  integration 
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„  ,  %  >AVkCERVC^  *M2I  +  Z(P  'ikA^Z0  +  Tp 

W^r  *>>  " - qjq - e 

i(k^./2zc)r2  -i(kv/2zI)rj  iCk^Zz^r^  -iG^/ZZjjJrJ 


f  i(M2-o)v2  4 

/  D(v)e  A  ^  e 


(W^cp  d2v 


(IV.  SS) 


■  /  n(xs/X)n!ys/Y)~d(rs)e 


t(^v£l/^2j)  '  (k^^2=*c)I 


2sirs*[{sk»r^/2^2R)  -  [y^^lzz^  -  (k .^v/ZsSq)] 


i-f(^r/2rzc)  -  (|^/2s2j)  ♦  (e^/^z^  - 


where  v  is  a  dunsy  variable  required  for  the  convolution  operaticn. 

The  quadratic  phase  tent  appearing  in  (IV.  55)  say  be  eliminated 
by  taking  the  image  plane  z  -  Zj  to  be  located  at 


1  s‘0  -  z€/zR)  -  uz^/z. 


(Tv.  S6) 


Equation  (IY.S6)  is  the  form  wuch  the  focusing  condition  (IY_4I)  takes 
for  the  approximate  calculation  xe  are  considering  here.  This 
condition  is  quite  useful  because  the  concept  of  focus  is  really  only  a 
first  cider  rhenctsena  anyway  - 
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Ft- 


Y 


Kith  elisisaticn  of  the  quadratic  phase  factor  is  the  rc 
Integration  £  _  beetles 


yipIcK  =  amt  e 


i(k*/2jbSro 


•  f  c(v; 


BMel(kA'  JVv*  e?71L  .  (IV.S5) 


■{siKlZlXtvJ  -  -  VM1- 


•{Sinc[2rV(v-  -  ^V'tajJj.AjCT;  -  k^V/21^3  d*v 


vhere 


v*  =  \£i/2tci  - 


sine  x 


sin  x 
x 


const  = 


4SE£PkiWcX!(2T)i  ^U!  *  V)  -*.(=0  * 


C_1 

S' 


*tri 


i(V.y2rr)r£  -iCw2cf)r‘  H^/2c_)rJ 
:  e  e 


As  .v  and  Y  beccse  large  the  sire  nsscciess 

peaked  around 


(IY.5S) 


increasingly 
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WZI  '  "  ‘^A^c/2©  +  ‘^A^R^R  ’  kA-^Z0  "  °* 


(IV.  59) 


Further,  if  the  angular  response  D(^  of  the  detector  elements  is 

sufficiently  broad,  D(v)  vail  be  nonzero  only  in  a  narrow  range  about 
* 

v  =  0.  Thus  for  sufficiently  large  holograms  and  for  sufficiently 
broad  detector  response  behaves  approximately  like  a  two- 
dimensional  delta  function  and  the  Tq  integration  (IV. 36)  can  be 
carried  out  to  give 


VjCrj)  «  hA^r^i  +  o) 


i(V2zo)  tetf*  +  -I2 

e 


(IV.  60) 


where 


M  = 


m(l  -  Zq/zj^)  +  (y/mKzQ/z^  * 


^(p/m)  (zn/zr)  -  Tf.Zn/Zj} 


CIV.  61) 


t  4KEKpVt4VC30rC2,,)2  ikV(zI  *  ZC5  -ikAtz0  + 
n  =  — - “ 


Z0ZI 


i(k^/2zc)r2  -idCy/Zzpr2  -i(kA/2zR)rR 
*e  e  *  e 


Thus,  from  Equations  (TV. 56),  (IV. 60) ,  and  (IV. 61)  we  see  that  to 
first  order,  and  provided  that  the  assumptions  stated  above  are  valid, 

2 

*  In  order  that  the  ideal  detector  response  D(v)  -  6  (v)  be 
obtained  the  detector  would  have  to  respond  equally  to 
propagating  and  to  evanescent  waves. 


the  reconstructed  image  is  a  visible  replica  of  the  object  boundary 

acoustic  field  but  magnified  laterally  by  M,  displaced  laterally  by  o, 

■? 

and  magnified  longitudinally  by  uM“.  Similar  expressions  are  obtained 

*  *  _ 
if  AqE^  is  substituted  for  AqEr.  The  ‘.relations  are 


2  .  ,  pzo _ 

z  ~2  , 

m  (z(/\  "  ^  +  v\ /ZC 

»  _  _ _ -1 

m(l  -  Zq/z^J  rlu7m5Tz^7z^T  * 

rn(zn/zp)  -  rr(u/m)(zn/zr) 

a  •  m  - V4>  - 


Relations  (IV. 56),  (IV.  61),  and  (IV. 62)  have  been  obtained  through 
different  techniques  by  Meier^4’^’^.  To  put  these  equations  in  his 
form  substitute  1/m  for  m,  1/y  for  u. 

We  now  consider  a  special  case  which  is  of  particular  importance 
in  acoustic  holography  because  of  the  ability  to  simulate  planar 
references  electronically.  If 


ZR  *  “  »  ZC  ■* 


(IV.  63) 


such  that 


rR  rr 

— -  -*■  const  =  tan  6D  ,  —  const  *  tan  6r  , 
Zr  k  zC  l 


(IV. 64) 


where  9r  and  6^  are  the  angles  which  the  wave  vectors  of  the  reference 
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ana  reconstruction  plane  /  uvc-s  make  with  the  normal  to  the  hologram 
plane,  then  the  preceeding  relations  reduce  to 


zT  = 

A 


$ 


M  = 


1 


m 


9 


(IV.  65) 


(rR  tan  eR  -  r^i/m  tan  8c)zQ 


The  +  sign  indicates  a  real  image,  the  -  sign  indicates  a  virtual 
image.  The  last  of  relations  (IV. 65)  is  valid  only  for  paraxial 
reference  and  reconstruction  plane  waves.  In  the  more  general  case 
substitute  sin  for  tan. 

Fran  (IV, 65)  it  is  clear  that  in  this  important  case  lateral 
magnification  does  not  depend  upon  the  wavelength  ratio,  and  we  see 
the  desirability  of  minimizing  the  reduction  m  of  the  hologram 
between  recording  and  reconstruction.  We  also  note  that  there  is 
distortion  of  the  reconstructed  object  scene  because  the  longitudinal 
magnification  is  considerably  different  from  the  lateral  magnification. 
For  example,  spherical  objects  will  appear  as  oblate  objects  in 
reconstruction.  Depending  upon  the  parameters  this  effect  can  be 
extreme, 

7.  EFFECT  OF  DETECTOR  DIRECTIVITY 

The  directional  characteristics  of  acoustic  detecting  elements 
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restrict 


1.  xesoiution 

2.  field  of  view. 

condition  (IV.  59)  the  effect  of  the  fora  of  D(<^)  upon  resolution 
xs  clearly  evident.  As  the  angular  spectrum  D(<^)  becomes  restricted, 
D(v)  becomes  less  sharply  peaked  acO”t  v  r  0.  This,  in  turn,  requires 
that  v  assume  various  non-2ero  values  in  (iv.59)  thereby  leading  to  a 
reconstruction  which  is  a  superposition  of  various  ina^.,  each  at  a 
slightly  different  lateral  position.  Thus,  the  reconstruction  appears 
to  be  blurred. 

In  extreme  cases  this  effect  can  lead  to  the  appearance  of 
double  images.  For  example,  in  the  one -dimensional  case  if  D(v)  were 
a  double-humped  pattern,  idealized  as  say 


D(v)  =  II(v/a)  +  It[(v  -  b)/a] 


(IV.  66) 


where  a  is  the  width  and  b  the  separation  of  narrow  rectangular  peaks, 
then  a  double  image  would  be  formed.  The  corresponding  angular 
spectrum  for  the  detector  is  of  the  fora 


(IV.  67) 


a  multi-phased  response.  Thus,  one  must  give  due  regard  to  the 
phase  characteristics  of  the  detector.  In  the  case  of  piezoelectric 
detectors  this  means  desinging  the  structure  to  minimize  internal 
reflections. 

The  dependence  of  field  of  view  upon  detector  pattern  is  obvious 
although  it  cannot  be  obtained  from  the  paraxial  equations.  Clearly, 


objects  lying  at  angles  outside  the  pattern  of  the  detection  elements 
vill  not  be  imaged.  This  effect  is  probably  sore  pronounced  in 
holography  than  in  conventional  imaging  because  in  the  latter  each 
image  point  is  formed  by  rays  converging  to  the  point  from  all  angles 
inside  the  solid  angle  subtended  by  the  lens  aperture  ct  the  point 
(provided  the  incasing  radiation  is  diffuse).  Thus,  if  a  lew  f  issber 
lens  is  used  surely  some  portion  of  the  acoustic  radiation  at  each 
isage  point  will  lie  within  the  detector  angular  response  pattern.  It 
wo-Xd  be  interesting  to  consider  to  what  extent  holographic  techniques 
can  be  used  to  compensate  for  the  detector  angular  response. 

8.  S?bjLH-MKDSIDrri  PRODUCT  FOR  ACOUSTIC  HOLGGRA.MS 

Space  bandwidth  products  have  been  calculated  for  visible 
holography54 .  The  technique  is  to  assume  that  the  hologram  size  2X,  2Y 
is  always  so  Orge  that  the  number  of  Fresnel  zones  recorded  is  limited 
by  the  film  resolution  rather  than  the  film  size.  Thus,  the  limit  to 
resolution  in  visible  holography  is  essentially  the  film,  not  the 
aperture  size.  In  this  sense  holography  provides  very  great 
resolution  simply  because  of  the  extremely  fine-grained  emulsions 
which  are  available.  For  film  limited  holograms  DeVelis  and  Reynolds 
have  calculated  that  the  one-dimensional  space  bandwidth  product  for 
coplanar  reference  holograms  should  be  approximately 

SBP  =  1.64s  X  (IV. 68) 

m 

where  is  the  raaximm  spatial  frequency  which  can  be  recorded  by  the 
film  and  X  is  the  film  size. 

In  acoustic  holography  a  different  situation  will  often  be 
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encountered.  Here  iv  is  quite  likely  that  the  resolution  of  the  "filrr*1 
or  sampling  procedure  will  be  sore  than  adequate  and  that  the  limit  to 


resolution  will  be  simply  the  finite  issuer  of  Fresnel  zones  which  can 
be  acccsodated  in  a  relatively  sail  aperture.  For  aperture  limited 
holograms  the  Rayleigh  criterion  gives  an  approximate  value  for  the 
resolution.  Thus,  for  acoustic  holograms 


SB?  ^ 


X‘ 

1.22Az, 


j 


(IV.  69} 


vhere  Zq  is  the  distance  between  hologram  and  object  planes.  Equation 
(IV.  69}  is  also  valid  for  acoustic  holograms  in  which  the  reference  is 
introduced  externally  and  the  interference  occurs  electronically  in  the 
display  rather  than  acoustically  in  the  propagating  seditea.  Of  course, 
if  the  resolution  of  the  display  device  is  a  limiting  factor  then 
(IV. 68}  applies  but  with  s^  replaced  with  s^/m  where  s^  is  the 
resolution  limit  of  the  display  and  m  is  the  demagnification  ratio. 

Certain  assumptions  relative  to  establishing  the  contrast  level 
which  constitutes  a  "resolved"  fringe  are  implicit  in  Equations  (IV. 68) 
and  (IV. 69).  Clearly,  spatial  coherence  and  film  or  display  MIF  are 
important  in  this  connection.  Here  it  las  been  assumed  that  the  fields 
exhibit  complete  spatial  coherence  and  that  the  MTF  is  uniform  up  to 
the  cutoff  sm  so  that  fringe  contrast  is  either  one (resolved)  or  zero 
(unresolved) . 

It  is  interesting  tc  determine  the  aperture  size  for  which  the 
resolution  becomes  detector  or  film  limited.  Equation  (IV. 68)  to 
(IV.  69)  we  obtain 
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X.  *  2s_Xa~ 

t  2  u 


Taking  s_  =  I/NX, 


?-• 

v  _^0 

xt  -  "N“ 


(iV.70) 


(IV.  71) 


which  is  independent  of  wave  length.  Typically  I  <  N  <  2  for  acoustic 
holograms  so  for  any  reasonable  range  the  aperture  required  for 
detector  limited  performance  will  be  enormous,  prohibitively  large. 

For  visible  holograms  1  <  N  <  1C  so  fils  United  performance  is  indeed 
possible. 

9.  SUSfSUlY 


We  have  presented  a  formulation  of  holographic  imaging  in  terms 
of  an  eigenvalue  equation  taking  the  fora  of  a  homogeneous  Fredhola 
equation  of  the  second  kind.  We  give  the  explicit  form  which  the 
equation  assumes  for  acoustic  holograms,  including  such  factors  as 
detector  directivity,  and  present  an  analysis  of  the  acoustic  equation 
under  paraxial  assumptions.  Even  in  this  case  there  is  longitudinal 
distortion  due  to  unequal  magnification  of  transverse  coordinates 
relative  to  longitudinal  coordinates.  If  the  paraxial  assumption  is 
invalid  then  the  problem  becomes  even  more  complicated  and  in  general 
all  linages  will  suffer  distortion  and  aberration.  It  remains  to  be 
seen  if  the  general  eigenvalue  equation  can  be  solved  in  the  sense  of 
finding  a  set  of  object  field  configurations,  the  eigenfunctions, 
which  are  perfectly  imaged;  or  even  to  prove  the  existence  of  such 
configurations . 


CHAPTER  V 


EXPERIMENTS  IN  ACOUSTIC  IMAGING  BY  H803USIS f 


iKiRoacriQi 


tie  have  presented  theories  of  scalar  i  raging  and  holography 
emphasiiing  acoustic  applications.  The  ultimate  goal  of  acoustic  imag¬ 
ing  is  to  obtain,  in  real  time,  visible  Irages  of  objects  imbedded  is 
media  opaque  to  all  but  acoustic  fields.  Ke  now  present  experissents 
demonstrating  attainment  of  visible  inages  of  underwater  objects  by 
acoustic  means  using  the  techniques  of  holography.  This  method  is 


known  as  acoustic  hole 


■irwiv 


Because  of  the  nature  of  the  experiments  and  difficulties  in 


performing  then,  it  has  not  been  possible  to  illustrate  each  aspect  of 
the  theory  with  an  appropriate  experiment.  Moreover,  this  research  vas 
not  intended  to  yield  an  izsediately  practical  Impiimentation  of  holo¬ 
graphy  into  acoustic  imaging.  Rather,  we  hope  for  the  core  limited 
goals  of  shewing  that  visible  images  of  underwater  objects  can  indeed 
be  obtained  by  acoustic  holography,  even  with  simple  equipment,  and  of 
studying  various  properties  of  acoustic  holograms. 

Before  proceeding  to  a  detailed  description  of  the  experiments , 
we  must  provide  a  motive  for  acoustic  holography.  Acoustic  imaging  is 
a  relatively  new  field,  the  first  significant  contributions,  due  to 
Sokolov,  coming  only  in  1936^.  Until  recently,  activity  in  acoustic 
imaging  has  been  restricted  to  an  occasional  effort  along  conventional 
lines;  i.e.,  using  acoustic  lenses  or  mirrors,  or  most  often  merely  an 
acoustic  shadow 58,59,50,61^  results  have  been  only  moderately 
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successful  22s,  is  any  case,  the  ranges  were  each  shorter  than  respired 
for  undersea  applications. 

Rapid  progress  in  acoustic  imaging  has  been  issodsd  by  t?io 
problscs:  (1)  lack  or  suitable  acoustic  detectors,  and  (2)  difficulty 
with  design  and  fabrics  tics  of  abb-statics  free  accustic  lenses  and 
sirrors  vita  lev  acoustic  loss  and  large  aperture.  Holography  provides 
a  way  around  the  seccsd  problem  sisnily  by  eliminating  the  need  for  a 
lens.  Ccrrventicnsi  imaging  la  obtained  by  intredneing  a  carefully  ccn- 
trolled  angular  dispersion  into  the  wave  field  by  causing  it  to  pass 
through  a  precisely  contoured  refractive  leas  material.  In  holography, 
the  dispersion  process,  which  presents  ?  serial  and  fabrication  preb 
less,  is  replaced  by  comparison  with  an  accurately  knesn  phase  refer¬ 
ence.  In  accustic  holography,  as  will  be  demonstrated,  this  reference 
can  easily  be  generated  electronically. 

Other  useful  properties  say  be  attributed  to  acoustic  holography. 
A  hologram  contains  information  about  all  cress  sections  of  the  object 
scene  at  right  angles  to  the  direction  of  wave  propagation,  even  if  the 
hologram  recording  surface  is  two  dimensional.  All  of  this  information 
is  retained  in  a  single  hologram.  Conversely,  a  lens  requires  refocus¬ 
ing  of  the  system  for  each  object  plane  outside  of  the  depth  of  field, 
and  only  one  depth  of  field  at  a  time  may  be  permanently  recorded. 

Echo  ranging  or  sonar  techniques  require  very  short  duration  acoustic 
pulses  in  order  to  accurately  determine  object  positions  along  the 
direction  of  propagation  whereas  holograms  may  be  made  using  continuous 
as  well  as  pulsed  illumination. 
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particularly  troublesome  because  there  is  not  likely  to  be  resolution 
to  spare  due  to  the  long  wavelengths  involved  with  acoustic  imaging. 

Real  time  acoustic  imaging  along  conventional  lines  using  the 
Sokolov  device  has  had  some  success  in  ultrasonic  flaw  detection  and 
other  short  range  applications*^,  but  this  device  does  not  appear  to 
have  sufficient  resolution  or  aperture  for  use  as  the  detecting  surface 
in  holographic  imaging.  Other  efforts  at  real  time  acoustic  holography 
have  been  made,  but,  for  the  most  part,  acoustic  holography  must  still 
be  considered  a  two  step  procedure  in  which  the  hologram  is  a  perma¬ 
nent  record  to  be  suitably  processed  for  reconstruction  at  a  later  time. 

The  great  need  in  acoustics,  which  has  yet  to  be  satisfied,  is  not 
for  sophistication  but  simply  for  any_ practical  scheme  which  enables  us 
to  "see"  acoustically.  Whether  acoustic  holography  will  ultimately 
prove  superior  to  conventional  or  other  techniques  for  this  purpose  will 
not  be  known  for  some  time. 

2.  PHYSICAL  ARRANGEMENT  OF  THE  EXPERIMENT 

The  experiment,  the  physical  arrangement  of  which  is  shown 
schematically  in  Figure  V.l  and  pictorially  in  Figures  V.2  through  V.5, 
is  performed  underwater  inside  an  anechoic  tank.  The  object  to  be 
visualized  is  '’illuminated"  by  a  5-cm-diam  barium  titanate  disk  driven 
at  1  MHz  (wavelength  X  =  1.5  mm  in  water)  (see  Appendix  D).  This  fre¬ 
quency  was  chosen  because  of  its  suitability  for  long  range  underwater 
acoustic  imaging.  It  appears,  upon  examining  Figure  III. 4,  that  1  MHz 
is  the  highest  frequency,  hence  smallest  hologram  aperture,  for  which 
the  attenuation  in  water  is  below  200  db/km.  Higher  frequencies,  where 


Figure  V.l.  Arrangement  for  Making  Underwater  Acoustic  Holograms 


C  Carriage 

H  Hologram  plane 

HC  Hydraulic  cylinder  (152  cm  stroke) 

L  Acoustic  illumination  source 

LB  Lathe  bed 

0  Object 

P  Probe  (457  cm  long  *  2.54  cm  dia) 

R  Receiver 

S  Tubular  supports 

T  Anechoic  tank  (3mx3mx7m  deep) 

WL  Water  level  (usually  5.5  m  deep) 
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Acoustic  Tank. Looking  Toward  Lathe  Bed  and  Hydraulic  Cylinder 


Figure  V.3.  View  of  Top  of  Acoustic  Tank 

Plane  of  hologram  coincides  with  opening  in  tank  cover. 
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Figure  V.S.  Detail  of  Scanning  Mechanism 

Position  sensing  potentiometers  are  at  upper  left  and 
lower  middle. 
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the  attenuation  is  much  greater,  wculd  be  more  suitable  for  biomedical 
and  flaw  detection  applications. 

The  scattered  acoustic  field  is  detected  by  a  3-nm-diam  bariwn 
titanate  disk  mounted  in  a  slender  probe  which  may  be  scanned  over  a 
130  cm  x  150  cm  vertical  plane  located  near  the  bottom  of  the  tank. 

(See  Appendix  D) .  The  hologram  plane  was  placed  at  the  bottom  jn  order 
to  minimize  possible  reflections  from  the  air-water  interface  at  the 
top  or  the  tank.  The  scanning  mechanism  consists  of  a  hydraulic 
cylinder,  controlled  by  a  linear  two-way  servovalve,  for  the  vertical 
motion,  mounted  on  a  lathe  carriage  and  bed  assembly  for  the  horizontal 
motion.  With  this  arrangement  any  hologram  size  and  spacing  of  the 
raster  lines  may  be  obtained.  To  minimize  probe  vibrations  and  water 
disturbance,  the  raster  lines  were  always  scanned  in  the  vertical 
direction  from  bottom  to  top. 

The  amplified  output  of  the  receiving  transducer  is  combined  with 
an  electric  reference  by  suitable  processing  (discussed  in  the  follow¬ 
ing  sections)  and  then  applied  directly  to  the  cathode  of  a  Tektronix 
type  561  oscilloscope  thereby  modulating  the  spot  intensity  of  the  CRT 
display.  Precision  potentiometers  connected  to  the  hydraulic  cylinder 
ram  and  the  lathe  carriage  provide  voltages  which  are  proportional  to 
receiver  position.  These  voltages  are  used  to  position  the  spot  of  the 
CRT  display,  the  modulated  spot  thereby  tracking  the  receiving  trans¬ 
ducer  position  at  reduced  scale.  The  resulting  intensity  modulated 
raster  is  time  photographed  and  then  photoreduced  approximately  5  to  1 
onto  high  resolution  glass  plates.  Visible  reconstructions  of  the 
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acoustically  "illuminated”  objects  are  then  obtained  from  the  reduced 
hologram  transparencies  in  the  usual  manner  by  passing  converging  or 
collimated,  coherent  light  from  a  laser  source  through  them. 

The  time  to  record  one  hologram  varied  between  four  and  nine  hours 
because  a  very  slow  scanning  rate  (5  cm/sec)  was  required  in  order  to 
minimize  vibrations  of  the  slender  probe  (2.5  cm  diameter  by  500  an 
long).  In  retrospect,  the  recording  process  could  have  been  speeded 
considerably  by  placing  the  hologram  apperture  near  the  surface,  thereby 
shortening  the  probe,  and  electrically  gating  the  signal  to  avoid  inter¬ 
ference  from  surface  induced  acoustic  reflections.  However,  this 
author  regards  mechanical  scanning  as  no  more  than  a  simple  and  in¬ 
expensive  expedient  for  gaining  information  directly  applicable  to  the 
design  of  holographic  systems  using  arrays  of  transducers,  possibly  a 
more  practical  approach,  but  certainly  a  much  greater  undertaking  in¬ 
volving  costly  apparatus.  Therefore,  the  lengthy  recording  time  was 
not  thought  to  be  other  than  a  nuisance. 

3.  ACOUSTIC  HOLOGRAPHY  WITH  ELECTRONIC  REFERENCE 

Imaging  by  holography  was  analyzed  in  Chapter  IV.  An  acoustic 
hologram  is  formed  by  interference  of  two  wavs  fields;  one,  A~(r), 
originates  from  the  illuminated  object,  the  other,  A^(r) ,  is  the 
reference. 

No  restrictions  are  placed  upon  the  object  illisaination  other 
than  adequate  spatial  coherence  (Section  II. 3)  and  that  the  object  be 
well  illuninated  (Section  III. 6).  However,  stringent  requirements  are 
placed  on  the  reference;  in  addition  to  spatial  coherence,  the  phase 


and  amplitude  variation  of  the  reference  wave  across  the  hologram  plane 

must  be  carefully  controlled  because  quality  reconstructions  are 

obtained  only  by  illuminating  the  recorded  and  developed  hologram  with 

an  accurate  replica  of  the  reference  (Section  IV.  4).  Simple  wave 

fronts,  such  as  aberration- free  spherical  or  planar  waves,  are  most 

easily  duplicated.  However,  within  the  limits  inpcsed  in  Section  IV. 4, 

the  reference  wave  may  be  as  complex  as  desired,  but  this  complexity 

47 

must  then  be  exactly  duplicated  in  the  reconstruction  illumination  . 

In  visible  optics  good  image  fidelity  is  obtained  for  either  case. 

In  acoustic  holography  although  in  principle  it  is  easy  to 
"illuminate"  the  object,  it  is  difficult  to  generate  true  spherical  or 
planar  acoustic  wavefronts  of  reference  quality;  i.e.,  without  severe 
aberrations  and  mode  effects  which  are  difficult  to  reprove  optically 
3nd  visual  reconstructions  from  such  holograms  say  be  severely 
distorted65.  However,  perfect  reference  wavefronts  may  be  simulated 
electronically  by  adding  the  amplified  output  E^fr)  of  the  receiving 
transducer  to  an  electronic  reference  signal  E^r)  which  is  obtained 
from  the  same  1  Mia  RF  source  which  drives  the  illisainating  transducer. 
That  is,  the  reference  is  derived  electronically  from  a  local 
oscillator  rather  than  acoustically  in  tie  tank.  Electronic 
simulation  cf  the  reference  is  possible  fceca?ise  piezoelectric 
acoustic  detectors  are  sensitive  to  instantaneous  amplitude  rather 
than  energy  (Section  II. 2). 

If  the  simulation  is  accomplished  sisply  by  adding  the  local 
oscillator  signal  directly  to  the  scattered  signal  without 
intervening  operations;  i.e.. 
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(V.  1) 


E^(r)  *  =  constant, 

the  coplanar  reference  is  obtained.  This  is  a  single  technique  which 
las  been  employed  by  other^6, 67,68 .  In  Figure  V.6  we  show  a  block 
diagram  of  the  apparatus  used  to  generate  coplanar  reference  holograms. 

To  simulate  the  inclined  planar  reference  requires  only  a  slight 
69 

complication  ;  the  reference  amplitude  remains  constant,  as  for  the 
coplanar  case,  but  the  reference  phase  varies  linearly  across  the  holo¬ 
gram  plane.  The  phase  of  a  planar  reference  wave,  propagating  with 
vector  wave  lumber  k  lying  in  the  xz  plane  and  making  angle  9  with  z, 
varies  across  the  hologram  plane  as  xk  sin  6.  To  simulate  this  wave 
electronically  the  phase  of  the  reference  signal,  derived  from  the 
local  oscillator,  must  vary  similarly  as  the  position  of  the  receiver 
changes.  However,  since  large  reference  angles  6  are  desired  in 
Leith-Upatnieks  holograms  sc  that  the  real  and  virtual  images  are 
completely  separated*6,  the  phase  variation  may  be  accomplished  in  a 
particularly  simple  way,  as  follows. 

In  order  to  uniquely  define  the  direction  of  the  reference,  the 
field  over  the  hologram  plane  must  be  ssspled  at  intervals  tx  at  least 
as  dose  as 


Lxk  sin  6  »s/2  .  (¥.2) 

However,  to  avoid  spurious  resolution  effects  the  raster  lines  rust  be 
separated  by  at  least  esse  receiver  effective  diameter  a.  Hence,  the 
aedass  reference  inclination  3s  consistent  with  both  criteria  is 


142 


%m3TCsmrn 


07-3) 


Haxisa  reference  inclinaticsi  for  scanned  holcgrsrs  is  plotted  in 
Figure  V.8  against  effective  receiver  disaster  expressed  in  wavelengths. 

In  our  experiments  the  receiver  «s  spfsrcxlsgteiy  two  scustic 
wavelengths  in  disaster.  Thus,  cssly  reference  inglinstic*  Engles  less 
than  about  8°  could  be  simulated  ursuhigucusly.  However,  this  provided 
adequate  separation  of  real  and  virtual  inages.  In  general  a  larger 


angle,  perhaps  30  ,  is  desirable;  but  the  sapli?f  resoiucisvs  nst  then 
be  less  than  an  acoustic  wavelength,  clearly  illustrating  the  is 
detector  resolution  required  by  holography  casar 
imaging  techniques. 

To  produce  scanned  acoustic  hologrscs  with 
at  B_  rerely  advance  or  retard  the  tisase  of  the  local  oscillator  signal 

•a 

added  to  da  scattered  signal  by  a/2,  relative  to  the  preceding  Use, 
for  each  sew  line  of  the  raster.  In  Figure  V.7  we  £sx  the  block  dia- 
grsa  or  the  apparatus  used  to  generate  inclined  reference  hoiograss  by 
this  technique.  The  required  phase  variation  was  easily 
with  sl  switch,  a  length  of  delay  line,  as d  the  polarity  inversion 
feature  of  an  csdllcsccpe. 

4.  SSRATISi  dr  THE  XSFQEiiSXKM  ST33S  IE3I 

Ha  key  eiesent  in  acoustic  is  fersatios  of  a  photo¬ 

graphic  trass parescy  of  bDlcgraa  ?(r)  patterned  after  t^e  acoustic 
standi^  wsve  field 


?(£>  « 


(7.4) 
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Figure  V.6.  Schematic  of  Copianar  Reference  Simulation 
0  1  Hiz  local  oscillator 

M  Modulator  and  power  amplifier  (50  W  RF  into  L) 

L  Acoustic  illumination  transducer 

AC  Acoustic  coupling  via  propagation  through  water 

R  Acoustic  receiving  transducer 

A  Amplifier 

S  Summing  network 

U  Electronic  processing 

C  Cathode  ray  tube  display 

T  Anechoic  water  tank 

g  Ground  through  water 
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Figure  V.7.  Schematic  of  Inclined  Reference  Simulation 

P  Pulse  generator  (0.8  msec  pulse  at  300  Hz  rr  typical) 
PD  Pulse  delay  (1.5  msec  typical) 

G  Gated  amplifier 

RS  Inclined  reference  simulation  network 

D  Delay  line  (0.25  usee) 


where  A^(r)  is  the  complex  propagating  acoustic  fieid  scattered  by  some 
coherently  illuminated  object  0  and  A^(r)  is  a  suitable  acoustic  refer¬ 
ence  field.  When  such  a  transparency  is  reilluminated  with  the  refer¬ 
ence  Apfr).  a  field  is  diffracted  by  F(r)  which  duplicates  A^fr)  in  some 
region,  hence  A^(r)  is  "reconstructed." 

The  procedure  by  which  the  informat ing  containing  term  is  obtained 
depends  upon  the  type  of  detector  used  to  form  the  hologram.  For  exam¬ 
ple,  when  the  detector  is  sensitive  to  field  intensity  rather  than 
amplitude 


P(r)  «  A(r)A  (r) 


(V.5) 


and  interferometric  techniques  are  required  to  produce  holograms.  To 
wit,  the  reference  and  object  fields  are  coherently  added  (hence  the 
strict  coherence  requirements)  and  then  applied  to  the  detector.  The 
output  of  the  detector  is 


PCD  «  A^(r)  +A^(r)  +2ReAQ(r) A^(r) 


(V.6) 


Ideally,  the  fi?n  processing  is  controlled  to  provide  perfect 
square  law  response.  Then,  if  either  the  reference  is  large 


4©  » 


(V.  7) 


or  the  object  is  diffusely  illuminated 


^CD 


constant 


(V.8) 
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Equation  01.6)  becomes 


P(r)  «  constant  +  2ReAQ(r)A^(r)  ,  (V, 9) 

the  desired  result. 

If  the  fields  are  acoustic,  amplitude  detection  is  possible  and 
generation  of  the  information  containing  term  is  not  restricted  to 
interferometric  techniques.  As  described  in  the  preceding  section,  the 
acoustic  fields  A^Cr)  and  A^(r)  are  first  converted  by  a  linear  process 
into  equivalent  electronic  signals  E^Cr)  and  E^Cr)  which  are  then 
combined  to  form  the  hologram.  If  they  are  added  and  their  sum  con¬ 
verted  by  electro-optical  means  into  a  photographic  transparency  accord¬ 
ing  to  a  square  law  or  at  least  power  law  transformation,  the  holograms 
thus  formed  are  equivalent  to  those  formed  from  visible  radiation. 

This  technique  was  introduced  in  the  preceding  section. 

Ihere  is  an  alternative  technique  which  offers  certain  advantages 
over  the  interferometric  method.  Instead  of  adding  the  object  and 
reference  signal  amplitudes  they  may  be  applied  to  a  modulator  which 
multiplies  chan.  This  technique  shall  be  called  heterodyne  detection. 

Consider  the  situation  presented  in  Figure  V.9  in  which  the 
object  and  reference  signals  are  multiplied  in  sane  sort  of  modulator  M. 
If  the  signals  are  monochromatic  (radial  frequency  w)  the  real  output 
E^Cr.t)  of  the  modulator  is 

*  EQ(r)  cos  [ut+^0(r) jE^(r)  cos  fwt+^fr)]  (V.10) 


where  Eq  and  E^  are  the  amplitudes  and  4>q  and  $R  the  phases  of  the  real 
signals.  Using  the  relation 


cos  a  cos  b  =  SSSiLt.  M  ♦  £°&  :  W 


CV.ll) 


the  modulator  output  becomes 


*  E0(r)ER(r)  cos(2ut  +  4>Q(r)  +  $R(r)]  + 
EqWErM  cos[4>0(r)  +  $>R(r)j  . 


(V.12) 


Clearly,  the  output  of  the  modulator  contains  a  d.c.  term  which  is 
identical  with  the  desired  standing  wave  term  (V.4^ 

The  second  harmonic  component  is  easily  eliminated  by  filtering. 
The  filtered  output  may  then  be  converted  into  a  transparency  according 
to  a  linear  transformation  such  that 

P(r)  «  constant  +  2&eA0(r)A^(r)  .  (V.13) 


Holograms  described  by  Equation  (V.13)  may  have  advantage  over  those 
described  by  Equation  (V.6).  First,  the  unwanted  A^  and  P^  terms  are 
automatically  eliminated,  regardless  of  the  ratio  A^/Aq  or  the  proper¬ 
ties  of  Aq.  Second,  greater  utilization  of  the  transparency  for  the 
storage  and  retrieval  of  object  information  is  obtained.  Third,  by 
performing  the  heterodyning  in  several  stages  high  signal  to  noise 
ratios  may  be  obtained. 


5. 


PHOTO-ACOUSTIC  CONVERSION  PROCESS 
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Optimization  of  hologram  quality  requires  that  certain  character¬ 
istics  of  the  hologram  forming  process  be  known,  particularly  the  de¬ 
pendence  of  amplitude  transmission  P  on  acoustic  field  amplitude  A  and 

* 

intensity  I  =  AA  .  Acoustic  hologram  formation  is  conveniently  sepa¬ 
rated  into  the  fan*  stages  shown  in  Figure  V.10:  (1)  detection  of  the 

acoustic  field  and  mixing  with  the  reference,  (2)  the  electro-optical 
conversion,  (3)  the  optical  to  photographic  conversion,  and  (4)  the 
scale  reduction.  Determination  of  the  composite  response  for  the 
acoustic  system  is  given  in  Appendix  E. 

6.  IMAGE  BRIGHTNESS 

An  acoustic  hologram  is  a  mosaic  of  alternating  light  and  dark 

fringes,  some  either  opaque  cr  clear;  others,  varying  shades  of  gray. 

* 

Ideally  P(r)  <*  A(r)A  (r),  but  the  less  restrictive  condition  that  con- 

* 

tours  P(r)  =  constant  be  geometrically  similar  to  contours  A(r)A  (r)  * 
constant  is  usually  sufficient  to  ensure  good  image  fidelity  in  recon¬ 
struction.  For  example,  by  increasing  fringe  contrast  it  is  often 
possible  to  brighten  the  image  without  causing  undue  distortion,  par¬ 
ticularly  important  when  the  hologram  aperture  is  small  as  in  the 
acoustic  case. 

In  order  to  gauge  image  brightness  for  interferometric  holograms 
*  2 

expand  P[A(r)A  (r)j  in  a  power  series  in  A  about  some  convenient  bias 

? 

or  ’’quiescent”  transmission  P(A£)  as  shown  in  Figure  V.ll  obtaining 
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Figure  V.10.  Schematic  of  the  Photo -Acoustic  Conversion  Process 
AC  Acoustic  signal  input  ;  PA  Optical  transparency  replica  output 
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Clearly,  issge  brightness  increases  isssotaaieally  with  visibility  as 
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multiplication  by  the  modulator.  Thus,  using  heterodyne  detection  very 
weak  signals  can  be  amplified  and  holograms  with  adequate  fringe  visi¬ 
bility  obtained. 

In  principle,  a  similar  effect  can  be  obtained  with  the  inter¬ 
ferometric  process.  However,  Gabor*  has  shown  that  interferometric 
magnification  is  limited  by  noise  in  the  large  reference  signal.  This 
noise  will  mask  extremely  weak  signals.  Furthermore,  in  order  to 
reduce  deleterious  effects  from  the  extraneous  product  term  A^,  the 
reference  amplitude  must  be  made  large.  In  so  doing,  the  camera  aper- 
ture  must  be  reduced  so  that  P(A^)  is  maintained  at  some  median  trans¬ 
mittance  lest  the  film  be  hopelessly  overexposed  by  the  reference.  How¬ 
ever,  closing  the  aperture  decreases  the  slope  of  the  characteristic 
curve,  hence  visibility  and  image  brightness  decrease. 

In  practice,  certain  departures  from  the  ideal  linear  situation 

described  above  were  found  to  be  desirable.  Two  considerations  influ- 

* 

enced  choice  of  the  functional  dependence  of  P  upon  AA  :  (1)  minimiza¬ 

tion  of  information  loss  in  the  reduction  stage,  ar.d  (2)  maximization  of 
image  brightness.  Resolution  of  fringes  in  the  reduced  transparency  is 
limited  by  the  demagnification  lens  and  is  maximum  for  high  contrast 
fringes.  However,  it  should  be  noted  that  it  was  found  possible  to 
maintain  a  gray  scale  with  fair  resolution  in  the  reduction. 

Maximum  image  brightness  demands  high  contrast  conversion  from 
acoustic  standing  wave  field  to  photographic  transparency.  However 
since  high  contrast  systems  have  limited  dynamic  range,  non-linear  con¬ 
version  in  the  form  of  clipping  will  usually  result.  For  interfero¬ 
metric  holograms  the  effects  of  non-linearity  are  minimized  without 
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sacrificing  contrast  by  reducing  the  operating  point  of  the  character¬ 
istic  curve  to  zero  transmission,  as  shown  in  Figure  V.ll. 

Over  most  of  the  hologram  plane  the  reference  to  sign  a?  amplitude 
ratio  Aj^/Aq  satisfies  A^/Aq  >  2,  at  the  signal  maxima  A^/Aq  >  1.  Thus, 
in  regions  of  signal  maxima  the  darkest  portions  of  the  fringes  will  be 

clipped.  However,  such  clipping  produces  negligible  image  distortion 
* 

because  P(AA  )  is  small  in  these  regions  even  if  linearity  is  strictly 
maintained. 

7.  EXPERIMENTS  IN  ACOUSTIC  IMAGING  BY  HOLOGRAfflY 

The  following  acoustic  imaging  experiments  were  performed  using 
the  techniques  and  apparatus  described  above  and  in  Appendices  D,  E, 
and  F. 

HOLOGRAM  #1  (Figures  V.12-V.16) 

The  simplest  possible  scene  was  chosen  for  hologram  number  one. 
The  scene  (Figure  V.12)  is  a  view  of  the  acoustic  radiation  pattern 
from  the  illumination  transducer  (see  Figures  D.l,  D.2,  D.3)  without 
any  intervening  object  and  with  electronically  simulated  coplanar 
reference.  Unfortunately,  a  3  inch  diameter  airfilled  standpipe 
protrudes  from  the  bottom  of  the  tank  inot  the  object  space  in  front  of 
the  hologram  plane.  However,  since  the  central  lobe  of  the  radiation 
pattern  diverges  narrowly  (only  about  10°)  it  was  anticipated  that  the 
standpipe  would  not  appear  in  the  reconstruction.  The  hologram  is 
shown  in  Figure  V.13  and  the  reconstruction  in  Figure  V.14. 

The  central  spot  in  the  reconstruction  is  clearly  the  image  of 
the  transducer,  approximately  a  point  source.  Occurrence  of  a  second 
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Figure  V.14.  Reconstruction  from  Hologram  #1 

Left  half  of  reconstruction  has  been,  cropped.  Spot  at  left 
corresponds  to  location  of  standpipe  SP  and  transducer  L. 
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spot  at  the  right  band  side  of  the  reconstruction  is  somewhat  perplex¬ 
ing.  (There  is  also  a  spot  on  the  left,  symmetrical  with  the  one  on 
the  right,  but  somewhat  weaker.  This  spot  is  not  shown  in  V.14.)  By 
selectively  illuminating  small  areas  of  the  hologram  at  a  time  it  was 
determined  that  the  spurious  images  arise  from  the  fringe  structure 
near  the  left  and  right  hand  edges  of  the  hologram. 

Two  possible  explanations  for  the  occurrence  of  the  edge  fringe 
structure  are  offered;  one  is  based  upon  Moire  fringe  interference, 
the  other  upon  multiple  reflections.  Of  the  two,  the  Moire  theory  is 
more  likely.  If  an  equally  spaced  line  grating  (Figure  V.lSa)  is 
superimposed  upon  a  symmetrical  Fresnel  zone  (Figure  V.l5b)  then, 
according  to  Moire  fringe  theory  ,  a  new  fringe  system  like  the  one 
showr,.  in  Figure  V.15c  will  be  formed.  Note  the  appearance  of 
symnetrical  Fresnel  zones  near  the  edges  of  the  pattern  in  addition  to 
the  original  zone  in  the  center.  These  spurious  zones  will  form 
diffraction  spots  to  the  right  and  left  of  the  spot  from  the  original 
zone.  Thus,  by  identifying  the  line  grating  with  the  scanning  raster 
and  the  central  Fresnel  zone  with  interference  between  the  transducer 
radiation  and  the  coplanar  reference  it  is  clear  that  the  spurious 
images  may  be  a  result  of  coarse  scanning. 

One  may  also  explain  the  existence  of  spurious  fringes  at  the 
edge  of  the  hologram  by  assuming  that  a  portion  of  the  acoustic  beam 
is  reflected  from  the  standpipe  to  both  walls  of  the  tank  and  then 
from  the  walls  to  the  hologram  plane  as  shown  by  the  dotted  rays  in 
Figure  V.  12.  These  doubly  reflected  rays  will  form  fringepattems 
similar  to  those  found  in  the  hologram.  Of  course,  if  either  the 
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(a)  line  grating 


(b)  Fresnel  zone 


(c)  Moire  fringe  system 

Figure  V.15.  Effect  of  Scanned  Sampling 

Superimposing  a  line  grating  upon  a  Fresnel  zone  yields  a 
spurious  Moire  fringe  pattern. 


walls  or  the  standpipe  is  non-reflecting  then  this  theory  does  not 
apply. 

Since  the  standpipe  is  filled  with  air  it  should  be  a  strong 
reflector  for  acoustic  radiation.  The  wall  is  a  different  matter.  The 
walls  of  the  chamber  are  composed  of  a  rubber  cone  structure  overlaid 
upon  redwood  paneling  attached  to  a  reinforced  concrete  foundation. 

The  cones,  approximately  2  cm  high  by  2  cm  base  diameter,  form  a  regular 
close  packed  square  array.  The  reflected  amplitude  of  normally  incident 
1  MHz  acoustic  radiation  was  found  to  be  15  db  below  incident  amplitude. 
However,  at  sufficiently  oblique  incidence  the  reflection  may  be  con¬ 
siderably  greater.  This  was  not  investigated. 

The  receiver  characteristics  introduce  additional  complications 
which  further  obscure  the  exact  origin  of  the  spurious  images  in  recon¬ 
struction  #1.  First,  cw  illumination,  was  used  during  the  early  holo¬ 
gram  experiments,  hence  coincidence  techniques  (range  gating)  could  not 
be  used  to  separate  reflections  from  direct  radiation.  Second,  an 
extremely  sensitive  non-linear  detector  was  employed  in  order  to 
maximize  the  number  of  fringes  recorded.  The  output  of  this  detector 
was  the  same  for  weak  secondary  reflections  as  for  strong  primary 
radiation.  Thus  fringe  contrast  could  not  be  used  to  separate  the  two. 

The  detector  departed  from  the  simple  linear  amplifier  system 
shown  in  Figure  V.6  because  it  was  felt  initially  that  the  number  of 
interference  fringes  recorded  must  be  maximized  in  order  to  insure 
adequate  resolution  and  signal  to  noise  ratio.  This  was  accomplished 

4 

by  passing  the  signal  through  a  high  gain  (10  )  amplifier  stage 
followed  by  symmetrical  clipping  before  going  to  the  CRT  grid.  See 
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Figure  V.16a.  The  large  amplification  insures  that  even  very  weak 
signals  will  produce  sufficient  signal  at  the  CRT  grid  to  turn  the  beam 
on,  clipping  prevents  overexpose  and  the  concomitant  loss  of  resolu¬ 
tion  and  fringe  contrast  with  signals  of  large  dynamic  range,  and  the 
symmetrical  design  avoids  problems  with  d.c.  restoration. 

The  effect  of  clipping  upon  the  transmission  contour  of  typical 
interference  fringes  in  illustrated  in  Figure  V.16a.  All  fringe  ampli- 

are  clipped  so  that  A  <  A^,.  Thus  amplitude 
information  is  no  longer  encoded  as  variations  of  fringe  contrast  since 
fringe  contrast  is  essentially  a  constant.  Amplitude  also  controls 
fringe  width  W,  as  may  be  seen  in  Figure  V.16a;  however,  this  is  a 
non-linear  dependence  (see  figure  V.16b)  so  that  basically  amplitude 
information  is  not  faithfully  recorded  in  this  type  of  hologram.  Non¬ 
linear  processes  of  this  type  are  useful  for  high  contrast  objects. 
However,  when  a  gray  scale  must  be  maintained  -  as  with  imaging  of 
biological  tissue  -  then  linearity  nay  be  required. 

HOLOGRAMS  #2  THROUGH  #6  (Figures  V.17-V.1S) 

The  holograms  of  this  set  are  almost  self-explanatory.  The 
object  scene  for  this  group  represents  the  next  logical  step  up  in 
complexity  from  the  object  scene  for  hologram  #1.  The  scene,  depicted 
schematically  in  Figure  V.17,  consists  of  an  acoustically  opaque  screen 
(3mm  thick  aluminum  sheet  bonded  on  both  sides  with  6  mm  thick  sponge 
rubber)  with  a  circular  hole  14.3  mm  diameter  in  the  center  illuminated 
from  behind.  Simple  transmission  objects  of  this  type  are  ideal  for 
evaluating  various  acoustic  holographic  techniques  because  the  hologram 


tudes  greater  than 


Figure  V,.  18(a).  Hologram  #2,  0°  Simulated  Reference  Angle 

Regions  of  high  acoustic  field  are  white. 

Factor  of  2  (m  =  2)  smaller  than  original  acoustic  field 


Phase  shift  per  line  =  tt. 


(2.8  mn)  and  the  reference  is  shifted  by  tt/2  for  each  new  line. 

Clearly,  the  ambiguity  in  Figure  V.18c  is  a  result  of  the  finite 
resolution  (.175  lines/nm)  of  the  acoustic  sensitive  "film."  This 
resolution  corresponds  to  a  photographic  emulsion  able  to  resolve  about 
420  lines/mm,  which  is  generally  considered  too  low  for  planar  refer¬ 
ence  visible  holography.  Figure  V.18e  was  obtained  using  the  same  line 
spacing  as  for  V.18d  except  that  the  phase  shift  was  m  between  lines 
rather  than  tt/2.  This  is  equivalent  to  doubling  the  reference  angle, 
however,  it  is  clear  that  ambiguity  is  again  introduced. 

Reconstructions  from  the  various  holograms  are  displayed  in 
Figure  V.19.  In  each  case  the  very  bright  spot  represents  an  image  of 
the  backlighted  hole.  It  will  be  interesting  to  compare  the  size  of 
the  spot  appearing  in  the  reconstruction  with  the  size  of  the  original 
hole.  The  hologram  transparency  is  1/26.4  the  size  of  the  original 
acoustic  field;  the  reconstruction.  Figure  V.19b  is  magnified  approxi¬ 
mately  12.  Thus  the  reconstructed  spot  should  measure  approximately 
7  ran  diameter.  Measurement  of  the  reconstruction  indicates  that  within 
experimental  error  this  is  indeed  the  case. 


In  these  experiments  linear  rather  than  clipped  detection  was 
employed  thus  only  the  central  diffraction  zone  of  the  hole  contributes 
to  the  hologram,  the  higher  zones  being  below  the  threshold  of  the  CRT 
display.  However,  it  is  important  to  note  that  approximately  correct 
spot  size  is  obtained  even  so.  The  higher  order  diffraction  rings  were 
not  required  in  order  to  obtain  proper  resolution  of  this  hole.  This 
result  is  in  agreement  with  visible  holography  experiments  quoted  by 
Develis  and  Reynolds.  If  clipped  detection  had  been  employed  then 


Figure  V.19(a).  Reconstruction  from  Hologram  #4 

Reconstruction  photographed  on  Polaroid  5S  PA  film  with 
.6328  micron  laser  illumination.  Positive  print  made  on 
Eastman  Kodak  Kodabromide  F-5  single  weight  paper. 


Figure  V- 19(b).  Reconstruction  from  Hologram  #5 

Zero  order  is  faint  area  in  center,  reconstruction  is 
bright  spot  near  right  edge. 
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Figure  V.  19(c).  Enlargement  of  Reconstruction  from  Hologram  #5 

Figures  V.19  (c) ,  (d),  and  (e)  demonstrate  that  recon¬ 
struction  is  at  least  20  times  brighter  than  background. 


Figure  V.  19(d).  Enlargement  of  Reconstruction  from  Hologram  #5 
Enlargement  exposure  5  times  less  than  for  Figure  V.  19(c). 


Figure  V.  19(g).  Reconstruction  from  Holograms  If 3  and  8 4 

TV  at  maximum  contrast,  minimum  brightness,  illustrating 
utility  of  'IV  displays  for  viewing  acoustic  holograms. 
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Figure  V.  19(i).  Reconstruction  from  Holograms  #2  and  #5 

TJ  at  maximum  contrast,  minimum  brightness,  illustrating 
utility  of  TV  displays  for  viewing  acoustic  holograms. 
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the  higher  order  diffraction  rings  would  also  have  contributed  to  the 

hologram  and  the  signal  to  noise  ratio  might  have  been  better  but  an 

erroneous  measure  of  the  hole  size  might  have  been  obtained.  Clipping 

appears  to  have  sane  relation  to  the  apodization  techniques  of  con- 
71 

ventional  optics  . 

HOLOGRAMS  #7  THR0UQ1  #10  (Figures  V.2G-V.28) 

Transmission  through  thin  planar  objects  is  useful  for  evaluating 
techniques  and  for  visual  displays.  However,  practical  applications 
particularlv  in  acoustic  imaging,  require  that  the  hologram  be  made  by 
acoustic  radiation  reflected  fran  the  object  scene  rather  than  trans¬ 
mitted  through  the  object  scene.  Thus,  the  remainder  of  the  holograms 
to  be  presented  here  are  of  reflecting  object  scenes  even  though  such 
scenes  are  considerably  more  difficult  to  image  primarily  because  of 
the  predominately  specular  nature  of  the  acoustic  reflections  from  most 
objects.  (See  Chapter  III.) 

Holograms  #7  through  #10  were  the  first  of  the  holograms  made  by 
reflected  acoustic  radiation;  the  objects  were  simple  reflecting  laminae 
fabricated  fran  steel  plate  (#7  and  #8)  or  aluminum  plate  (#9  and  #10). 
Holograms  #7  and  #8  (Figures  V.21  and  V.22)  were  made  with  clipped 
detection  and  cw  illumination.  The  reconstruction  is  shown  in  Figure 
V.23. 

Holograms  #9  and  #10  (Figures  V.25  and  V.26)  were  made  with 
linear  detection  and  pulsed  illumination.  The  reconstructions  are 
shown  in  Figure  V.27.  Pulsed  illumination  was  required  in  order  to 
prevent  spurious  signals  radiated  by  the  illumination  transducer  from 
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Figure  V.20.  Layout  of  Object  Scene,  Holograms 
0  Object,  15  on  x  20  an  x  .1  cm  thick  s 


Regions  of  high  acoustic  field  are  white. 

Bright  spot  in  center  of  hologram  is  the  central  lobe  of 
transducer  L  reflected  from  the  object  0. 

Raster  line  spacing  =  5.7cm,  230  lines  total. 
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Figure  V. 22.  Hologran  #3,  7.7°  Sibilated  Reference  Angie 

Raster  line  spacing  =  5.7  zxs. 

Phase  snift  per  line  = 
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Figure  V. 23(b).  Recons true t i on  fret.  Hslogras  #7 


Viewed  directly  by  TV  at  sedius  contrast  and  brightness. 
Photograph  of  TV  monitor  on  Polaroid  46L  fils,  positive 
print  on  Eastman  Kodak  Kodabrcnide  F-5  single  weight  paper. 
Figures  V. 25(a) -(c)  illustrate  the  ’utility  of  viewing 
acoustic  holograss  with  television  svsts®.  The  variable 
contrast,  brightness,  and  threshold  of  these  systems  allow 
instantaneous  and  continuous  adjustment  of  the  display  for 
op  tin  25  appearance  of  the  ioage  against  the  high  noise 
background  characteristic  of  acoustic  feolcgrsss.  The  jitt 
inherent  in  TV  rasters  rethices  bothersome  stickling  vithou 
serious  loss  of  resolution. 
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Figure  V.  25(c).  Reconstruction  free  Kologrss  #7 

Viewed  directly  by  TV  at  saxisaa  contrast,  siniss  brigh 
ness,  clearly  cesonstTating  mat  tns  recons true  tea  i^age 
is  actually  tiich  brighter  than  background.  Viewing 
directly  by  eye  or  contirticas  terse  photography  obscures 
this  fact. 
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V.24,  Layout  of  Object  Scene,  Holograms  #9  and  #10 
>  Object,  20  an  *  22,  !>  an  k  .2  cm  thick  bolted  to  fixed  support  tube  S 


being  recorded  on  the  hologram.  Spurious  isolation  was  eliminated  by- 
range  gating  the  received  acoustic  pulses;  a  typical  pulse  reflected  by 
the  target  is  shown  in  Figure  V.28.  The  structure  at  the  pulse  edges  is 
due  to  the  difference  in  path  length  between  the  front  and  back  of  the 
transducer  and  between  the  various  positions  on  the  object  and  hologram 
plane.  Care  was  taken  to  adjust  the  coincidence  circuitry  so  that  only 
a  narrow*  portion  from  the  center  of  the  acoustic  pulse  was  used  to  fore 
a  hologram. 

The  triangular  shaped  notch  used  as  the  object  for  holograms  #9 
and  #10  was  picked  so  that  some  measure  of  the  resolution  of  the  system 
could  be  obtained.  The  object  was  located  approximately  1.7  m  from  the 
hologram  plane  and  it  appears  that  object  points  separated  greater  than 
12  ran  could  be  resolved.  This  corresponds  to  3.6  milliradians  angu7 r<r 
resolution.  Conversely,  since  the  hologram  plane  measured  approxi¬ 
mately  79  *  63  cm  the  Rayleigh  criterion  predicts  that  the  angular 
resolution  could  be  no  better  than  2.3  milliradians-  Thus,  the 
resolution  seems  reasonable  for  the  aperture  available.^ 

HOLOGRAMS  #11  AND  #12  (Figures  V.29-V.31) 

The  next  logical  step  in  an  elementary  study  of  acoustic  holo¬ 
graphy  is  to  image  objects  with  some  measure  of  three-dimensionality, 
the  simplest  of  these  being  multiplanar  arrays.  In  Figure  V.29  we  show 
an  arrangement  of  two  laminae,  one  square,  the  other  circular,  situated 
in  non-parallel  planes  with  the  square  approximately  9  cm  behind  the 
circle.  In  Figure  V.30  we  show  a  conventional  coplanar  reference 
acoustic  hologram  (#11)  of  this  scene.  In  Figure  V.31  we  shew  a 


(a)  In  tipper  half  of  hologram. 


(b)  In  lower  half  of  hologram. 


Figure  V.28.  Typical  Acoustic  Pulses  Scattered  by  an  Object 

Horizontal  scale  (time)  =  .2  msec/division . 

Vertical  scale  (voltage  input  to  amplifier  G)  =  .5  volts/div 

Note  that  the  acoustic  signal  requires  approximately  .2msec 
for  rise  and  fall  due  to  difference  in  path  length  from 
various  portions  of  the  object  to  the  hologram  plane.  Only 
center  .4  msec  of  pulse  was  used  for  the  holograms. 


Figure  V. 29. 


M _ Mill 


hologras  (#12)  identical  Kith  hoiGgrss  #11  sccept  that  Nereis  #11  was 
generated  using  conveaticgal  interference  seteciits,  #12  «s  generated 
using  the  heterodyne  detection  scheme  outlined  is  Section  V.4.  Sr 
asparing  the  two  hslsgrss  it  is  cleat  Sat  they  are  equivalent  except 
Sat  the  fringe  contrast  of  #12  is  less  than  #11.  However,  fringe 
contrast  is  easily  adjusted  by  changing  system  gain  and  the  zero  signal 
intensity  of  the  CHf  spot. 

The  reconstructions  frog  holograms  #11  and  #12  are  rather  poor 
for  several  reasons.  First,  because  of  the  narrow  bean  angle  or  the 
acoustic  illumination  it  was  not  possible  to  properly  illuminate  the 
edges  of  the  laminae,  hence  the  edges  were  not  well  defined  in  the 
reconstruction.  Moreover,  because  of  the  narrow  best  angle  it  was  not 
possible  to  separate  the  laminae  further  hence  the  depth  of  focus  of  the 
hologram  overlaps  the  depth  of  field  of  the  scene  and  the  two  objects 
appear  coplanar  in  the  reconstruction.  Finally,  it  appears  that  for 
the  small  apertures  of  acoustic  hologrsss  one  say  wish  to  sacrifice 
accurate  linear  display  of  the  hologram  fringes  for  a  binary  display 
which  improves  signal  to  noise  ratio  at  the  expense  of  a  completely 
true  gray  scale  rendition  of  the  object. 

It  is  interesting  to  note  that  the  difficulties  of  specular 
reflection  are  already  clearly  evident  even  in  simple  situations  like 
Figure  V.29.  Note  that  the  acoustic  field  in  the  hologram  plane  is 
composed  of  essentially  two  ‘bright  patches".  That  on  the  left  is  due 
primarily  to  reflection  from  the  rank,  object  element,  that  on  the 
right  is  due  primarily  to  reflection  from  the  square  object  element. 

Now  the  angle  of  inclination  between  the  two  object  planes  was  only 


S’,  However,  if  sngle  had  bees  increased  it  ‘«gi!c  then  aot  have 
been  possible  to  cents  in  fce'h  natives  sinui  reneges  iy  within  the 
hologram  aperture,  hence  it  would  not  have  been  possible  to 
sissiltmaousiy  inage  both  coject  deserts  . 


HHJDSSS  *15  .A3©  #14  {Fibres  Y-32-Y.35) 

*e  have  presented,  in  the  foregoing,  acoustic  hole-grans  of 


various  object  scenes  of  progress  ibely  none  complicated  structure. 


inis  procedure  was  adopted  as  the  nest  logical  way  in  vhich  tc  conduct 


an  elementary  study  of  acoustic  inaging-  By  beginning  with  extrsaely 


sirple  objects,  accurate  reasuresents  can  be  ccnpared  with  theoretical 
calculations  without  midue  difficulty-  Then,  the  special  problems 
posed  by  core  realistic,  but  also  sore  complicated,  scenes  can  be 
readily  isolated  a:  t  studied.  We  conclude  this  logical  progression 
with  holograms  213  and  *14,  holograms  of  a  very  general  object  scene 
representative  of  those  likely  to  be  encountered  in  practical  undersea 
acoustic  imaging  situations.  These  holograms  illustrate  a  serious 
obstacle  to  useful  acoustic  imaging  namely  the  specular  nature  of 
acoustic  reflection  from  cost  objects.  (See  Chapter  III.) 

The  physical  arrangement  for  holograms  *13  and  *14  is  shown  in 
Figure  V.32.  The  object,  which  is  viewed  in  reflection,  is  a 
cylindrical  underwater  light  globe  housing  with  one  end  spherical  and 
the  other  end  capped  by  a  flat,  bronze  fitting.  Hologram  #13  is  shown 
in  Figure  V.33,  and  the  reconstruction  in  Figure  V.34.  Unlike  the 
preceeding  holograms,  instead  of  a  complete  image  of  the  object  only  a 
diffracted  highlight  is  observed.  This  is  a  characteristic  displayed 
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Figure  V.33.  Hologram  #13,  0°  Simulated  Reference  Angle 

Fringe  system  at  left  due  to  cylinder,  fringe  system  at 
right  due  to  spherical  end  cap. 
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by  most  objects  illuminated  with  acoustic  radiation.  The  explanation 
is  quite  simple. 

In  Chapter  III  the  role  of  surface  texture  in  determining  the 
acoustic  scattering  properties  of  material  surfaces  was  discussed.  For 
most  surfaces,  the  characteristic  depth  of  the  texture  is  considerably 
less  than  X  (X  =  1.5  mm  for  1  MHz  acoustic  radiation  in  water),  the 

V 

characteristic  correlation  length  of  the  texture  much  greater  than  X. 
Accordingly,  most  objects  observed  by  acoustic  means  appear  to  be 
smooth;  i.e.,  the  acoustic  scattering  is  specular  rather  than  diffuse. 

In  Chapter  III  we  also  discussed  the  consequences  of  viewing 
through  small  apertures  specularly  reflecting  curved  surfaces  illumi¬ 
nated  with  plane  or  spherical  waves;  only  the  central  lobe  of  the 
diffraction  pattern  will  be  observed.  Simple  geometrical  arguements 
show  that  this  lobe  arises  mainly  from  that  portion  of  the  reflecting 
surface  directly  opposite  the  center  of  the  viewing  aperture,  hence 
the  object  appears  as  a  diffracted  highlight.  The  glass  housing 
(Figure  V.35a)  used  as  the  object  for  hologram  #13  is  acoustically 
smooth  in  the  sense  outlined  above,  hence  ony  would  expect  that  the 
reconstruction  from  hologram  #13  would  show  only  a  diffracted  high¬ 
light.  This  is  indeed  the  case. 

To  circumvent  the  highlight  problem  one  usually  rescrts  to 
diffuse  illumination.  However,  generation  of  diffuse  but  coherent 
acoustic  radiation  is  not  a  simple  task.  Instead  we  choose  to  make 
the  object  quasi-diffusely  reflecting  by  attaching  to  it  airfilled 
glass  spheres  approximately  2.5  cm  diameter  as  shown  in  Figure  V. 35(b). 
This  artifice  simulates  a  surface  texture  with  a  characteristic  depth 
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of  several  wavelengths. 

The  modified  glass  housing  was  replaced  in  exactly  the  same 
position  as  for  hologram  #13  and  a  hologram  was  then  made  of  the  new 
object  scene  using  the  same  hologram  parameters  as  before.  The  new 
hologram  (#14)  is  shown  in  Figure  V.36.  Note  that  it  bears  a  certain 
resemblance  to  #13  although  many  of  the  fringes  are  now  considerably 
mottled  or  even  washed  out  completely.  Of  course,  this  is  precisely 
what  one  would  expect  as  a  specular  object  is  gradually  "roughened" 
into  a  diffuse  object.  Many  areas  of  the  hologram  exhibit  very  fine 
fringe  structures  again  indicitive  of  diffuse  reflection. 

Because  the  resolution  of  the  CRT  display  was  limited,  the  fringe 
detail  in  the  physical  hologram  plane  was  inadequately  preserved  in  the 
recorded  hologram.  By  watching  the  spot  during  a  typical  scan  it  was 
apparent  that  there  was  definite,  high  contrast  fringe  detail  in 
the  hologram  at  spacings  as  small  as  one  fifth  the  resolvable  spacing 
of  the  CRT.  To  obtain  the  required  display  resolution  would  require 
the  use  of  a  magnetically  focused  CRT.  Thus,  a  useful  reconstruction 
from  hologram  #14  was  not  possible.  However,  simply  by  ccsnparing 
holograms  #13  and  #14  it  is  apparent  that  diffuse  acoustic 
illumination  of  rough  surfaces  will  indeed  allow  imaging  of  the  entire 
object,  not  just  the  highlights. 

8.  SENARY 

A  simple  scanning  technique  for  generating  acoustic  holograms  of 
underwater  objects  in  the  laboratory  has  been  described.  Using  this 
system  acoustic  holograms  have  been  recorded  which  show  angular 
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resolution  of  3.6  milliradians,  approximately  1.5  times  the  Rayleigh 
limit.  These  experiments  were  performed  at  1  f-fiz,  a  frequency  suitable 
for  long  range  underwater  imaging  and  detection. 

Because  instantaneous  amplitude  is  an  acoustic  observable,  the 
reference  field  can  be  simulated  electronically.  Moreover,  by  resorting 
to  heterodyne  or  phase  detection  the  cross  product  term  between  object 
and  reference  signals  may  be  generated  without  the  undesired 
extraneous  terms  which  occur  in  conventional  holography.  Holograms  of 
the  Leith-Upatnieks  type,  with  the  inclined  planar  reference  simulated 
electronically,  and  heterodyne  detection  holograms  are  presented. 

By  using  a  CRT  for  the  acousto-optic  conversion  it  is  possible  to 
synthesize  electronically  a  wide  variety  of  conversion  characteristics 
from  linear  to  hard  limiting.  Because  of  the  rather  small  apertures 
involved  it  appears  that  hard  limiting  (high  fringe  contrast)  or 
clipping  is  often  desirable  in  order  to  maximize  reconstructed  image 
brightness. 

Various  other  aspects  of  acoustic  imaging  are  presented.  In 
particular,  the  severity  of  specular  acoustic  reflection  is  clearly 
evident  in  these  experiments . 

Finally,  the  utility  of  variable  contrast  television  displays  for 
viewing  acoustic  holograms  is  demonstrated. 
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APPENDIX  A 

FOURIER  TRANSFORM  CONVENTIONS 


In  optics  we  shall  have  occasion  to  calculate  Spatial  Fourier 
Transforms.  Consider  the  general  function  of  position  A(r)  where 


r  =  xx  +  yy  +  zz 


(A.  la) 


and  x,  y,  z  are  cartesian  coordinate  unit  vectors  with  z  lying  along 
the  general  direction  of  propagation.  Position  in  a  specific  trans¬ 
verse  plane  ?  =  z^  =  const,  is  denoted  where 


Ic  =  xx  +  yy  +  z,z 


CA.lbl 


The  Spatial  Fourier  Transform  A(r^)  of  A(r^.)  is  taken  as  the 
two-dimensional  Fourier  Transform  in  the  transverse  plane  z  *  z^. 


ACrJ)  =  A(x,zc)  =  JJ  ACr^e 


■/ 


-Znir^s  7 
AC^e  *  -d2rc 


(A.2) 


with  the  inverse 


A(s,zc)  - 


=  j  A(s,zc) 


2lTi  S’T-  7 

}e  “  “Vs  =  AQ^) 


(A.  3) 


Here  J  denotes  integration  over  the  entire  transverse  plane  z  -  z^7 
C 

J  denotes  integration  over  the  entire  £  plane  and 


v  +  y 


(A.  4) 


is  the  spatial  frequency  in  cycles/ length.  In  general,  the  Spatial 
Fourier  Transform  depends  upon  which  plane  z  =  const,  is  chosen  for  the 
transformation;  this  z  dependence  is  explicitly  stated.  For  example, 
in  the  plane  z  =  0  the  Spatial  Fourier  Transfora  of  A(r>  is  denoted 
A(s,0). 

Temporal  Fourier  Transforms  are  also  required.  The  Temporal 
Fourier  Transform  A(r,t)  of  a  function  of  time  and  position  A(r,t)  is 
taken  as 


ACryt)  =  A(r,v)  =  J  A(r,t}e"2~1Vtdt 


(A.  5) 


with  the  inverse 


A(r,v)  =  J~  A  (r ,  v)  e2"  ivtdv  =  A(r;t) 


(A.6) 
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and  v  is  the  temporal  frequency  in  cycles/sec.  In  general,  the 
Temporal  Fourier  Transform  depends  upon  position;  this  r  dependence  is 
explicitly  stated.  Positional  dependence  will  always  be  denoted  prior 
to  temporal  dependence. 
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APPEXDEC  B 


coK^HiincK  a?©  jgmMm  ussss 


The  convolution  A*B  of  ts*o  functions  A{x)  and  3(x)  is  tsusES  as 


A(x)  *  B(x)  =  J  ACu)3(x-u)du 


fs  r» 


The  correlation  A#3  of  tvo  functions  A(x)  sad  B(x)  is  tarm  ss 


A(x)  ★  B(x)  =  J  A{u)3(xnj)d3 


(B.2) 


In  cities  two-diggnsicaal  convolutions  and  correlations  are  also 
needed.  They  are  taken  as 


aid 


*  ” 

A(x.y)  *  S(x,y)  *  /^ACu.v)B(x-ii,y**3© 


A(x,y)*3(x,y)  =  JJ  .\{usv}B{x^^< 


05.3} 


CB.*3 


For  functions  of  position  in  transverse  planes 
convolution  and  corr elation  bears 

s 

AQr^}  *  B(tv)  *  y  ACr-lBC^-iyi^ 


(B-5) 


sns 


* 

^  3v£^J  *  ^*AC*w5».y^is  *0 


(B.S} 


Convolution  and  correlation  in  both  space  and  time  are  encountered. 

If  the  meaning  is  not  clear  then  *r  and  ^  will  denote  spatial  opera 
tions,  and  ^  will  denote  temporal  operations,  and  *ks  will 
denote  spatial  frequency  operations,  *v  and  will  denote  temporal 
frequency  operations. 

The  following  properties  of  convolutions  will  prove  useful. 


A(x)  *  B(x)  -  A(s)B(s) 


(B.7) 


A(x)  *  B  (x)  =  A(5)B  C-s) 
AQ^B^x)  =  A(s)B*(s) 


(B.8) 

(B.9) 


A(s)  *  B(s)  -  A(x)B(x)  . 


(B.10) 


The  following  properties  of  correlations  will  prove  useful. 


A(x)^rB(x)  -  A(-s)B(s) 


A(x)^B*(x)  =  A(-s)B  (-s) 


(B.ll) 

(B.12) 


A(-x)  ★  B(x)  =  A(s)B(s) 


(B.13) 


A(s)  ★  B(s)  *  A(-x)B(x) 


(B.14) 


In  addition,  convolutions  and  correlations  are  distributive, 
associative,  and  comnutative.  From  the  preceding  relations  we  obtain 


|A(x)  | 2  =  A(s)  *  A* C-s)  -  A(-s)  "k  B*(-s) 


(B.15) 


APPENDIX  C 


DECOMPOSITION  INTO  PLANE  HAVE  SPECTRUM 


Any  field  5(r^)  satisfying  the  Helmholtz  equation  can  be  decom- 
posed, into  a  spectrum  of  plane  waves  thus  “ 


ao 

CCic)  «  5(p,zc)  =  / 


i[K-rc+kz(K)zc]  2 


(C.l) 


where  W(k)  is  a  weighting  function  giving  the  amplitude  of  the  kth 
directed  plane  wave  component  of  5(1^.)  and  jc  is  the  transverse  wave 


number 


k  =  k  +  kzz 


,2  2  ,2 

1 C  =  k  +  If 


(C.2) 


Alternatively,  the  field  can  be  decomposed  into  a  spectrum  of  spatial 
frequency  components  by  taking  a  Spatial  Fourier  Transform  thus 


5(rc)  ■  j  Sts. 


2»i£-2c,2 

zc)e  d  s 


(C.3) 


where 


s  =  sxx  +  Syy 


(C,4) 


represents  a  spatial  frequency  (number  of  periodic  variations  of  any 
physical  feature  per  unit  length) ,  the  direction  of  £  indicating  the 
orientation  of  the  spatial  variations  in  the  transverse  plane. 
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Equating  (C.l)  with  (C.3)  demands 


K  -  2rrs 


CC.5) 


|(s,zc)  =  We 


o  ik  (2rrs)zr 
ue  2  SffCais) 


(C.6) 


Evidently 


5(s,0)  =  (2w)2K(2rr5) 


i(s,zc)  =  £(s,0)e 


ik^C^TfS) 


(C.7) 


Thus,  the  plane  wave  spectrum  W(<)  is  preserved  in  propagating 
between  transverse  planes  but  the  spatial  frequency  spectrum  is  not. 
Even  so,  it  is  usually  more  convenient  to  obtain  the  plane  wave  spec¬ 
trum  from  the  spatial  frequency  spectrum  rather  than  the  converse.  The 
reason  is  that  the  fields  of  interest  usually  arise  from  sources 
located  on  boundaries  of  known  geometry.  These  boundary  fields  are 
easily  decomposed  into  known  spatial  spectra  to  which  the  resulting 
plane  wave  spectrum  is  simply  related.  Thus,  writing 
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the  plane  wave  spectral  decomposition  becomes 
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Substituting 
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into  Equation  (C.9)  we  obtain 
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Finally 
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It  is  inviting  to  consider  the  equivalence  of  the  plane  wave 

expansion  (C.12)  with  the  Rayleigh-Scnraerfeld-Green  expression  for 

diffraction  by  a  planar  apperture  located  in  the  transverse  plane 
.  .  73 
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(C.13) 


thereby  emphasizing  the  non- paraxial  nature  of  the  latter.  The 
integral  (C.13)  can  be  written  in  mere  compact  form  as  a  convolution 


SCtr)  =  £(r£)  *  G(r£) 


(C.14) 


where 


/  ikllr‘lA!\  rz  -Z')  ,  ikilc"2c« 

■  8Z  V  lie'll  /  lic’^c1  \  Tic^)  lie'll  ■ 


(C.15) 


219 


iiIim  ttaUli!'.'*  -j,  itto  i'’1'  oli!|  *;i 


Equivalence  of  the  Rayleigh  theory  with  the  plane  wave  expansion 


demands 


ikz(2ts)(zc-z£)  (zc-zy 
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(C.16) 


If  this  assertion  is  valid  then  the  two  theories  are  equivalent,  aid  we 
may  use  the  plane  wave  expansion  instead  of  the  nen-paraxial  Rayleigh 
theory  whenever  it  proves  convenient. 

To  prove  (C.16)  we  must  show 
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i  yk2-(2ti)2(s^sJ)C2c-z^  i2-(sxx*Syy) 

(C.17) 


ds  ds 
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which  is  equivalent  to  showing 


iki^i 


7 


i  y'  k2-(2ir)2(s2+s2)(2c-zy  i2-(sxx+Syy) 
i  /  kZ-(2Tr)2(sx+s2) 

(C.18) 


The  relation  (C.18)  has  been  demonstrated  by  several  authors 
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APPENDIX  D 


APPARATUS  DETAILS 

1.  ACOUSTIC  ILUMNATICN  SOURCE 

The  source  of  acoustic  illumination  was  a  barius  titanate  disk 
(6.3  ca  dia  x  1.3  ca  thick)  driven  in  a  higher  order  thickness  mode  at 
i  Mi z,  by  a  crystal  controlled  oscillator  and  power  amplifier  deliver¬ 
ing  40  watts  of  r.f.  power.  In  order  to  saxiaize  conversion  of  r.f. 
power  to  acoustic  power  the  disk  was  air-backed,  however,  it  is  doubtful 
that  conversion  efficiency  exceeded  30*.  The  disk  was  stunted  in  an 
altsairags  housing  with  O-rings  to  seal  against  approximately  1/2  ataos- 
phere  pressure.  The  transducer  assesbly  is  shown  in  Figures  D.l  and 
D.2. 

76 

Standard  calculations  indicate  that  the  angular  beam  width 
(first  lobe)  for  a  simple  piston  421  diameter  should  be  approximately 
1°.  In  order  to  increase  this  beat  width  a  horn  was  placed  in  front  of 
the  disk.  The  resulting  bean  width  was  about  10°.  He  bean  pattern  is 
shown  in  Figures  D.3  and  D.4.  The  lack  of  circular  symmetry  is  thought 
to  be  due  to  the  clasping  bolts  which  were  located  in  a  regular 
hexagonal  pattern. 

2.  ACOUSTIC  RECEIVER 

The  device  used  to  detect  the  scattered  acoustic  radiation  was 
a  bariisa  titanate  disk  (.304  ca  dia  *  .203  ca  thick)  operating  in  the 
fundamental  thickness  mode.  Because  of  mechanical  problems  associated 
with  the  ambient  water  pressure  it  was  not  possible  to  a-r-back  the 
disk  although  this  would  have  been  desirable.  The  disk  was  mounted  in 
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3  6061  T6  alissiTEB.  housing  U  6061  T6  alusiiss  boro 
BT  birius  titanate  disk  0  neoprene  "u”  ring 

C  PS  62/U  cable 


Figure  D.l  Illuaiiution  Transducer  Assembly 
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a  slender,  tapered  aluminum  housing  attached  to  the  end  of  a  long  steel 
probe  (2.5  cm  dia  x  460  cm  long).  The  acoustic  sensing  element  was 
coupled  to  an  amplifier  chain  providing  a  voltage  gain  of  approximately 
10*\  Because  of  the  high  gain  and  because  of  the  proximity  of  receiver 
and  transmitter  special  precautions  were  required  in  order  to  eliminate 
stray  r.f.  pickup  in  the  receiver  chain.  Elimination  of  stray  r.f. 
pickup  is  especially  important  in  acoustic  holography  because  such 
spurious  signals  will  be  manifested  as  a  coplanar  reference  signal. 
Accordingly,  the  first  two  stages  of  amplification  were  housed  inside 
the  probe  immediately  adjacent  to  the  sensing  element  and  were  powered 
by  self-contained  batteries.  Thus  the  low  level  stages  were  isolated 
by  at  least  250  cm  of  water.  The  receiver  assembly  is  shown  in 
Figures  D.5  through  B. 8. 

Standard  calculations  indicate  that  the  angular  beam  pattern 
width  (first  lobe)  for  a  simple  piston  2X  diameter  should  be  approxi¬ 
mately  35°.  This  basic  pattern  will  be  modulated  by  the  response 
characteristic  of  the  crystal  lattice  and  by  internal  reflections.  The 
pattern  for  the  receiver  used  is  shewn  in  Figure  D.9. 
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!  D.6.  Receiving  Transducer  Body,  Assembled 

Detector  element, 3  mm  diameter,  located  in  the  relief  between  the  2  and  3  inch  marks 


Receiving  Transducer  Body*  Disassembled 

sft  to  right:  27. S V  battery  pack,  detector  element  assembly,  amplifier 


APPENDIX  E 


CHARACTERISTICS  OF  TOE  PHOTO-ACOUSTIC  CONVERSION 

The  process  by  which  a  photographic  transparency  P(r)  is  formed 
from  an  acoustic  field  A(r)  is  conveniently  separated  into  three  stages: 
1)  the  electro-optical  conversion,  2)  the  optical  to  photographic 
conversion,  and  3)  the  scale  reduction.  See  Figure  V.10.  Dependence 
of  P  upon  A  is  most  easily  controlled  electronically  in  the  electro- 
optical  conversion  stage  rather  than  photographically  in  the  following 
stages. 

The  barium  titanate  receiving  transducer  generates  an  electric 
potential  E  proportional  to  the  incident,  acoustic  field.  Since  E  is  a 
low  level  signal  it  is  amplified  100  by  preamplifier  A  located  inside 
the  probe,  deep  underwater  where  stray  RF  fields  are  small.  The 
amplified  signal  is  combined  with  the  reference  and  then  passed  to  the 
cathode  of  the  CRT  display  through  gated  amplifier  G.  By  appropriate 
choice  of  G  and  electronic  processing  U  a  wide  range  of  characteristics 
P(A)  may  be  obtained. 

A  Tektronix  type  561  oscilloscope  with  P-31  phosphor  was  used  for 
the  CRT  display  because  the  electrostatically  focused  CRT  of  this  unit 
has  a  small  spot  size  (.25  nan  diameter)  which  is  relatively  uniform  in 
size  and  intensity  over  an  8  *  10  an  display  area  (320  *  400  resolution 
cells).  The  P-31  phosphor  produces  a  bright  green  spot  (spectral  peak 
at  .52  micron,  spectral  range  from  .4150  to  .6000  micron)  which  decays 
in  32  msec,  to  0.1%  of  peak  intensity,  long  enough  to  adequately  time 
average  1  MHz  signals  but  fast  enough  to  follow  fringe  variations 
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during  scanning.  In  retrospect,  a  magnetically  focused  CRT  with  P-15 
phosphor  (.05  msec,  decay  time)  would  have  been  of  great  value  since 
these  devices  may  have  resolutions  almost  10  t-mes  greater  than  the 
best  electrostatically  focused  CRT,  particularly  important  when  diffuse 
scattering  is  involved  (holograms  #13  and  #14). 

Typical  characteristics  of  the  electro-optical  conversion  are 
shown  in  Figure  E.l  plotted  as  mormalized  CRT  spot  intensity  I/ctl^  as  a 
function  of  signal  amplitude  input  E  to  amplifier  G;  two  linear  gains 
are  shown.  By  introducing  non-linear  amplification  it  is  possible  to 
alter  the  shape  of  1(E).  For  example,  the  Child'sLaw  response  of  the 
CRT  electron  gun  may  be  compensated  to  give  an  overall  square  law 
response  for  1(E). 

The  CRT  spot  intensity  was  obtained  by  measuring  the  open  circuit 

voltage  of  an  International  Rectifier  type  A5M  selenium  photocell 

placed  in  the  camera  film  plane  with  aperture  set  at  f/1.9.  The  open 
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circuit  voltage  of  this  cell  is  given  by 


-  "  G  +  Voi0) 


(E.l) 


where  Iq  is  the  dark  current  of  the  cell,  a  is  a  conversion  constant 
relating  junction  current  to  incident  light  intensity  I,  T  is  cell 
temperature ,  K  is  Boltzmann's  Constant,  and  e  is  the  electronic  charge. 
Normalized  intensity  I/aIQ  is  sufficient  information  since  proper  film 
exposure  is  most  easily  determined  by  trial-and-error. 


The  scanned  CRT  display  is  time  photographed  with  a  Tektronix 
type  C-12  camera  equipped  with  f/1.9  lens.  Choice  of  film  was  based 


upon  obtaining  high  resolution  and  contrast  of  the  scanned  fringes. 

This  demanded  a  fine  grain  emulsion  (100-150  lines/nm  resolution)  with 

anti-halation  backing  and  a  photographic  speed  approximately  ASA  50.  It 

was  anticipated  that  an  orthochromatic  spectral  response  would  be 

required  in  order  to  prevent  any  cathode  glow  penetrating  the  CRT 

phosphor  layer  from  exposing  the  film  during  the  long  exposure  period 

(4-9  hours) .  This  precaution  proved  unnecessary. 

Two  film  types  were  tried:  Kodak  Contrast  Process  Ortho,  a  high 

contrast  emulsion,  and  Kodak  Panatotnic-X,  a  continuous  tone  emulsion. 

Both  provide  gray  scale  rendition  but  the  relative  dynamic  range  of  the 

high  contrast  emulsion  is  much  smaller.  Polaroid  type  55  P/N  is  also 

suitable  and  is  convenient  tc  use.  but  its  self -developing  feature 

allows  little  lattitude  in  developing. 

Photographic  exposure -development  characteristics  are  usually 

plotted  as  photographic  density  (log^l/Pj  where  Pj  represents 

intensity  transmittance  of  the  developed  emulsion)  against  log^  where 

E  is  the  exposure  (total  energy) .  Such  plots  are  known  as  Hurter- 

Dr if field  curves.  In  this  work  it  is  more  convenient  to  plot  Pj  versus 

E  for  incoherently  illuminated  transparencies  but  P*  (amplitude 

transmittance  of  the  developed  emulsion)  versus  E  for  coherently 

78 

illuminated  transparencies  .  Typical  H-D  curves  for  Contrast  Process 

Ortho  film  processed,  in  a  high  contrast  developed.  Kodak  typ>e  D-ll,  and 

Panatomic-X  film  processed  in  a  continuous  tone  developer,  Kodak  type 
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D-76,  are  shown  in  Figure  E.2  . 

After  processing,  the  transparency  (on  one  of  the  above  emulsions) 
obtained  from  the  CRT  display  is  reduced  in  size  approximately  5  to  1 
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Amplitude  Transmittance  as  a  Function  of  Exposure  for  Kodak  High 
Resolution  Glass  Plates 

D-X9  developer  for  5  minutes  at  68°  F.  (Exposure  scale  x  100) 


onto  Kodak  High  Resolution  Glass  Plates  processed  wit h  £od*k  Type  B-19 
developer.  Typical  versus  F  characteristics  for  this  esalsica  are 
given  in  Figure  E.3^®.  The  final  reduced  transparency ,  the  acasstic 
hologram,  is  a  '‘"positive"  of  "die  sccaistic  fielfi,  that  is,  CTirges  of 
large  acoustic  amplitude  are  recorded  as  rrisges  of  high  iigat 
amplitude  transmittance  in  the  hologram.  The  overall  response 
characteristic  for  the  acoustic  hologram  transparency  is  obtained  by 
canbining  the  three  component  responses  proceeding.  Typical 
characteristics  are  shown  in  Figure  t.4. 


APPENDIX  F 

RECONSTRUCTION  FROM  COPLANAR  REFERENCE  HOLOGRAMS 


We  present  here  the  special  techniques  required  in  order  to 
obtain  quality  reconstruction  from  coplanar  reference  holograms  (often 
encountered  in  acoustic  holography) .  The  basic  problem  with  coplanar 
reference  holograms  is  that  the  zero  order  diffracted  wave  (geo¬ 
metrical  shadow  of  the  hoiogram)  is  superimposed  upon  and  interferes 
with  the  reconstructed  images.  In  order  to  obtain  useful  reconstruct¬ 
ions  this  interference  must  be  eliminated.  (There  is  also  interference 
between  the  two  images  but  this  is  of  relatively  less  importance.) 

The  process  for  eliminating  zero  order  interference  can  be 
understood  by  examining  relations  (TV. 56)  and  (IV. 62)  for  the  distance 
Zj  between  reconstructed  image  planes  and  the  hologram  plane.  These 
equations  are  sunsnarized  as 


!i  ±m2(l  -  Zq/Zj)  .  MZo/zc 


cr-i) 


where  m  is  the  ratio  of  original  acoustic  field  dimension  to 
corresponding  hologram  transparency  dimension,  y  is  the  ratio  of 
acoustic  recording  wavelength  to  visible  reconstruction  wavelength,  Zq 
is  the  distance  between  object  point  and  hologram  plane,  z^^  is  the 
distance  between  reference  focus  and  hologram  plane,  and  zc  is  the 
distance  between  reconstruction  illumination  focus  and  hologram  plane. 

For  the  important  case  z^  ®  (valid  for  all  holograms  described 
in  Chapter  V)  Equation  F.i  reduces  to 
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If  the  demagnification,  the  wavelength  scale,  and  the  distances  z^  and 
Zq  are  such  that 


(F.3) 


then  Zj+  corresponds  to  a  real  image  located  between  z^,  and  the  hologram 
plane,  and  z^  corresponds  to  a  virtual  image  located  behind  the 
hologram  plane.  If  condition  (F.3)  is  not  satisfied  then  both  images 
are  real,  z^+  lying  beyond  zc,  Zj_  lying  between  zc  and  the  hologram 
plane.  The  later  condition  is  illustrated  in  Figure  F.l. 

It  is  now  clear  how  to  eliminate  the  zero  order  field;  simply  use 
convergent  illumination  during  reconstruction.  Then  the  zero  order 
will  cane  to  focus  at  zc  which,  as  we  have  already  seen,  lies  between 
the  two  images.  By  placing  a  minute  opaque  stop  S  at  as  shown  in 
Figure  F.l,  the  zero  order  can  then  be  removed  without  appreciably 
effecting  the  image  fields.  Zero  order  removal  can  also  be  accomplish¬ 
ed  by  placing  the  converging  lens  directly  in  front  of  rather  than 
behind  the  hologram. 

From  Equations  (IV. 61)  and  (IV. 62)  the  lateral  magnification  M, 
the  fraction  by  which  the  image  is  smaller  than  the  object  in  planes 
parallel  to  the  hologram  plane,  is 


ire  F.l.  Viewing  Acoustic  Holograms 
Closed  Circuit  Television 


Under  the  condition  z^  ® 

M±  *  5  ±  '  (F-S) 

Thus,  when  both  images  are  real,  the  one  lying  closest  to  the  hologram 
will  be  smaller  and  reversed  from  the  other.  Equation  (F.5)  shows  that 
by  suitable  adjustment  of  the  parameters  large  magnifications  can  be 
obtained.  However,  as  the  magnification  increases  so  does  Zj  and,  as  a 
matter  of  practicality,  the  images  will  be  restricted  to  rather  small 
sizes  as  a  result  of  the  large  demagnification  m.  Thus,  in  order  to 
view  the  images  magnification  will  be  required,  for  example  by  placing 
a  lens  behind  the  zero  order  stop. 

We  now  briefly  discuss  a  unique  method  for  display  of  the  images 
reconstructed  from  acoustic  holograms.  As  we  have  shown  these  images 
will  usually  be  real  rather  than  virtual  and  will  be  of  small  size. 
Hence,  the  display  of  acoustic  hologram  reconstructions  is  a  problem  in 
real  image  visualization.  Conventionally,  acoustic  holograms  are  view¬ 
ed  by  magnifying  and  then  projecting  the  reconstructed  real  images  onto 
a  diffusely  reflecting  or  transmitting  screen.  However,  there  is 
another  technique  which  is  often  more  useful. 

We  show  in  Figure  F.l  a  scheme  for  viewing  acoustic  holograms  in 
which  the  real  image  is  projected  directly  onto  the  photosensitive 
surface  of  a  vidicon  television  camera  tube  and  the  reconstruction  is 
then  viewed  on  a  TV  monitor.  There  are  several  advantages  to  this 
system  of  viewing: 

1)  the  image  is  magnified  electronically  (typically  between  12 
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and  25  to  one  depending  upon  the  relative  size  of  vidicon  and  monitor), 

2)  the  slight  jitter  inherent  in  the  TV  raster  tends  to  average 
out  the  speckling  due  to  coherence  effects, 

3)  saturation  effects  tend  to  eliminate  the  mottling  character¬ 
istic  of  acoustic  images, 

4)  the  contrast,  brightness,  and  threshold  sensitivity  of  the 
display  may  be  varied  instantaneously  and  continuously  over  a  wide 
lattitude. 

The  last  is  most  important  because  acoustic  reconstructions 
appear  as  bright  images  against  a  uniform  background  and  it  is  often 
difficult  to  distinguish  one  from  the  other  using  conventional 
methods.  With  the  variable  display  it  is  quite  easy  to  separate  the 
bright  image  from  the  dim  background.  Numerous  examples  of  the  utility 
of  this  viewing  technique  are  presented  in  Chapter  V;  for  example, 
compare  Figures  V. 23(a)  and  V. 23(c). 
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conventional  holography.  A  scanning  techniqu:  for  generating  acou.s.i'*.  holograms  of 
underwater  objects  in  the  laboratory  is  described  in  detail.  Using  this  systenj  acousti< 
hologra-ns  have  been  recorded  which  show  angular  resolution  of  3.6  miiliTadians,  approxi¬ 
mately  1.5  times  the  Rayleigh  limit.  A  variable  contrast  television  display  was  used  t< 
view  the  acoustic  holograms.  To  limit  the  attenuation  of  acoustic  wave  in  sea  water  to 
a  tolerable  value,  only  acoustic  signal  with  frequencies  below  1  MHz  should  be  used.  Wi 
used  1  MHz  signal  for  our  experiment.  Consequently,  the  quality  of  any  acoustic  image  : 
degraded  by  poor  resolution  and  specular  reflection.  It  is  suggested  that  diffuse  or 
incoherent  illumination  bt  used  to  overcome  the  defect  that  acoustic  images  often  appeal 
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